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ABSTRACT 


THE FEASIBILITY OF TROPOSPHERIC AND TOTAL OZONE 
DETERMINATION USING A FABRY-PEROT INTERFEROMETER AS A 
SATELLITE-BASED NADIR-VIEWING ATMOSPHERIC SENSOR 

by 

Allen Maurice Larar 


Chairperson: S. Roland Drayson 

Monitoring of the global distribution of tropospheric ozone (O3) is desirable for en- 
hanced scientific understanding as well as to potentially lessen the ill-health impacts 
associated with exposure to elevated concentrations in the lower atmosphere. Such 
a capability can be achieved using a satellite-based device making high spectral res- 
olution measurements with high signal-to-noise ratios; this would enable observation 
in the pressure-broadened wings of strong O3 lines while minimizing the impact of 
undesirable signal contributions associated with, for example, the terrestrial surface, 
interfering species, and clouds. The Fabry-Perot interferometer (FPI) provides high 
spectral resolution and high throughput capabilities that are essential for this mea- 
surement task. Through proper selection of channel spectral regions, the FPI opti- 
mized for tropospheric O3 measurements can simultaneously observe a stratospheric 
component and thus the total O3 column abundance. Decreasing stratospheric O3 
concentrations may lead to an increase in biologically harmful solar ultraviolet radi- 
ation reaching the earth’s surface, which is detrimental to health. 


In this research, a conceptual instrument design to achieve the desired measure- 
ment has been formulated. This involves a double-etalon fixed-gap series configura- 
tion FPI along with an ultra-narrow bandpass filter to achieve single-order operation 
with an overall spectral resolution of approximately .068 cm -1 . A spectral region 
of about 1 cm " 1 wide centered at 1054.73 cm " 1 within the strong 9.6 /xm ozone 
infrared band is sampled with 24 spectral channels. Other design characteristics 
include operation from a nadir-viewing satellite configuration utilizing a 9 inch (di- 
ameter) telescope and achieving horizontal spatial resolution with a 50 km nadir 
footprint. 

A retrieval technique has been implemented and is demonstrated for a tropical 
atmosphere possessing enhanced tropospheric ozone amounts. An error analysis 
assessing the impact on retrieved O 3 amounts of the most significant uncertainties 
associated with this particular measurement has been performed for several different 
types of atmospheres. Results show the proposed instrumentation to enable a good 
measurement of absolute ozone amounts and an even better determination of relative 
changes, with a range of accuracy to within 7.55 to 20.6% for integrated tropospheric 
amounts (and 1.99 to 4.02% for total 0 3 column abundance) and a corresponding 
range in precision to within 7.73 to 10.4% (and 3.30 to 3.95% for total 0 3 column 
abundance), for the atmospheric conditions considered. 


This thesis is dedicated to my parents for setting the foundation for me to make my 

goals such as this degree attainable. 
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CHAPTER I 


INTRODUCTION 


In recent years, there has been an ever increasing need for improvements in the 
capabilities to remotely sound ozone (0 3 ) in the earth’s atmosphere. This deepen- 
ing interest has evolved from advances in the understanding of the significance of 
atmospheric ozone, as well as from recognition that anthropogenic activities may 
potentially affect the ozone spatial distribution and total column abundance. 

Approximately 90 percent of ozone present in the earth’s atmosphere is located in 
the stratosphere with maximum mixing ratios of around 10 parts per million by vol- 
ume (ppmv) near 35 km. Most of the remaining ozone is situated in the troposphere 
where clean, background concentrations can range from 20 to 80 parts per billion 
by volume (ppbv) due to natural variability; however, tropospheric ozone concentra- 
tions can get as high as 500 ppbv in polluted urban areas (Seinfeld, 1986). Despite 
its relatively small abundance, atmospheric ozone plays a vital role in shielding life 
on earth from biologically harmful solar ultraviolet (UV) radiation; hence the term 
“good ozone” which is sometimes used to refer to stratospheric ozone. Conversely, 
tropospheric ozone is often called “bad ozone” since elevated surface concentrations 
can cause detrimental health effects, as well as damage to crops and other vegetation. 

Many chemical species must be considered in order to understand the current-day 
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atmospheric ozone distribution and how it might change in the future. The man- 
made and natural chemicals of interest include the nitrogen oxides ( NO x ) from urban 
surface emissions and high-altitude aircraft, nitrous oxide ( N 2 O ) from agricultural 
and combustion practices, chlorofluorocarbons (CFC’s) used primarily as aerosol pro- 
pellants and refrigerants, bromine compounds, carbon monoxide (CO) and carbon 
dioxide (CO 2 ) from combustion processes, methane ( CH 4 ) and other hydrocarbons 
which come from a variety of natural and anthropogenic sources. It has become 
well known that atmospheric concentrations of these species as well as other anthro- 
pogenically emitted gases have been increasing over recent years. This is of concern 
due to a potential climatic impact since these species in addition to O 3 and water 
vapor (H 2 O) are significant in the energy budget and resulting thermal character 
of the atmosphere, and due to the potential of increased penetration of biologically 
harmful UV fluxes to the earth’s surface resulting from reduced stratospheric ozone. 
Stratospheric ozone depletion occurs primarily through catalytic cycles involving hy- 
drogen ( HO x ), nitrogen ( NO x ), chlorine ( ClO x ), and to a lesser extent bromine 
(. BrO x ) species. 

Pollution produced by urban areas in the Northern Hemisphere is causing in- 
creased tropospheric ozone on a regional and possibly even global scale. Levels of 
tropospheric O 3 have been increasing and will continue to increase as concentrations 
of the precursor gases (oxides of nitrogen, methane and other hydrocarbons) neces- 
sary for the photochemical formation of tropospheric O 3 continue to rise. There is 
evidence suggesting that average surface 0 3 concentrations have more than doubled 
in the past century at middle to high Northern Hemisphere latitudes (Volz and Kley, 
1988). The finding of enhanced tropospheric ozone over the south tropical Atlantic 
Ocean (Fishman et al., 1990), in a region remote from industrialized areas, shows 
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that the problem is not unique to urban areas; this situation appears to result from 
transport of precursors emitted from wide spread vegetation burning primarily in 
Africa, and the subsequent photochemical production of ozone. 

Information on the polar ozone distribution has become increasingly desirable 
during the last several years. The unique meteorological conditions present in polar 
regions enable heterogeneous chemistry, involving ice crystal surfaces within Polar 
Stratospheric Clouds (PSC’s), to enhance the catalytic destruction of stratospheric 
ozone (primarily by man-made chlorine species) leading to massive depletions in 
the polar springtime; this was first realized with the initial reporting of the “ozone 
hole” over Antarctica in 1985 (Farman et ah, 1985; Stolarski et ah, 1986). In addi- 
tion, global total ozone is also decreasing, with rates having a geographical location, 
altitude, and seasonal dependence (Watson et ah, 1988). In fact, there is strong ob- 
servational evidence from many different measurement techniques showing that, over 
the past decade, annual average ozone has decreased by approximately 10 percent in 
the mid- and high-latitude lower stratosphere (WMO, 1991). The interest regarding 
ozone depletion over sparsely-populated polar regions is not simply due to scientific 
curiosity but, more importantly, the concern of whether similar depletions can occur 
in the more densely populated lower latitudes; this might be possible from dynamical 
transport of ozone depleted air from the polar regions (through the vortex or after 
the vortex breaks down during the polar springtime), or possibly from an analogous 
heterogeneous chemical mechanism involving stratospheric aerosol surfaces. 

Knowledge of the global scale distribution of atmospheric ozone and its temporal 
variability can only be achieved through implementation of a satellite-based remote 
sensing technique. Many satellite-based passive instrumentation methods (having 
the radiation source external and thus not provided by the instrument) have been 
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previously implemented primarily to infer stratospheric ozone amounts and, to a 
lesser extent, total column abundances of atmospheric ozone (Grant, 1989). These 
instruments have included absorption techniques involving direct and indirect solar 
radiation as well as emission measurements, and have operated in both limb and 
nadir viewing satellite geometries. For example, the Stratospheric Aerosol and Gas 
Experiment (SAGE) I&II instruments (McCormick et al., 1979; McCormick, 1987) 
are satellite-based spectrometers that measure the absorption of solar UV radiation 
by ozone, and operate in a limb-viewing solar occultation geometry by acquiring 
data only during spacecraft sunrise and sunset events. The Halogen Occultation 
Experiment (HALOE) instrument (Russell et ah, 1993; Baker et ah, 1986) also in- 
fers stratospheric ozone vertical profiles by observing absorption of solar radiation 
from a solar occultation viewing geometry, however, this is accomplished using a 
broadband filter radiometer operating in the IR portion of the spectrum. There has 
been a series of instruments employed to deduce ozone amounts through measure- 
ment of backscattered solar UV radiation observing the terrestrial atmosphere from 
a nadir-viewing geometry. Most recently, the Solar Backscatter Ultraviolet (SBUV 
2) instrument (Frederick et ah, 1986) has been used to infer total ozone as well as 
for ozone profiling. The Total Ozone Mapping Spectrometer (TOMS) instrument 
(Heath et al, 1975) also measures backscattered solar UV radiation similar to an 
SBUV instrument, but it allows for off-nadir viewing by scanning cross-track and it 
provides integrated ozone amounts rather than profiles. Alternatively, atmospheric 
ozone amounts can be deduced from observation of molecular emission. The Infrared 
Interferometer Spectrometer (IRIS) was a Michelson interferometer instrument em- 
ployed in a nadir-viewing configuration to measure IR emission by ozone in the 9.6 
jim band spectral region (Hanel et ah, 1971). Operating in a limb-viewing geometry, 
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the Limb Infrared Monitor of the Stratosphere (LIMS) instrument (Gille and Russell, 
1984) employed a broadband filter radiometer to also measure emission by ozone in 
the IR (9.6 //m band) spectral region. 

Many different types of instrumentation have been designed for the measurement 
of stratospheric and total ozone, however, development of a satellite-based tropo- 
spheric ozone remote sensor has remained a difficult task to accomplish due to the 
challenges associated with such a measurement. Several factors contribute to this 
difficulty. For instance, clouds usually restrict a limb-viewing technique from sensing 
below the upper troposphere, as the long horizontal path is very likely to contain 
cloud interferences in the instrument’s field-of-view (FOV). This effect can be mini- 
mized by the small horizontal footprint (and thus better horizontal resolution) asso- 
ciated with a nadir-viewing configuration, however, then one usually has to contend 
with a large surface contribution to the measured signal. In addition, when view- 
ing in the nadir, the necessary vertical resolution is usually achieved by raising the 
spectral resolution; although, this directly negatively impacts the maximum attain- 
able signal-to-noise of the measurement. Another consideration in sensing the lower 
atmosphere is that the unfavorable contribution of interferant species can be worse, 
as the increased atmospheric pressure causes their overlap effects to be spectrally 
broader or possibly continuum in nature. 

Although a satellite-based instrument capable of directly measuring tropospheric 
ozone has never been developed, data from existing instruments have been used 
to deduce characteristics associated with the global tropospheric ozone field. The 
TOMS instrument is used to measure total column ozone, however, its signal can 
be sensitive to low altitude ozone variations (Fishman and Browell, 1988); therefore, 
its measurements can be used to infer variations in tropospheric ozone (Fishman et 
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al., 1987). In addition, TOMS data have been used in conjunction with vertically 
integrated SAGE profile data to deduce tropospheric “residual” amounts (Fishman 
and Larsen, 1987). This method requires the spatial mapping of coincident data and 
associated measurement errors from two different instruments. 

An instrument capable of observing tropospheric ozone must be of high spectral 
resolution in order to isolate the pressure-broadened wings of strong 0 3 lines and have 
sufficient signal-to-noise to enable removal of effects introduced by the surface and 
interferant species. A Fabry-Perot Interferometer (FPI) has a high optical through- 
put capability while operating at high spectral resolution and has the possibility of 
being very advantageous when compared with conventional spectrometers (Hernan- 
dez, 1986; Jacquinot, 1954; Vaughan, 1989). In addition, a multiplex advantage is 
achievable as all necessary frequencies can be measured simultaneously using a mul- 
tichannel configuration. These characteristics make the FPI a suitable candidate for 
attempting this difficult tropospheric remote sensing measurement. Instruments uti- 
lizing Fabry-Perot interferometry have been used to monitor the earth’s atmosphere 
from the surface for some time and have also been involved in previous satellite- 
based applications (Vaughan, 1989). A spherical Fabry-Perot interferometer was 
used to measure the 630 nm airglow oxygen emission line aboard the Orbiting Geo- 
physical Observatory (OGO- 6 ) during 1969 and 1970 (Blamont and Luton, 1972); 
upper atmospheric temperatures of neutral components in the thermosphere were 
then deduced from these measurements. A plane Fabry-Perot interferometer was 
launched on the Dynamics Explorer (DE) satellite in 1981 to measure the tempera- 
ture, meridional wind and density of several species of excited atoms and ions in the 
thermosphere (Hays et al., 1981). The High Resolution Doppler Imager (HRDI) is 
a triple etalon Fabry-Perot interferometer (Skinner et al., 1987), aboard the Upper 
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Atmosphere Research Satellite (UARS) launched in 1991, designed to infer upper at- 
mospheric wind vectors by observing the Doppler shift of atmospheric absorption and 
emission lines in the O 2 atmospheric bands (Hays, 1982). Also on the UARS plat- 
form is the Cryogenic Limb Array Etalon Spectrometer (CLAES) instrument (Roche 
and Kumer, 1989) which is a spectrometer utilizing solid Fabry-Perot etalons; the 
primary objective of CLAES is to retrieve stratospheric profiles of temperature and 
trace species related to ozone layer photochemistry (including ozone itself) from mea- 
surement of limb IR spectral emissions. Both the DE-FPI and the HRDI instrument 
were designed to perform a multichannel measurement by recording fringes on image 
plane detectors consisting of concentric ring anode elements positioned in their fringe 
planes. These detectors enable simultaneous monitoring of all spectral channels and 
utilize the entire fringe to achieve maximum signal. The satellite-based Fabry-Perot 
instruments thus far which have utilized the entire fringe pattern have been limited 
to the visible portion of the spectrum. Similar applications at longer wavelengths 
have not been exploited due in part to lack of detector technology in the IR spectral 
region. The IR-FPI detection difficulty will be overcome in this research through 
implementation of the new Circle to Line Interferometer Optical (CLIO) system 
(Hays, 1990; U.S. patent #4893003, 1990), developed by researchers at the Space 
Physics Research Laboratory (SPRL) of the University of Michigan. This system 
converts the circular interferometric fringes into a linear pattern which then can be 
detected by conventional linear array detectors; and thus opening the door for a host 
of possible IR-FPI applications. 

The intent of this research is to assess the feasibility of determining tropospheric 
and total ozone using a single-order double-etalon Fabry-Perot interferometer op- 
erating from a satellite-based nadir-viewing geometry. The basic characteristics of 
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ozone will be presented in chapter II, along with the radiative transfer equation for 
upwelling radiation reaching a satellite instrument configured in the nadir-viewing 
geometry. Chapter III and supporting appendix A contain the necessary theory of 
Fabry-Perot interferometry, starting from a simple, ideal instrument and progressing 
to real instrument characteristics in a multi-etalon configuration; chapter III con- 
cludes with a brief discussion of tuning techniques while having a concentration on 
the new CLIO system. A development of the proposed double-etalon FPI instrument 
design to achieve a single-order high-resolution measurement of atmospheric ozone 
is presented in chapter IV. In chapter V, a retrieval technique and error analysis 
procedure are formulated and then applied to realistic atmospheric test cases to sim- 
ulate instrument performance capabilities. A summary of this research as well as 
suggestions for future related efforts are given in the concluding chapter VI. 


CHAPTER II 


SATELLITE-BASED PASSIVE DETECTION 

OF OZONE 


2.1 Introduction 

The only practical means of obtaining knowledge pertaining to the global distribu- 
tion of trace atmospheric species and their associated temporal variability is through 
implementation of a satellite-based remote sounding measurement technique. Ozone 
is the atmospheric constituent of interest for this research endeavor and, accordingly, 
the emphasis of this chapter will be directed toward space-based detection of ozone. 
The chapter will discuss ozone characteristics, including atmospheric spatial distri- 
butions and spectral bands, as well as present the radiative transfer equation for the 
desired nadir-viewing satellite geometry. In concluding, this chapter will formulate 
expressions for measured signal, weighting, and energy contribution functions while 
incorporating the instrumental function to account for the measurement process. 

2.2 Atmospheric Ozone Vertical Distribution Characteris- 
tics 

The atmosphere is conventionally divided up into layers based upon the verti- 
cal structure of its temperature distribution. The approximate positioning of these 
layers, the troposphere, stratosphere, mesosphere, and thermosphere is illustrated 
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Figure 2.1: Atmospheric layers and temperature profile from 0 to 100 km, based 
upon the 1976 U.S. Standard Atmosphere. 
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Figure 2.2: The standard ozone profile from 0 to 100 km, of the 1976 U.S. Standard 
Atmosphere. 
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in Figure 2.1. The average temperature distribution shown is based upon the 1976 
Standard Atmosphere data set (Anderson et al., 1986), and has its vertical character 
resulting from radiative processes involving both solar and terrestrial radiation. The 
rise in temperature in the stratosphere to a peak at the stratopause is due to the 
absorption of solar ultraviolet radiation by ozone. The vertical profile of average 
ozone for the standard atmosphere (U.S. Standard (1976), Anderson et al., 1986) is 
plotted in Figure 2.2 in units of molecular concentration and mass mixing ratio. As 
this Figure illustrates, most of the mass of atmospheric ozone is contained in the 
lower stratosphere with a maximum concentration found at an altitude of around 20 
km; whereas a maximum in ozone mixing ratio occurs much higher, near 35 km. The 
plots shown correspond to average vertical structure; the altitude of maximum ozone 
for an individual profile varies with season, latitude, slightly with longitude, and me- 
teorological conditions. Approximately 90 percent of the ozone mass resides in the 
stratosphere, around 10 percent in the troposphere, and the residual trace amount 
can be found above the stratopause. Ozone extends over a large vertical range. How- 
ever, as even the magnitude of the peak ozone mixing ratio demonstrates, ozone is 
indeed a trace constituent; nevertheless, it is the most significant trace species in 
the middle atmosphere. Its extremely efficient absorption of UV radiation is realized 
when one considers if all the ozone in a column of air extending from the surface to 
the top of the atmosphere were reduced to standard temperature and pressure (0°C, 
1013.25 mb), the column height would be only around 3mm! Total ozone is defined 
as being the total amount of ozone present in a vertical column and is commonly 
measured in Dobson units (DU). A Dobson unit is defined by 1 DU = 10 cm (or 
1 DU = 2.69 x 10 16 molecules 0 3 cm -2 ); therefore, the global mean abundance of 
ozone is approximately 300 DU. The main isotopic species of ozone 16 03 accounts for 
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Figure 2.3: The normal vibrational and rotational modes of the nonlinear triatomic 
ozone molecule. 

over 99 percent of atmospheric ozone. Ozone can also contain an atom of heavy oxy- 
gen 18 0 yielding the minor isotopic species 16 <9 16 0 18 <9 and l6 0 18 0 l6 0 which make 
approximately .4% and .2%, respectively, of atmospheric ozone (Khrgian, 1975). 

2.3 Selection of Spectral Band and Viewing Geometry 

Molecules can have rotational (due to rotations of the nuclei), vibrational (due 
to nuclear vibrations), or electronic (due to motion of the electrons about the nuclei) 
energies, which individually and thus as a sum total can exhibit only discrete values; 
thus they are quantized. In addition a molecule can have translational energy, a 
continuous form of energy due to motion of the molecule’s center of mass in space. 
A loss or gain in energy E = hui occurs when a molecule emits or absorbs a photon 
of frequency o>(s -1 ) 5 respectively. This forms the basis for the existence of line spec- 
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tra, since a molecule can absorb or emit only at discrete frequencies corresponding 
to a change from one quantized energy level to another. Although these molecu- 
lar energies can only change by discrete amounts, numerous energy transitions are 
possible corresponding to different combinations of changes among the electronic, vi- 
brational, and rotational energy levels. However, for a particular wavelength region 
certain transitions may not be possible. For example, electronic energy transitions 
typically correspond to large energy differences occurring at wavelengths less than 
1/xm and are therefore found in the visible and UV spectral regions. Whereas small 
energy differences are associated with pure rotational transitions and thus occur at 
longer wavelengths (20 /im < A < 1000 /im) located in the far-IR portion of the 
spectrum. And the infrared region A < 20 ^m is where one would expect to find 
vibrational transitions. Of interest in explaining the infrared radiative transfer of 
atmospheric ozone are the so called “vibrational-rotational bands.” As its name 
implies, this type of band results from a simultaneous change in rotational energy 
that often occurs when a molecule undergoes a change in vibrational energy level. 
The occurrence of these transitions and thus the corresponding spectral bands can 
be better appreciated by considering molecular structure. Ozone is a nonlinear tri- 
atomic molecule with its atoms situated at vertices of an isosceles triangle (Khrgian, 
1975). It possesses a relatively strong rotational structure, as well as three funda- 
mental vibrational bands, iq,i/ 2 , and 1/3 at 9.066, 14.27, and 9.597 ^m, respectively. 
The normal vibrational and rotational modes and resulting spectral bands for the 
ozone molecule are illustrated in Figure 2.3. The very strong v 3 and fairly strong 
fundamentals combine to form the important 9.6 /im band of ozone. The 1/2 funda- 
mental at 14 /im is well-masked by the intense 15 /xm band of CO 2 . Another strong 
ozone band near about 4.7 /xm is located in a spectral region of low available atmo- 
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Figure 2.4: Schematic illustrating configuration for nadir-viewing satellite geometry. 

spheric radiation, positioned in between the solar and terrestrial emission curves. It 
is also worth noting that ozone’s absorption bands in the UV part of the spectrum 
are the result of electronic transitions of the ozone molecule, as expected, since elec- 
tronic transitions usually occur at shorter wavelengths (A < 1 /xm) in the visible and 
ultraviolet portions of the spectrum. 

The optical properties of ozone, namely its characteristics of absorption and emis- 
sion of radiation, are responsible for shielding the biosphere from potentially harmful 
solar UV energy which results in formulating the stratospheric temperature profile 
shape; however, most importantly from a measurement viewpoint, these optical prop- 
erties provide the basis for remote detection of its presence. Atmospheric sensing 
measurements by satellite instruments utilize two distinct kinds of viewing geome- 
try: the limb view where radiation leaving the atmosphere on a nearly tangential 
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trajectory is observed, and the near-nadir view in which the instrument’s field-of-view 
(FOV) observes radiation leaving the atmosphere in directions almost perpendicular 
to the surface. The limb view generally enables high vertical resolution (determined 
by the instrument’s FOV), however, a coarse horizontal resolution (dictated by con- 
centration distributions and radiative characteristics of the target species as well as 
atmospheric temperature) is usually obtained. In addition, sensing below the up- 
per troposphere is usually not possible due to FOV obstructions by clouds. The 
smaller horizontal footprint of a nadir view (determined by the instrument’s FOV) 
can minimize the limiting effect of clouds in viewing the troposphere. However, the 
vertical resolution in nadir viewing can be somewhat coarse as, like the limb view- 
ing horizontal resolution, it is determined by (distribution, thermal, and spectral) 
properties of the target species. Also, knowledge of surface characteristics (such 
as temperature and emissivity) are usually required to permit removal of the lower 
boundary signal contribution. To enable sensing the troposphere with a relatively 
small horizontally footprint while minimizing the impact of clouds, a nadir-viewing 
geometry (as illustrated in Figure 2.4) will be assumed for this research application. 
When applying this remote sensing technique, one utilizes measurements of radiation 
emerging from the top of the atmosphere in several spectral intervals (or channels) 
ranging from high absorption to low absorption within the band of interest. Ther- 
mal energy observed in the high absorption region originates from high altitudes, 
as emission from lower levels gets absorbed by the atmosphere; conversely, observa- 
tions in lower absorption spectral regions enable viewing lower in the atmosphere. 
However, relatively high stratospheric ozone concentrations further complicate the 
task of sensing tropospheric ozone. Nevertheless, with knowledge of the atmospheric 
temperature profile, measured spectral variation of outgoing radiation can be used 
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Figure 2.5: The ozone 9.6 /im absorption band line strength distribution. 

to infer information on the absorbing gas distribution. 

Ozone possesses absorption bands at various positions in the electromagnetic 
spectrum. While reflected solar UV radiation could be measured, continuous obser- 
vations are possible during day (independent of solar zenith angle) or night by sensing 
in the infrared spectral region. As already discussed, ozone has several vibration- 
rotation bands in the infrared. In particular, some of its stronger bands have absorp- 
tion maxima centered about 4.75, 9.57, and 14.2 /im. As mentioned by Houghton 
et al. (1984), several factors should be considered prior to selecting an absorption 
band for a remote sounding application. The effect of interferant species should be 
minimized by choosing a band that is not overlapped by bands of other atmospheric 
constituents; this would push in favor of the 9.6 fx m band due to the presence of 
strong CO 2 bands that occur near wavelengths of 4.3 and 15 /zm. It is desirable 
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for local thermodynamic equilibrium (LTE) to apply, so that the source function for 
thermal emission from the band will equal the Planck function, and thus simplify the 
radiative transfer calculations. LTE is said to exist for a particular transition when 
the ratio of number of molecules in the excited energy state to those in the ground 
state is given by a Boltzmann distribution at the local kinetic temperature; this will 
be the case if the probabilities of excitation and de-excitation by collisional processes 
are significantly larger than the probabilities of excitation and de-excitation by any 
radiative or chemical processes, and therefore decreases in likelihood at higher alti- 
tudes. From LTE considerations, it is best to view emission in a spectral region near 
terrestrial blackbody emission (8-12 /im), since terrestrial radiation is driving the 
molecules to excitation temperatures near their local atmospheric kinetic tempera- 
tures. Conveniently, ozone has a strong absorption/emission band in this spectral 
window region. It has been estimated (Mlynczak, 1989) that LTE is a valid assump- 
tion for the 9.6 pm ozone band for heights up to slightly lower than 60 km and 70 
km for daytime and nighttime, respectively; since signal contributions at or above 
these altitudes are negligible for a nadir-viewing instrument sensing in the 9.6 pirn 
spectral region, it would be safe to neglect non-LTE effects for such an application. 
It is also beneficial for the wavelength of the chosen band to be long enough such 
that scattered solar radiation is insignificant and thus can be neglected relative to 
thermal emission; this is the case for wavelengths greater than approximately 4 to 
5 pm, which makes the A = 4.75 /rm band a little less desirable. In addition, it is 
preferable to observe a strong signal from the target species of interest. This signal is 
formed by a complex coupling between vertical concentration profiles of the radiating 
species and the vertical temperature field, characteristics of which are determined 
by the particular atmospheric vertical region of interest and associated spatial and 
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temporal variability, a 5 well as the strength of the absorption band. The distribu- 
tions of species and temperature are not controllable, however, some freedom may 
exist in band selection. Actually, from a maximum signal standpoint it is desirous 
for one to choose a band having strong absorption lines located spectrally near the 
peak blackbody emission for the ambient temperature of the atmospheric region to 
be sensed; this roughly corresponds to 10 to 14 //m for the warm to cold temperature 
range, respectively, of the U.S. Standard atmosphere average temperature within the 
troposphere and stratosphere. In the IR portion of the spectrum, the strongest ozone 
bands are near 10 /jm, where absorption results primarily from the 1/3 and v 1 funda- 
mentals of 16 0 3 (Flaud et al, 1990). All of these factors provide a propensity toward 
basing this measurement on the A = 9.6 /xm (wavenumber 1043 cm -1 ) ozone band, 
which is a narrow but strong infrared absorption band consisting of a number of 
closely situated spectral lines. The spectral line distribution of this band (Rothman 
et al., 1987; Flaud et al., 1990) are illustrated in Figure 2.5. 

2.4 The Radiative Transfer Equation 

The radiative transfer equation (RTE), describing the transfer of radiant en- 
ergy through the atmosphere, can be formulated uniquely for the particular satellite 
viewing geometry. For a nonscattering atmosphere in LTE, the RTE expressing clear 
column (cloud free) radiance reaching a nadir-viewing instrument at wavenumber v 
and zenith angle 6 may be written as 

r — oo 3 t[u 0 p) 

1(^6) = e s (u,0)B(u,T s )r(^^Ps) + B(v,T(p)) - dlnp; (2.1) 

e s (p, 6) represents the surface emissivity, T(p) is the atmospheric vertical temperature 
profile as a function of pressure, t(p, 6,p) corresponds to the atmospheric spectral 
transmittance between pressure level p and the instrument, and the subscript s de- 
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notes a surface parameter. In addition, the radiation emitted by a blackbody at 
wavenumber v and temperature T is described by the Planck function B(v,T) and 
is given by 


B(v,T) = 


au 


exp(T)- 1 ’ 


( 2 . 2 ) 


with the constants defined by a = 2 he 2 and 6 = c is the speed of light, while h 
and k are Planck’s and Boltzmann’s constants, respectively. Equation 2.1 assumes 
reflected solar energy to be negligible, since the solar energy reaching the earth’s 
atmosphere (solar Planck function with “dilution” factor) is insignificant relative 
to that of terrestrial origin for the spectral region of application. In addition, the 
contribution from surface reflected downward thermal flux is ignored for this study; 
this assumes the bottom boundary to be well represented by a blackbody, which is 
not a bad assumption since IR surface emissivities are usually close to unity. How- 
ever, note that there are exceptions to the validity of this assumption; for example, 
over desert regions where the surface emissivity may be e, <C 1, for which cases 
inclusion of the surface reflected downward thermal flux may be necessary. The 
first term in the RTE (equation 2.1) represents the contribution from surface emis- 
sion while the second term is due to atmospheric emission. The quantities 
and B(v, T(p )) are frequently referred to as weighting functions and (energy) 
contribution functions, respectively, as they indicate the vertical dependence of the 
information content contained within the upwelling radiance leaving the top of the 
atmosphere. The monochromatic transmittance t(v,0,p) can be given by 


-(v, 0, p) = exp [- J Y, K u ,pids} = exp l~ J K»,dui], 


(2.3) 


where and pi are the monochromatic absorption coefficient and density, respec- 
tively, for species i, ds is an increment in path length, and dui = p{ds is the com- 
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monly used increment of masspath. The summation over the index i accounts for 
the decrease in transmission due to multiple species in the path and arises from 
monochromatic transmittances being multiplicative; for example, the transmittance 
through a layer of N absorbing species (r = T\Ti ■ • • tn) can be expressed as: 


= exp [- J K Ul dui] exp [- J K^du 2 ] • • • exp [- J K^duw] = exp [~ J 

(2.4) 


The transmittance through a series of consecutive homogenous paths can be obtained 
by multiplying the monochromatic transmittance for each path. Transmittances 
through inhomogeneous paths (such as through the earth’s atmosphere) may be 
approximated by a series of thin homogeneous paths (i.e., Weinreb and Neuendorffer, 


1973). 

The physical structure of molecules determines their energy levels associated with 
possible transitions, and thus the central wavenumber about which the molecules ab- 
sorb and emit radiation. Whereas, the applicable broadening mechanism dictates the 
spectral distribution of the particular absorption/emission feature about its central 
wavenumber. Correspondingly, the absorption coefficient of a spectral line can be 
expressed as (Kuhn and London, 1969): 


k v — SP{y — i/q), 


(2.5) 


where S is the line strength (or integrated absorption coefficient) which is defined 
by S = / 0 °° K^dv, vq is the central (or line center) wavenumber, and P(i / — v 0 ) is 
the normalized broadening profile function (/ 0 °° P(v — vo)dv = 1). The shape of the 
profile function is usually determined by two broadening mechanisms for radiation in 
the earth’s atmosphere. Molecular collisions during the absorption/emission process 
give rise to collisional or pressure broadening, and the corresponding profile function 
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Figure 2.6: Approximate altitude dependence from 0 to 60 km of Lorentz, Doppler, 
and Voigt half-widths corresponding to ozone lines within the 9.6 fim 
band for the 1976 U.S. Standard atmosphere. 
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can be represented by the Lorentz line shape as 


Pl{v - i> 0 ) = - 


&L 


(2.6) 


7 r a L 2 + (u - t-o) 2 ’ 

ol is the Lorentz half-width, and represents half the width at half the maximum 
of the profile function. The Lorentz half-width is proportional to pressure and is 
approximately inversely proportional to the square root of temperature. Actually, 
&L °c where n ranges from .5 to 1.0 (Rothman et al., 1987). The other significant 
mechanism for broadening in the terrestrial atmosphere is Doppler broadening which 
originates from a Doppler shift in the frequency of radiation associated with the 
absorption/emission feature due to thermal motion of the radiating molecules. The 
Doppler broadening profile function is of Gaussian form and can be given by 

Pd{v ~ J'o) = —J^exp [- In 2(- — — ) 2 ], (2.7) 

cud V &£> 


where the Doppler half-width can be represented analytically as an = in 

this expression, k is Boltzmann’s constant, T is the temperature in Kelvin, m is the 
mass of the radiating molecule, and c is the speed of light. Unlike the Lorentzian 
half-width, the Doppler half-width does not have a pressure dependence and there- 
fore its change with altitude is due to temperature alone. As discussed by Kuhn 
and London (1969), at low altitudes < 20 km (high pressures) pressure broadening 
dominates and at high altitudes > 50 km (low pressures) Doppler broadening is the 
most significant. However, for intermediate altitudes oti and ac can be comparable 
in size, at which both broadening mechanisms are significant and must be considered 
simultaneously. The resulting broadening profile function is formed by the convolu- 
tion of the Lorentzian and Doppler line shapes yielding the so called Voigt profile 
which can be expressed as (see for example, Armstrong, 1967; Drayson, 1976): 

1 /In 2 y f°° exp(— t 2 ) 


i \ 1 [nZy r 

'(" - Vo) = \ / 

a/) V 7r 7 r j- 


-oo y 2 + (x — t) 2 


dt, 


( 2 . 8 ) 


24 


where the ratio of Lorentz to Doppler widths is defined by y = ^-\/ln 2 and x = 
\/ln 2 corresponds to a wavenumber scale in units of Doppler width. The Voigt 
profile is formed from the combination of two independent broadening formulas, 
and in the high pressure ( y — > oo) or low pressure ( y — ► 0) limits it approaches 
the Lorentz or Doppler profile, respectively. Figure 2.6 illustrates the approximate 
Lorentz, Doppler, and Voigt half-widths as a function of altitude for 0 3 lines within 
the 9.6 /im band; ai 8 = .08 cm" 1 and a d 8 = 2.89 x 10“ 3 cm” 1 (Kuhn and London, 
1969) were used for surface (T = 288 K and p = 1 atm) values of the Lorentz and 
Doppler half-widths, respectively, and n = .76 was used for the Lorentz half-width 
temperature dependence coefficient (Rothman et al., 1987). Note that regardless of 
the dominant broadening mechanism, the far-wing line shape is often represented 
simply by a Lorentzian line shape. The actual far-wing shape may depart from 
Lorentzian, but this will not be of concern for the current application since measure- 
ment channels will not extend too far from line centers in the wings of strong lines 
(i.e., for the tropospheric 0 3 observations). 

2.5 Instrumental Function 

The RTE (equation 2.1) gives an expression to evaluate upwelling atmospheric 
and surface radiance 7(z/, 0), the radiation leaving the top of the atmosphere and 
reaching a satellite sensor, and contains the quantities and 7?(^, T{p)) dr ^^ , 

the weighting and energy contribution functions, respectively, which indicate the 
vertical dependence of the contributions to this radiance. However, 7(//, 0), %^; P ^ 
and B{v , T(p)) dT ^*) do not incorporate effects due to the instrument observational 
process. In a measured signal spectrum S(//, 0), the upwelling radiance within a 
narrow frequency interval £ 7 /, which is given by /(*/, 0)6 v' , gets redistributed by the 
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instrumental function finst( jy — iy )• Therefore, the measured signal can be represented 
by the well known convolution integral 


S(i/,0)= / finst{v - u)I{u\0)du . 
Jo 


(2.9) 


And the limits of integration can usually be replaced by cutoff frequencies associated 
with a bandpass filter. Analogously, the measurement process may be reflected in 
the weighting and energy contribution functions by the formulations 

I? 


W{v,0,p) - 


-v'W 


( 2 . 10 ) 


and 


C(i/,0,p) = 


fo°° finstiv - T(p)) d -^^dr' 


/o°° finst^ - v'W 


( 2 . 11 ) 


respectively. The denominators in equations 2.10 and 2.11 serve to normalize the 
functions with respect to instrument transmission. It is often useful to express the 
measured signal also normalized to instrument transmission, as 


Sn{v,0) 


Sc ?fin,t(v-v)l(v\0)dv 
Jo°° finsti" - v')dv 


( 2 . 12 ) 


to enable direct comparison with atmospheric radiance or signals from other instru- 
ments in the same spectral region. Equations 2.9-2.12 can be used in an iterative 
process varying the instrumental function to establish the optimum instrument type 
as well as its associated parameters for the particular measurement of interest. 


CHAPTER III 


APPLICATION OF FABRY-PEROT 
INTERFEROMETRY 


3.1 Introduction 

The periodic nature of the Fabry-Perot instrument function can be advanta- 
geous when observation of periodic spectra is desired (Houghton, 1961; Wang, 1990). 
However, for most applications, additional optical elements are necessary to reduce 
the effect of unwanted passbands. This is frequently accomplished using additional 
Fabry-Perot etalons (Cocks, et al., 1980; Skinner et al., 1987) in conjunction with 
a bandpass filter. Appendix A contains an overview of Fabry-Perot interferometry 
theoretical concepts fundamental to this research application, including a procedure 
to formulate the associated instrument function including all significant broadening 
mechanisms. This chapter will introduce the concept of a multiple etalon configu- 
ration with particular emphasis placed upon the implementation of a double-etalon 
design. In addition, this chapter will contain a discussion of focal plane fringe detec- 
tion techniques including a detailed summary of the new Circle to Line Interferometer 
Optical (CLIO) system (Hays, 1990; U.S. patent #4893003, 1990) technology. 
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3.2 Spectral Isolation Using Multiple Etalons in Series 

A single order measurement of high spectral resolution is a hard task to achieve 
using an individual Fabry-Perot etalon. Application of a single etalon is restricted 
by the free spectral range uj sr of the interferometer; with the overall effective finesse 
being TV#, a single etalon FPI can observe at most TV# spectral elements of spectral 
width The transmission passbands become narrower as the etalon plate spacing 
is increased (vj S r decreased) or as the plate reflectivity is increased, but this also 
causes the undesirable neighboring passbands to move closer together and decreases 
the instrument throughput, respectively. However, two or more etalons may be used 
together in a series configuration to narrow the passband, and/or extend the effective 
free spectral range, and/or increase the contrast, depending upon the relationships 
between plate spacings and reflectivity coatings of the different etalons. Therefore, 
with the proper selection of etalon parameters, a multiple-etalon configuration can 
eliminate unwanted passbands and improve sideband rejection; and thus enable a 
successful single order, high spectral resolution measurement. 

Application of two etalons in series can be traced back to 1897 (Hernandez, 1986) 
in work done by the pioneers of the FPI instrumentation, Charles Fabry and Alfred 
Perot; however, this was not for spectroscopic work. Use of multiple etalons for 
spectroscopic applications occurred somewhat later in time by other researchers; for 
example, in 1927 Houston (1927) successfully used two etalons in series. And since 
then, multiple etalons have been employed to handle many tasks. A double etalon 
in series configuration has been utilized by many researchers to distinguish closely 
spaced emission lines (i.e., Bradley and Kuhn, 1948; Daehler and Roesler, 1968; 
Daehler, 1970; Cannell and Benedek, 1970). Double etalon systems have been used 
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to record atomic spectral line emissions of 0( l D) at A = 630 nm in the daytime upper 
atmosphere from the earth’s surface (Bens et ah, 1965; Cocks et al., 1980). Triple 
etalon systems have also had many successful applications. The PEPSIOS (Poly- 
Etalon Pressure Scanning Interferometric Optical Spectrometer) device (Mack et al., 
1963) ha 5 been the most popular of the triple etalon configurations; it consists of 
a pressure scanned triple etalon system arranged in a series configuration, having 
nearly equal etalon plate spacings. The PEPSIOS instrument has been applied 
to measurement of the 630nm day glow (Barmore, 1977) as well as for many other 
atmospheric investigations (Hernandez, 1986; Vaughan, 1989). Rees et al. (1982) 
used a triple etalon configuration different from the PEPSIOS design to measure 
atmospheric O 2 absorption lines from a balloon platform. A successful satellite-based 
application of a multiple etalon has been achieved by the High Resolution Doppler 
Imager (HRDI) instrument (Hays, 1982; Skinner et al., 1987); it is a triple etalon 
Fabry-Perot interferometer designed to measure upper atmospheric winds from the 
Upper Atmosphere Research Satellite (UARS). 

Since a double etalon configuration is the simplest representation of a multiple 
etalon system, and since such a configuration is intended for implementation in this 
research, its properties will now be considered in more detail. As mentioned earlier, 
various combinations of etalon parameters may be desired depending on what is ex- 
pected from the multiple etalon system. Plate separation (or spacer) ratios can be 
grouped roughly into three categories (Roesler, 1974): multiple ratios corresponding 
to spacers having quite different lengths (^ > 2), vernier ratios in which the spacers 
have approximately the same length ^ 1 or ^ < 2), or half-multiple ratios with 
= n + where n is an integer. The etalon with the largest plate separation (or 
smallest vj sr ) is usually referred to as the resolving etalon or high-resolution etalon 
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Figure 3.1: Transmission functions corresponding to two different double-etalon con- 
figuration scenarios, a) shows a case of equal plate spacings (-£ = 1). b) 
corresponds to a case of etalon plate spacers having a half-multiple ratio 
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(HRE), as the width of its passbands pretty much determine the multiple etalon 
resolution. Conversely, the etalon with the smallest plate spacing (or largest vj sr ) 
is commonly called the low-resolution etalon (LRE) or suppression etalon, with its 
purpose being to block some of the undesired or parasitic orders of the HRE and thus 
yield an extension of the effective free spectral range. Transmission functions result- 
ing from a couple of different double-etalon configuration scenarios are illustrated 
in Figure 3.1. Figure 3.1a) shows a case of equal plate spacings = 1). Ignoring 
coupling between etalons and broadening effects, the multiple etalon transmission 
function can be given approximately by the product of the individual Airy functions 
(E = A\ • A • • • An); as shown in the Figure, for a double-etalon configuration 
this yields a transmission function of higher resolution and contrast. The broad 
wing characteristic of a single-etalon transmission function is greatly suppressed us- 
ing another etalon in the vernier spacer ratio case, enhancing the purity of spectral 
isolation. However, using equal plate spacings (Figure 3.1a)) will not extend the 
free spectral range through blockage of passbands. A Gaussian function is assumed 
for the shape of the ultra-narrow bandpass filter. As depicted in Figure 3.1a), a 
multi-order transmission function results from the double-etalon/filter system, since 
the bandpass filter width (HWHM) is greater than the effective free spectral range of 
the double etalon. Figure 3.1b) corresponds to a case of etalon plate spacers having 
a half-multiple ratio with ^ = 1.5. Extension of the effective free spectral range 
is achieved for this case, as is obvious from the double-etalon transmission function 
plot. However, due to the finite contrast of etalons, small parasitic (or ghost) peaks 
result where the blocked passbands once existed. Incorporating the bandpass filter 
further reduces the magnitude of the parasitic peaks and nearly eliminates unwanted 
passband orders; thus yielding a single-order transmission function for the double- 
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etalon/filter configuration. It should be noted that various pairs of spacer ratios 
can produce a transmission function with the same extension of free spectral range, 
having ghost peaks of approximately the same strength but in possibly different 
locations. 

The method of coupling Fabry-Perot etalons in series can be quite simple. Etalons 
having the same working aperture diameter may be placed in a series with their axes 
parallel without use of any intermediate optics. The etendue is limited by the etalon 
with the largest plate separation. Alternatively, etalons of different diameters may be 
coupled using magnification between the etalons without loss of etendue (Vaughan, 
1989). As demonstrated by the Airy function (equation A.l), the angular positioning 
(in the fringe plane) of a transmitted wavenumber does not depend upon etalon plate 
separation; this implies that the wavenumber dependence on angle for a multiple- 
etalon system is the same as that for a single etalon, when the multiple-etalon con- 
figuration consists of a simple series without employing intermediate magnification. 
Therefore spectral isolation apertures for a simple series polyetalon can be chosen in 
the same way as for a single etalon. Conversely, etalons of different diameter coupled 
with intermediate magnification each have a different wavenumber dependence on 
angle, since the angle subtended at each etalon by a point on the intermediate aper- 
ture is different for each etalon (Roesler, 1974). Consequently, only a single-channel 
measurement can be made at the central fringe, where the etalons of different plate 
spacing are “in register.” 

As mentioned earlier, the transmission of a simple series of etalons can be given 
by the product of the individual Airy functions, ignoring reflections between etalons 
and broadening effects. Broadening alone can be incorporated by replacing the Airy 
functions with the plate defect broadened Airy functions, and then convolving the 
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Figure 3.2: A simple double-etalon series configuration. The transmission through 
the attenuating medium is defined as r m , while the etalon transmittances 
and reflectances are denoted by the Ei and pi, respectively. 

aperture function L with the resulting product (i.e., E = L*[(Ai *£ ) i)-(«42*£ ) 2)] i for 
a double-etalon configuration ignoring reflections). The inter-etalon reflections can 
be handled in several ways: they can simply be ignored if their effect is shown to be 
insignificant; an attempt can be made to eliminate them by introducing a slight tilt 
of etalon axes with respect to each other, which can interfere with alignment for a 
multichannel configuration, or by using a system of polarizers (Hernandez, 1986) at 
the cost of severely reducing throughput; alternatively, their impact can be reduced 
to acceptable levels by placement of an attenuating medium in between the etalons. 
This last method involving placement of a loss factor in between the etalons can be 
the simplest and most effective approach (Vaughan, 1989). 


33 


FPI double etalon bandpass 



0.60 0.70 0.80 0.90 1.00 

Loss factor of medium between etalons (r m ) 


Figure 3.3: Transmission through a double-etalon system incorporating inter-etalon 
reflections and with an attenuating medium placed in between the 
etalons. Defect broadening effects are not included, a) shows trans- 
mittance through the system including reflections for several loss values 
of the intermediate medium. The solid curve represents the transmis- 
sion with reflections ignored and without any intermediate attenuation 
(r m = 1). b) illustrates the minimum achievable transmittance within 
one order of the HRE central order as a function of loss factor between 
the etalons, with and without the inclusion of inter-etalon reflections. 
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Mack et al. (1963) derived a transmittance function for a triple etalon system 
possessing attenuating media between adjacent etalons while incorporating reflec- 
tions between etalons. Their effort was associated with the PEPSIOS pilot model 
instrument, however, their resulting transmission expression is applicable to other 
triple-etalon systems (Skinner et al., 1987). An analogous expression can be derived 
for a double-etalon system. First one can consider a simple double-etalon series con- 
figuration as illustrated in Figure 3.2. The transmission through the attenuating 
medium is defined as r m , whereas the etalon transmittances and reflectances are de- 
noted by the F, and pi , respectively. The etalon transmittances may be expressed 
by etalon transmission functions (equation A. 17) as F, = A, * Di, and the etalon 
reflectivities can be defined simply by pi = 1 — Ei (Born and Wolf, 1975; Skinner et 
al., 1987; Vaughan, 1989). This simple relationship serves as a very good approxima- 
tion since absorption is small in dielectric reflecting film coatings, leaving scattering 
principally responsible for reduction of intensity; and the scattered ray travels very 
close to the reflected and transmitted rays (Skinner et al., 1987). The expressions 
to the right of the HRE in Figure 3.2 show the contributions to the total intensity 
of transmitted radiation from direct transmission and from the successive reflections 
between the etalons. With an initial intensity of radiation /o entering the LRE, the 
total intensity of radiation exiting the double etalon system is obtained by summing 
the direct transmission and successively reflected components, and can be given by 

I = IoEiE2T m (\ + T m 2 piP2 + T m 4 pi 2 p2 2 + •••)• (3-1) 

Therefore, the system transmittance would be Ep = -f^, or 

OO OG 

E d = FiF 2 r TO [l +'52{ T ,n 2 p\p2)'} = E 1 E2T m J2( T m 2 P 1 P 2 )' • (3.2) 

1=1 i =0 

The Maclaurin expansion of i s the following geometric series: = Y^L o^ n * 
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Utilizing this series expression, the double-etalon transmission function can be ex- 
pressed as 


Ed = 


E x E 2 t„ 


(3.3) 


1 - T m 2 pip2 

Figure 3.3 shows how inter-etalon reflections affect a double-etalon system transmis- 
sion function for different values of attenuation for an intermediate medium r m . As 
r m is reduced from 1.0 to .95, the passband wing suppression is greatly improved 
with the minimum system transmittance dropping by more than one order of mag- 
nitude; thus illustrating how intermediate placement of a low attenuation medium 
can improve the degrading transmission effects introduced by incoherent reflections 
between etalons. 

The effect of a finite sized aperture can be accounted for in a double-etalon 
system similar to its incorporation into a single-etalon system, and analogous to 
equation A. 21 the double-etalon FPI instrumental transmission function can be rep- 
resented by 

T D (V) = E D {V) * L(V). (3.4) 


The construction of a double-etalon FPI instrumental transfer function including 
broadening effects and incoherent reflections between etalons is summarized schemat- 
ically in Figure 3.4. Figure 3.5 qualitatively depicts the formation of a double-etalon 
instrumental transfer function from the Airy (.4), defect (D), aperture (L), and 
bandpass filter ( Tp ) functions. 

It was mentioned earlier (Section A. 3) that for optimum FPI transmittance with 
minimal loss of finesse (or for optimum luminosity-resolution product (LRP)) the 
reflectivity finesse should be matched to the defect finesse of the plates (Nr ~ Ap). 
However, the requirements are more severe for a poly-etalon system (Hernandez, 
1986; Mack et al., 1963; Roesler, 1969, 1974). Plate defect tolerances are smaller 
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Figure 3.4: Schematic representation of the construction of a double-etalon FPI in- 
strumental transfer function including broadening effects and incoherent 
reflections between etalons. The loss factor of an attenuating medium 
(r m ) is included to reduce the impact of inter-etalon reflections. 
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Figure 3.5: Qualitative diagram summarizing the formation of a double-etalon in- 
strumental transfer function from the Airy (.A), defect (D), aperture 
(L), and bandpass filter (7V) functions. The Airy functions simulated 
correspond to an etalon plate spacer ratio of jL = 1.5. All plots shown 
except those for the defect and aperture functions correspond to nor- 
malized transmission functions having the same horizontal and vertical 
scales. However, the defect and aperture functions are plotted with an 
abscissa having one-seventh the range of the abscissa scale used in the 
other plots; this was done to better illustrate their functional forms. 
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for a series of etalons with fixed reflectivity finesse than those acceptable for a single 
etalon system in order to maintain optimum transmittance. Conversely, the plate 
reflectivities in a multiple etalon system with fixed defect finesse must be signifi- 
cantly lower than those for a single etalon to achieve the same peak transmittance 
(Roesler, 1969); and the impact of lower plate reflectivities is lower spectral purity 
(or resolution) for the multi-etalon configuration. A complete treatment for a series 
of etalons would involve specification of each individual defect function and relative 
orientations of the etalons. However, the analysis can be simplified by making some 
basic assumptions while maintaining the essential characteristics of the multi-etalon 
transmission function. This research will assume the same symmetric defect distri- 
bution function to be valid for each etalon ( D\ = = D ('k)), and that correlations 

among defects of the different etalons can be ignored. In addition, the etalons are 
assumed to be aligned perfectly with respect to one another. 

3.3 Methods of Tuning and the CLIO System 

The discussion in Section A. 2 illustrated that the Airy function reaches a local 
maximum in transmission wherever 2fitv cos 9 = n, the integer order of interference. 
It is obvious from this equation that spectral information can be inferred through 
variation of t, or 6. This forms the basis for the three conventional types of 
spectroscopic operation of a Fabry- Perot interferometer. The spectroscopic variable 
can be the index of refraction of the medium between the plates /z (pressure scanning), 
the spacing between the plates t (mechanical scanning), or the angle of incidence of 
the incident radiation 9 (spatial scanning). A schematic representation of these 
spectroscopic variables was given in Figure A.l. 

In pressure scanning, the tuning or scanning is accomplished by varying the 
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pressure within the etalon cavity, thus changing the refractive index of the medium 
between the plates fi. This is usually used in conjunction with a single detector 
centrally located within the fringe plane; thus allowing observation of a single spectral 
“channel” at any given time. 

Mechanical scanning usually has one plate mechanically translated to vary the 
spacing between the plates t. As with pressure scanning, this method typically 
incorporates usage of a single centrally placed detector providing a single-channel 
measurement. 

Scanning spatially in the fringe plane is the foundation for many “spatial scan- 
ning” techniques. This method is sometimes referred to as “angular scanning” since 
the angle 9 changes with spatial position in the fringe plane, or as “tilt scanning” 
since 6 can be changed by simply tilting the etalon. An aperture/detector can be 
scanned across fringes in the image plane; this would yield a single-channel mea- 
surement. Conversely, a photographic plate or imaging photodetector can be fixed 
in the fringe plane; this would constitute a multichannel measurement, since multi- 
ple spectral channels of information can be observed simultaneously. Alternatively, 
many detectors or a diode array may be positioned in the fringe plane to also yield 
a multichannel measurement capability. 

The scanning of multiple Fabry- Perot interferometers is the same as for single 
etalons for pressure scanning and, in the case of a simple series configuration (without 
intermediate magnification), spatial scanning; since the tuning is independent of 
plate spacers. The spatial scanning technique is not applicable to etalons coupled 
with magnification since the angular dependence of transmittance is altered by the 
coupling. Tuning by mechanical scanning is dependent upon plate spacing; therefore, 
any technique to mechanically scan multiple etalons must enable separate adjustable 
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scanning mechanisms for each etalon. 

The intended application of this research is for a satellite-based instrument re- 
motely monitoring atmospheric ozone through utilization of a nadir-viewing geom- 
etry. Accordingly, a multichannel spectral measurement capability is required to 
enable any sort of ozone vertical distribution discrimination while maximizing signal- 
to-noise of the measurements; thus necessitating implementation of a spatially scan- 
ning configuration, from the above discussion of spectroscopic tuning techniques. 
Additionally, a satellite-based instrument further restricts the potential choices of 
tuning techniques (i.e., implementation of a photographic plate would not be feasi- 
ble for a space-based application). Historically, many approaches have been pursued 
to develop multichannel optical devices. Several instrument designs made in the past 
have had the different spectra] channels of information detected simultaneously at 
separate detectors. For example, Hirschberg and Platz (1965) employed a technique 
using concentric mirrors in the fringe plane (each with a slight tilt with respect to the 
optical axis) to reflect light onto separate photodetectors. However, recent advances 
in detector technology have led to achieving simultaneous multichannel operation 
using a single detector unit (Hernandez, 1986; Vaughan, 1989); for example, position 
sensitive devices such as the imaging photodiode (Rees et al., 1981) and more re- 
cently the image plane detector (Hays et ah, 1981; Killeen et ah, 1983), and a charge 
coupled device (CCD) as an imaging device (Abreu and Skinner, 1989). One of the 
most ingenious devices with demonstrated success for space-based applications, the 
image plane detector (IPD) developed for the Dynamics Explorer Fabry-Perot Inter- 
ferometer (DE-FPI) satellite mission provided a multichannel spectral measurement 
through spatially scanning the interference fringe plane using shaped anodes on a 
single detector unit. The DE-IPD is an electro-optical device that converts energy 
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reaching the twelve concentric anodes in the image plane (which mimic the fringe 
pattern of the etalon) into a set of discrete electron pulses. This type of approach 
has some significant advantages over conventional scanning techniques (Killeen et 
al., 1983): the need for mechanical scanning is eliminated, and it provides a signal 
multiplex advantage by detecting all the energy in the imaged ring pattern. In ad- 
dition, a longer dwell time per spectral channel may be achieved since all channels 
are measured simultaneously; this is of particular interest for satellite applications 
where the instrument field-of-view (FOV) is in motion relative to the target scene. 
The successful implementation of an IPD on the DE-FPI led to a similar multiple 
anode detector being developed for the High-Resolution Doppler Imager (HRDI) in- 
strument aboard the Upper Atmosphere Research Satellite (UARS) (Skinner et ah, 
1987; Hays and Wang, 1991). However, this more advanced multichannel anode de- 
vice developed for the HRDI instrument has thirty-two channels instead of the twelve 
used with the DE-FPI. 

Image plane detectors have a proven track record in space-based applications, 
however, lagging detector technology at longer wavelengths prevents their develop- 
ment in the IR region of the spectrum. To circumvent this IR-FPI detection difficulty, 
this research will incorporate the new Circle to Line Interferometer Optical (CLIO) 
system (Hays, 1990; U.S. patent #4893003, 1990) technique to accomplish the spatial 
scanning of the fringe plane, and thus enable multichannel operation. A successful 
laboratory demonstration of this technique has been carried out in the Space Physics 
Research Laboratory (SPRL) at the University of Michigan (Hays, 1990). In addi- 
tion, the CLIO system has been included as a vital ingredient in the Multi-Order 
Etalon Sounder (MOES) instrument design concept (Wang, 1990). 

The CLIO system converts circular fringes of an FPI into a linear pattern which 


42 



Figure 3.6: Conceptual diagram of the CLIO system optical configuration used to 
convert a circular Fabry-Perot fringe pattern into a linear pattern along 
the axis of a 45 degree half-angle internally reflecting transforming cone. 
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Figure 3.7: View from fringe plane of the CLIO system without the circle to line 
converter in place. 
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can be detected by a conventional linear array detector; thus resembling a standard 
spectrometer output, while enabling utilization of the high throughput of an FPL 
In spectra] regions where an IPD can be utilized, application of a CLIO system may 
be more advantageous; the CLIO system can put the same amount of energy on a 
smaller detector surface area thus improving the signal-to-noise, and implementation 
of a CLIO system should be more cost beneficial than the corresponding IPD (Hays, 
1990). The circle to line transformation is accomplished with a 45 degree half- 
angle internally reflecting cone segment of 90 degree azimuthal extent as the primary 
optical element; it is centered above the detector plane while its vertex is centered in 
the fringe plane. In addition, an optical kaleidoscope is located in the optical path 
prior to the cone to enable transformation of the entire circular ring pattern (instead 
of just one- fourth of the ring energy being transformed by the 90 degree cone segment 
alone). The basic optical configuration of the CLIO system is illustrated in Figure 3.6. 
The mirrored kaleidoscope concentrates the entire fringe pattern onto the 90 degree 
cone sector. The cone then reflects the incident energy onto the detector array 
centered on the cone axis below. The overall effect of the kaleidoscope and 45 degree 
half-angle cone combination is to transform energy that would have intersected a 
circle of radius r 0 in the fringe plane onto a detector element whose center intersects 
the cone axis at a distance ro from the cone vertex. The kaleidoscope converts 
all of the energy that would be incident on a complete circle in the fringe plane 
into a small angular segment, corresponding to the azimuthal extent of the cone 
segment. For a 90 degree cone sector, this can be accomplished using a simple 90 
degree “telekaleidoscope” (Hays, 1990): this optical device consists of a segment of 
a telescope objective enclosed by two mirrors joined at the optical axis of the lens 
system and terminated with the transforming cone segment. If viewed from fringe 
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plane without the cone present, it would appear as if there were four identical etalons 
and fringe forming lenses; the real image would be above three virtual images in a 
crossed pattern (see Figure 3.7). Thus implying an etalon “effective” diameter D e jj 
greater than the actual etalon diameter D e . From the simple geometrical relationship 
depicted in Figure 3.7 it can easily be shown that 

D e fj = D e [\ + \/2) « 2.414£) e . (3.5) 

This expression will be of use when design considerations are addressed in the next 
chapter. 


CHAPTER IV 


DOUBLE-ETALON INSTRUMENT DESIGN 
TO ACHIEVE SINGLE-ORDER 
HIGH-RESOLUTION MEASUREMENT OF 
ATMOSPHERIC OZONE 

4.1 Introduction 

The water vapor continuum absorption provides the most significant impact of 
all interfering species within the 10 fim spectral region, since its effect cannot be 
minimized through channel spectral positioning as can be done for spectral lines 
of interfering species. In addition, water vapor is most significant for atmospheres 
within the tropics, geographical regions for which tropospheric ozone measurements 
are highly desirable. Consequently, this chapter will begin with a discussion of the 
water vapor continuum absorption present in the 10 ^m spectral region. A method 
for calculating the upwelling radiance reaching a satellite sensor, including ozone 
and all significant interfering species, will then be presented. Simulations assessing 
the radiance sensitivity to atmospheric ozone perturbations will be shown to suggest 
the necessary instrument spectral resolution and approximate placement of the mea- 
surement band. Double-etalon design specifications to achieve the desired spectral 
isolation will then be developed. Demonstration of ozone information content of the 
measured signals will follow a justification for spectral placement of the instrument 
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band and corresponding channels. The chapter will then conclude with addressing 
additional design considerations leading up to a signal-to-noise ratio analysis. 

4.2 Water Vapor Continuum Absorption in the 10 /tm Spec- 
tral Region 

Absorption by water vapor (HiO) in the 800-1200 cm -1 spectral window region 
of the lower atmosphere is dominated by continuum absorption. An accurate calcu- 
lation of the absorption in this spectral region must include the contribution from 
continuum absorption as well as the local line absorption. However, the contin- 
uum absorption remains poorly understood and difficult to accurately model due 
to the inability of current theories to adequately explain existing observations, the 
large scatter present among laboratory measurements, and the uncertainty associated 
with extrapolating limited laboratory observations to the pressure and temperature 
ranges present in the open atmosphere. Nevertheless, many measurements of wa- 
ter vapor absorption in the thermal IR atmospheric window have been performed 
over the years by a variety of instrumentation techniques (Grant, 1990), from which 
empirical water vapor continuum absorption models can be deduced. 

Water vapor continuum absorption is spectrally unstructured absorption by HiO 
vapor, excluding calculable contributions from the wings of lines within some defined 
proximity. However, different cutoff distances for line wings have been suggested in 
the literature, ranging from 1 cm -1 (Burch and Alt, 1984) to 25 cm -1 (Clough et al., 
1980) with several values in between. Although these undesirable differences in wing 
cutoff definitions may cause some confusion, they are probably insignificant relative 
to the laboratory measurement uncertainties (Grant, 1990). Two main theories have 
been developed to explain the poorly understood water vapor continuum absorption: 
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1) The water vapor absorption line shape model, where the continuum is assumed to 
be caused by extreme wings of absorption lines of H 2 0 molecules (monomers); and 2) 
the water vapor cluster model, where clusters consisting of multiple H 2 0 molecules 
bound together are assumed to result in the observed absorption continuum. The 
water vapor cluster model is also known as the water dimer model since the fraction 
of and corresponding contribution by water dimer dominate over those of all the 
other water clusters. 

Five different theoretical line shapes are conventionally used for collision-broadened 
water vapor absorption lines (Burch, 1981); based upon different sets of assumptions 
regarding the nature of collisions, they all accurately represent H 2 O absorption near 
line centers and degrade to inadequate representation toward the far wings. In fact, 
none of the theoretical line shapes can accurately predict the extreme wing ab- 
sorption without application of an empirically derived line shape correction. Some 
characteristics of the continuum not well explainable by monomer theory, such as 
the pressure square dependence and negative temperature dependence of absorption, 
have supported the water vapor dimer explanation of the observed continuum ab- 
sorption (Suck et al., 1982). However, current H 2 0 dimer theories fail to support 
the water dimer as being the only or even dominant contributor to this continuum 
absorption in the IR region; accordingly, the corresponding models fail to explain 
the spectral character of the absorption. Therefore it is common for one to assume 
that water vapor continuum absorption is caused by absorption from the wings of 
distant strong lines (Thomas and Nordstrom, 1985), which are significant due to col- 
lisions among water vapor molecules (self-broadening) and collisions between H 2 0 
vapor and other species (foreign-broadening); observations show that the absorption 
coefficient increases with the water vapor partial pressure (self-broadening) and with 
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the total pressure (foreign-broadening). Prior to knowledge that different broadening 
is caused by different molecules, the combined broadening effect was considered as 
pressure broadening. However, water vapor self-broadening is now known to usually 
be much stronger than foreign-broadening. 

Considering only water vapor absorption, the absorption coefficient can be ex- 
pressed as 

AC = AC /(i/,r,p) + AC c (l/,J», (4.1) 

where ac/ is the water vapor line absorption coefficient (as discussed in section 2.4) 
which can be accurately computed using the line-by-line method (see section 4.3), 
and ac c is the water vapor continuum absorption coefficient. 

As a commonly implemented first-order approximation, ac c can be expressed as a 
linear function of the water vapor partial pressure and the partial pressure of foreign 
broadening molecules, 

*c(v,T,p) = C,(v, T) P h 2 o + £C/>,7> /i , (4.2) 

t 

where C s and Cj t are the coefficients representing self-broadening and foreign-broadening 
due to the ith species, respectively, in units of cm 2 molecules -1 atm -1 ; pu 2 o and pj x 
are the water vapor partial pressure and partial pressure of the ith broadening com- 
ponent, respectively, in units of atmosphere. Since dry air is composed primarily of 
molecular nitrogen (N 2 ) and molecular oxygen ( O 2 ), the foreign-gas broadening can 
be accurately estimated by considering only N 2 and O 2 . Equation 4.2 can then be 
expressed as 

k c (v,T,p) = C s (v, T)ph 2 o + Cn 2 {v,T)pn 2 + Co 2 {v,T)po 2 , (4.3) 

or by 


«c(^,r,p) = C s (v, T)pn 2 o + C air (v,T)\p- p H2 o]; 


(4.4) 
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p is the total atmospheric pressure, so that p — ph 2 o is the partial pressure of dry 
air; C a ir is the foreign broadening coefficient for dry air, which, upon considering the 
proportions of N 2 and 0 2 in dry air, can be given by 

C aiT = C N2 [ 0.79 + 0.21^1]. (4.5) 

The self- and foreign-broadening coefficients derived by Zhao (1992) for the 500- 
1200 p m region will be utilized for modeling the water vapor continuum absorption in 
this research. Zhao (1992) applied linear regression to the Burch and Alt (1984) data 
using boundary values for good connection to nearby spectral regions in the fitting 
process; in addition, LOWTRAN6 (1020 to 1200 cm -1 ) data were used in deducing 
the foreign-broadening coefficients due to lack of Burch and Alt (1984) data. Zhao 
(1992) arrived upon the following set of expressions: 

C Sw {u) = [15218exp (-0.004216//) - 20.42] x 10" 24 , (4.6) 

and 

C s 296 (//) = [34003 exp (-0.005884//) + 47.56] x 10" 24 , (4.7) 

representing the self-broadening coefficients (in units of molecules -1 cm 2 atm -1 ) for 
284 K and 296 K, respectively, as a function of wavenumber v (cm -1 ); and accounting 
for the self-broadening coefficient temperature dependence, 

C a (u,T) = C S296 (//)[^4^] (23 - 67 - Z ^ ) ; (4.8) 

L/ $284 \ 1 ' ) 

the nitrogen foreign-broadening coefficient wavenumber dependence at a temperature 
of 296 K can be given by 


C Ni96 {v) = [1 .3167 x 10" 6 exp (. 01472//) + 16. 338+1. 6219 xlO 6 exp (— .012383//)] x 10 26 ; 

(4.9) 
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and assuming the temperature dependence of the foreign-broadening coefficient to 
be the same as that for the self-broadening coefficient at a particular wavenumber, 
gives 

Cn(»,T) = ( 4 . 10 ) 

^*284 W 

assuming the ratio of broadening coefficients to be the same as that for broadening 
efficiencies and that the 0 2 to 7V 2 broadening efficiency ratio in the 10 pm region is 
approximately equal to 0.75, equation 4.5 can be used to represent air broadening 
by 

C air (v,T) = Q.%C n (v,T)-, (4.11) 

in addition, the self- and foreign-broadening coefficients are assumed to be constant 
for temperatures lower than 273 K, that is 

C.(v,T) = C s (j/,T = 273K) (4.12) 

and 

C N (v, T) = C N (u, T = 273 K), (4.13) 

when T < 273 K. 

4.3 Atmospheric and Total Upwelling Radiance Calcula- 
tions 


The expression for upwelling radiance given in equation 2.1 of section 2.4 can 
be reformulated into a more direct form which lends itself to higher computational 
accuracy as 


I(v,e) = e,(v,e)B(u,T a )TMp,)+ P B(v,T(p))dT(v,e,p), (4.14) 

J T(u,e,p s ) 


where the integral over transmittance corresponds to along the path of radiation. 
Simulations for this research will be restricted to a nadir (zenith angle 6 — 0°) view- 
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Figure 4.1: Representation of an inhomogeneous atmosphere by a sequence of homo- 
geneous slabs. 
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18 

260.92990 

10.10 

41 

0.20424 

60.50 

19 

237.78906 

10.70 

42 

0.14409 

63.00 

20 

182.18594 

12.40 

43 

0.10126 

65.50 

21 

141.70000 

14.00 

44 

0.06905 

68.10 

22 

110.21021 

15.60 

45 

0.05220 

70.00 

23 

85.76624 

17.20 

46 

0.00032 

100.00 


Table 4.1: Stratification used to represent the 1976 U.S. Standard atmosphere. 
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ing geometry. Since off-nadir views, although usually necessary to achieve the desired 
global coverage using a small horizontal footprint, will not add any further insight 
into this instrument feasibility analysis; this removes the zenith angle dependence of 
the radiative transfer equation (RTE). In addition, the definite integral term present 
in the RTE is usually estimated by a numerical approximation (or quadrature) tech- 
nique. Dropping the zenith angle dependence and utilizing the trapezoidal rule for 
integral approximation, equation 4.14 can be expressed as 

flayers 

I{v)*e,(if)B(i/,T s )T(v,p t )+ J2 £.(",7 1 (p))£'7-«(2',p)- (4.15) 

t=l 

The « sign is shown since implementation of a quadrature technique provides an 
estimate of the definite integral, the accuracy of which depends directly on the num- 
ber of subintervals used in the evaluation; in theory, the quadrature induced error 
— ► 0 as the number of subintervals — ► oo. Radiation trajectories through the non- 
homogeneous (real) atmosphere are usually treated by dividing the atmosphere into 
many vertically successive homogeneous layers (i.e., layers for which the pressure and 
temperature can be represented as constants). The monochromatic transmittance 
through the atmosphere can then be obtained from the product of the transmittances 
for each of the homogeneous layers along the path. The number of layers must be 
> the number of levels for which atmospheric parameter estimation is desired, and 
it must be large enough to achieve sufficient accuracy in the radiance calculations 
(i.e., an iterative process whereby the number of layers is successively increased until 
convergence of the calculated radiance is achieved can be used). In general, finer 
atmospheric stratification produces higher accuracy in the radiance calculations at 
the cost of computational speed. A schematic diagram illustrating the representation 
of an inhomogeneous atmosphere by a sequence of homogeneous slabs is shown in 
Figure 4.1. For the N layer atmosphere depicted in Figure 4.1, equation 4.15 can be 
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expressed as 

I(y) « €,{v)B{v y T a )Ti{y,p s ) + '52B(i/,T i )[T i+ i( v) - t;(i/)], (4.16) 

» = 1 

where Tj is the monochromatic transmittance from level j to the satellite. A 45 layer 
atmosphere was used to represent the real atmosphere in calculations for this study. 
The exact stratification varied slightly depending on the particular atmosphere’s 
thermal characteristics; Table 4.1 shows the layering used for the 1976 U.S. Standard 
Atmosphere (Anderson et ah, 1986) which was used in the majority of calculations. 

The atmospheric parameter retrieval process for remote sounding applications 
usually requires the upwelling radiance to be accurately calculated at high spectral 
resolution. Accordingly, calculations of atmospheric transmittances must also be 
made with high accuracy and at high spectral resolution. This is customarily accom- 
plished using a line-by-line (LBL) method for transmittance calculations (Drayson, 
1966), whereby the contribution to the absorption coefficient from all significant 
spectral lines is accounted for at each frequency of interest. At the expense of com- 
putational speed, the numerical accuracy of the LBL method can be made as high 
as desired, making the uncertainty of spectral line parameters (and their associated 
atmospheric variations) the limiting factor in the overall transmittance calculation 
accuracy. LBL calculations for this research were performed using a slightly modified 
(for adapting to this particular research application) computer program developed 
by Professor S. Roland Drayson at the Space Physics Research Laboratory (SPRL) 
of the University of Michigan; this LBL program originates from earlier work by 
Drayson (1966). The necessary atmospheric parameters (pressures p t -, temperatures 
Ti, and masspaths u t ) required to uniquely describe the vertical sequence of hori- 
zontally homogeneous slabs were calculated using a modified version of the air mass 
computer program FSCATM (Gallery et ah, 1983) contained within the Fast At- 
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Figure 4.2: Spectral lines located within the 9.6 fim ozone band according to the 
1992 AFGL compilation. 
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Figure 4.3: Atmospheric, transmission and upwelling radiance computed at .01 cm -1 
resolution for the 1976 Standard Atmosphere. Calculations include spec- 
tral lines for all species contained within 1992 AFGL line parameter com- 
pilation except OOF? and H 2 S. 
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mospheric Signature Code (FASCODE). The 1992 Air Force Geophysics Laboratory 
(AFGL) atmospheric absorption line parameters compilation (Rothman, 1993) was 
utilized; it is based upon earlier releases (Rothman et al., 1987) and contains im- 
proved O 3 parameters for the 10 fim spectral region (Flaud et al., 1990). Absorption 
lines in the 10 fim spectral region are shown in Figure 4.2. Considering both line 
strengths and atmospheric concentrations, the most radiatively active molecules in 
this spectral region consist of: O3, H*0, CO*, N*0, and C H 4 . Since vertical profiles 
of SO*, NO*, NH 3 , H*0 2 , and PH 3 were available in the 1976 Standard Atmo- 
sphere database (Anderson et al., 1986) used for this study, these species were also 
included in all radiance calculations. However, profiles of COF* and H*S were not 
readily available and these two species were not included in the calculations since 
they are relatively insignificant in this spectral region. In addition, the radiance 
computations assume the surface to radiate as a perfect blackbody with an emissiv- 
ity e s = 1. Subsequent calculations of transmittances and corresponding radiances 
for the 1976 U.S. Standard Atmosphere showed a quite favorable comparison be- 
tween the programs developed for this research and the FASCODE2 program unit 
alone. Figure 4.3 shows the atmospheric transmission and total upwelling radiance 
(including contributions from surface and atmosphere) calculated for the 1976 U.S. 
Standard Atmosphere within the 10 fim spectral region. 

Through perturbing the vertical distribution of atmospheric ozone in regions 
about which information content is desired and calculating the corresponding change 
in upwelling radiance, one can determine spectra] regions desirable for instrument 
band placement as well as the approximate instrument spectral resolution necessary 
to observe the calculated radiance residuals. To this end, radiance calculations were 
performed for two additional atmospheric scenarios: First a case illustrating sensi- 
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Figure 4.4: Effect of atmospheric ozone perturbations on upwelling radiance within 
the 9.6 /im band region. 
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tivity to a tropospheric ozone perturbation was simulated with a 10% increase in tro- 
pospheric ozone from the nominal (1976 U.S. Standard Atmosphere) vertical profile; 
analogously, calculations with a 1% increase in stratospheric ozone were performed 
to demonstrate sensitivity to stratospheric ozone changes. The ozone perturbation 
magnitudes were chosen to roughly correspond to a 1% change in total column ozone, 
an occurrence of which would be desirable to detect. Although the vertical unifor- 
mity and magnitude of these perturbations may deviate from reality, thus making 
the utility of the radiance residual magnitudes suspect, useful information pertain- 
ing to instrument spectral band placement and resolution can still be obtained. This 
simulation is summarized in Figure 4.4 with plots of the radiance residuals resulting 
from the above described ozone perturbations. As expected, maximum sensitivity 
to tropospheric changes can be found in the wings of the strongest ozone lines, for 
example, the 1048-1058 cm -1 spectral region. The negative radiance residuals result 
from the increased absorption reducing the large surface contribution to a greater 
extent than the radiance supplement from increased emission by atmospheric layers 
cooler than the surface temperature. The positive radiance residuals in the strato- 
spheric perturbation case are due to the dominance of the change in atmospheric 
emission in spectral regions centered about strong lines which therefore minimize 
the surface contribution; since the atmosphere is essentially opaque in these spectral 
regions, tropospheric perturbations cannot be seen and thus yield radiance residu- 
als « 0 in these same spectral locations. Additionally, the sensitivity to changes 
in atmospheric ozone goes to zero as one approaches the wings of the band where 
the absorption lines are less intense. From Figure 4.4, one can conclude that some- 
where in the 1048-1058 c?n -1 spectral region would be ideally suited for placement 
of the instrument band (or channels) to maximize information content for both the 
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tropospheric and stratospheric ozone measurements. 

The character of the radiance residual plots can be better appreciated by expand- 
ing a small spectral region along with the corresponding ozone spectral lines. This 
is shown in Figure 4.5 for a 2 cm” 1 interval near the center of the 1048-1058 cm” 1 
spectral region of interest; transmission through the nominal standard atmosphere is 
also illustrated for a reference of atmospheric opacity as well a s an indicator of rela- 
tive importance of the surface contribution to the upwelling radiance. As Figure 4.5 
demonstrates, regions in the wings of strong absorption lines have the largest sensi- 
tivity to changes in tropospheric ozone and would thus be the best spectral locations 
for placement of tropospheric channels; for example, 1051.8 cm” 1 , or 1052.1 cm” 1 for 
a higher tropospheric sensitivity along with a higher stratospheric sensitivity con- 
tribution. These plots show that an instrument resolution ~ .07 cm” 1 should be 
sufficiently narrow for the tropospheric ozone measurement in addition to enabling 
determination of stratospheric ozone content. 

4.4 Selection of Double-Etalon Instrument Characteristics 

As discussed in appendix A, one is desirous of obtaining the best quality etalon 
plates possible to maximize the overall defect finesse Nd- Many manufacturers (i.e., 
Imperial College Optical Systems (ICOS), LTD) are currently able to produce large 
etalons with working apertures up to 150 mm in diameter of very high plate qual- 
ity with Nd % 50 for the IR spectral region. However, the overall finesse may be 
severely degraded during the mounting process in preparation for space flight (Hays, 
private communication, 1993). An overall defect finesse of 25% below the nominal 
value given above, or Nd ^ 37.5, will be arbitrarily assumed in this research. This 
author feels that if further degradation of Nd is anticipated in preparation for space 


Transmission (%) 


62 






Frequency (cm' 1 ) 


Figure 4.5: Effect of atmospheric ozone perturbations on upwelling radiance in a 
narrow spectral region of strong ozone lines. The transmission through 
an unperturbed atmosphere and the ozone spectral line distribution are 
also illustrated. 
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flight, then an alternative optical component integration procedure should be devised 
and/or the application of solid etalons (in place of the presently assumed fixed-gap 
etalons) might be of consideration for this particular research application. Due to 
lack of apriori knowledge regarding the finesses corresponding to individual defects, 
it will be assumed that random surface imperfections, spherical curvature, and lack 
of parallelism all contribute to the total defect of the plates with equal defect fi- 
nesses. Using equation A. 18, the individual defect finesses can be calculated to be 
n Dg = N Ds = N Dp 65.0. In addition, the defect distribution functions (see, for 
example, Section A. 3) will be assumed to be symmetrical about their means and 
to have mean values of zero (//£> = 0). To maximize throughput of a single etalon, 
the reflectivity finesse should be chosen such that Nr < Nd, which (for Np & 37.5) 
implies plate reflectivities R < .92. However, the peak transmittance in a multiple 
etalon system is lower than that of a single-etalon instrument, unless the plate re- 
flectivities are lower in the multi-etalon configuration. In addition, the multi-etalon 
transmittance decreases sharply and becomes more sensitive to plate defects for high 
plate reflectivities. Roesler (1969) discussed the effect of plate defects on the trans- 
mission function for a poly-etalon FPI. For a double-etalon series configuration, he 
suggested a reflectivity finesse equal to the defect finesse lowered by a factor of .644 
(i.e., Nr = .644A(d) in order for the system transmittance to not be degraded from 
the multi-etalon configuration; for the design under consideration, this implies plate 
reflectivities R & .88. However, the use of lower plate reflectivities can decrease the 
spectral resolution and increase the amount of parasitic or off-band light transmitted. 

As mentioned in section 4.3, the desired spectral resolution is 6 r e ~ .07 cm" 1 . 
Since the double-etalon spectral resolution is determined primarily by the high- 
resolution etalon (HRE), we can use equation A. 24 to estimate the necessary free 
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spectral range for the HRE; the effective finesse can be computed from equation A. 25 
by assuming the aperture finesse to be equal to the reflectivity finesse (i.e., Ni % Nr). 
If we consider, for example, the coating reflectivities for both etalons to be R = .90 
(Nr ~ 29.8), this would imply an effective finesse Nr ~ 18.4 and a free spectral 
range of vj sr « 1.29 for the HRE. The uj sr of the low-resolution etalon (LRE) must 
be chosen such that single-order passband transmission is achieved in addition to 
minimizing the amount of unwanted energy through the optical system. This un- 
wanted energy is frequently referred to as white or parasitic light and can be defined 
by the ratio of all energy transmitted through the system outside the HRE central 
order to the energy transmitted within the HRE central order (Skinner et al., 1987). 
The white light W can then be expressed mathematically as 


W = 


HREv f , r 

fo°° T F (u)T D (u)du - r + HRE > T F (u)T D (v)dv 

t/n ~ J — 


, HRE» fsr 

r **£,„ T F {v)T D {v)dv 

VQ 2~ 


(4.17) 


where Tp(v) is the bandpass filter transmission function, and Td(j/) is the double- 
etalon instrument transfer function which includes broadening effects and reflections 
between etalons (equation 3.4). As illustrated in Figure 3.3, placement of a loss 
medium in between the etalons can significantly reduce the undesirable increased 
white light effect from inter-etalon reflections at the cost of slightly decreasing peak 
transmission; accordingly, a loss medium transmission value of r m = .95 will be as- 
sumed in this research. To optimize the instrument’s spectral isolation capability 
and thus minimize the amount of white light, it is desired that the bandpass filter 
be as narrow as technologically feasible; this also implies a reduction in peak trans- 
mission. For this spectral region, ultra-narrow bandpass filters with a full-width at 
half-maximum on the order of .4% of the central wavenumber (FWHM ~ A%vo) 
have been produced (i.e., by the Optical Coatings Laboratory, Inc., (OCLI) for the 
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Figure 4.6: White light transmitted through a double-etalon Fabry-Perot system. An 
ultra-narrow bandpass filter of Gaussian shape is included. A medium 
with transmission equal to .95 is assumed to be in between the etalons 
to reduce the effect of incoherent reflections. 

CLAES instrument aboard the UARS platform); the corresponding peak transmis- 
sion is T Frnax « 65%. An ultra-narrow bandpass filter with FWHM= .5%^o and 
Tfmax = 65% be assumed in this study; such a filter can be modeled by a Gaus- 
sian shape. Thus the filter transmission will be represented by 

Tf{ v) = exp [ V 0) ( 4 - 18 ) 

where the standard deviation a of a Gaussian can be given by a = • 

Evaluating equation 4.17 using the instrument specifications noted above and 
a range of LRE free spectral range values yields white light amounts depicted in 
Figure 4.6. Although near-vernier gap ratios enable the highest spectral resolution, 
they also correspond to high white light transmission since the HRE parasitic peaks 
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R 

7/W 

w 

YV mxn 

%(cm _1 ) 

"/snlcm" 1 ) 

i > 2 (cm'‘) 

0.850 

0.783 

0.329 

0.067 

0.898 

3.223 

0.860 

0.758 

0.260 

0.068 

0.959 

3.441 

0.870 

0.735 

0.201 

0.068 

1.027 

3.685 

0.880 

0.702 

0.152 

0.068 

1.103 

3.958 

0.890 

0.668 

0.111 

0.068 

1.189 

4.267 

0.900 

0.634 

0.080 

0.066 

1.286 

4.615 

0.910 

0.580 

0.055 

0.068 

1.396 

5.009 

0.920 

0.536 

0.037 

0.068 

1.520 

5.455 

0.930 

0.471 

0.024 

0.067 

1.660 

5.957 

0.940 

0.405 

0.015 

0.066 

1.816 

6.515 

0.950 

0.329 

0.010 

0.066 

1.985 

7.123 


Table 4.2: Simulation demonstrating effect of varying plate reflectivity while main- 
taining fixed spectral resolution on peak double etalon transmission To max 
and minimum white light W m i n through the system. 
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Instrument Design Specifications 

Ultra Narrow Bandpass Filter: 

Gaussian shape assumed; central wavenumber, uq= 1054.73 cm -1 ; 
FWHM = .5 % i/q = 5.274 cm -1 ; maximum transmittance, Tp maI = 65% 


Etalons: 



HRE (etalon #1) 

LRE (etalon #2) 

gap t 

0.3289 cm 

0.09157 cm 

Vfsr 

1.52 cm -1 

5.46 cm -1 

n d 

37.53 

>37.53 

reflectivity R 

0.92 

0.92 

Nr 

37.66 

37.66 

Tetalorimax 

0.732 

0.732 


Additional parameters: 

Aperture finesse Ni = 37.66 

Transmission through medium between etalons T m — .95 

Transmission through remaining optical elements r 0 = .85 

Overall effective system finesse Np ~ 22.35 

Effective resolution Sp E = 0.068 cm -1 

White light W= .037 

Instrument peak transmittance r,= 0.28 

Etalon plates diameter D e ss 13.3 cm 

Telescope diameter dj = 9.0 inches 

Nadir footprint size d^p = 50.0 km 


Table 4.3: Double-etalon instrument design specifications. 
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Change From HRE Central Order (An) 


Figure 4.7: Illustration of ultra-narrow bandpass filter function, broadened double- 
etalon transmission function, and overall instrumental transfer function 
of the central wavenumber for an FPI with the chosen specifications. 
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are not blocked. The effect of the bandpass filter can be seen in Figure 4.6 as 
the local maxima decrease in intensity for increasing integer spacer ratios. In this 
simulation, the white light minimum (W m i n « .08) occurs for a free spectral range 
ratio « 3.59. A design compromise exists between peak throughput and white 


light; since white light can be reduced significantly at the expense of a relatively 
small drop in throughput, if the plate reflectivities are increased while maintaining 
the overall resolution approximately constant (and thus allowing the free spectral 
ranges to get larger). This is demonstrated in Table 4.2. Based on these results, 
a plate reflectivity of R = .92 (Nr « 37.7), with free spectral ranges j// 5ri = 1.52 
cm -1 and Vf ST2 = 5.46 cm -1 will be implemented into this design and subsequent 
instrument simulations. Table 4.3 summarizes the instrument design specifications 
discussed above, with an additional loss factor of r Q = .85 incorporated to account for 
decrease in transmission from the remaining optical elements not yet considered (i.e., 
beamsplitters, lenses, etc.). Note that the choice of central wavenumber vq will be 
justified in the next section, while diameters of the etalon plates, entrance telescope, 
and nadir footprint will be arrived upon in section 4.6. The spectral character of 
the resulting instrument function is illustrated in Figure 4.7; the product of the 
bandpass filter and double-etalon transmission functions (TrTd) shows the effective 
central wavenumber instrument transfer function. 


4.5 Band Placement and Choice of Spectral Channels 


The 1048-1058 cm" 1 spectral region was selected in section 4.3 as having po- 
tential for placement of the instrument measurement band in order to enable the 
desired atmospheric ozone observation. Spectral locations for the individual chan- 
nels will be determined in this section based upon simulations performed using the 
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instrument function designed in section 4.4. However, some of the signal related equa- 
tions introduced in section 2.5 must first be re-expressed to reflect the nadir-viewing 
geometry and the Fabry-Perot instrument function formulation. Accordingly, the 
measured signal S(v) (equation 2.9), the measured signal normalized to instrument 
transmission Sn(u) (equation 2.12), and the energy contribution functions C(i/,p) 
(equation 2.11), can be expressed as 


S(u) = f 9 T d {v - v)TF(v)l{v)dv , 

Jviow 

SvM = -^TTTTT-. < and 


C(v,p) = 


cr t d (v - smw 

QT T - u)T F {u)B{v\T{p)) d -^dv 
r^ h T D {v-v')T F {u')du‘ 


(4.19) 

(4.20) 

(4.21) 


respectively; the ui ow and Vhigh represent the low and high wavenumber cutoffs, re- 
spectively, corresponding to the bandpass filter transmission. To ensure the existence 
of necessary information content, ozone energy contributions functions (O3 ECFs) 
with their wavenumber and pressure dependence must be closely examined. The 
spectral lines of interfering species and their corresponding effect on the O 3 ECFs 
must also be considered. In addition, the signal sensitivity to perturbations in atmo- 
spheric ozone amount needs to be incorporated into the band and channel selection 
process. 

Figure 4.8 shows all spectral lines located within the 1048-1058 cm -1 region. For 
this same spectral region, the O 3 ECFs and how they are affected from including 
interfering species are illustrated in Figure 4.9; the effect on the O3 ECFs from 
including interfering species can be calculated from 


A[C(/',p)lo 3 ] = 


(C(y,p) \aii species ~C(^P) \o 1 only) 


C{v,p)\ 


(4.22) 


O3 only 
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Frequency (cm ) 

Figure 4.8: Spectral lines located within a portion of the 9.6 //m ozone band accord- 
ing to the 1992 AFGL compilation. 
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Figure 4.9: 03 energy contribution functions and how they are affected from includ- 
ing interfering species for the 1048-1058 cm ~ 1 region. 
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Figure 4.10: Signal sensitivity to a 10 percent increase in tropospheric ozone. 
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Figure 4.11: Signal sensitivity to a 1 percent increase in stratospheric ozone. 
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The O3 ECFs in Figure 4.9 show that spectral regions can be chosen to maximize 
stratospheric or tropospheric information content in the measured signal; the effect 
of interfering species’ spectral lines and the water vapor continuum can easily be 
observed in the corresponding A[C(^,p)| 03 ] contour/image color display product. 
Applying equation 4.20 to the upwelling radiance data sets generated during the 
atmospheric ozone perturbation simulations discussed in section 4.3 (see, for exam- 
ple, Figure 4.4), one can calculate the corresponding instrument signal sensitivity to 
changes in atmospheric ozone content. For the 1048-1058 cm -1 spectral region being 
considered, the signal sensitivities to a 10% increase in tropospheric ozone and to a 
1% increase in stratospheric ozone are illustrated in Figures 4.10 and 4.11, respec- 
tively. As the plots demonstrate, convolution with the instrument function serves 
to smooth out the radiance response to the ozone perturbations. It is also inter- 
esting to note that peak sensitivities to the different perturbations occur spectrally 
near each other, due to the mutual close proximity to strong absorption lines; this 
implies that both the tropospheric and stratospheric measurements can probably be 
achieved with the instrument band covering a narrow spectral region. 

The information contained within Figures 4.9, 4.10, and 4.11 was scanned 

thoroughly to determine regions suitable for instrument band placement. Decision 
criteria were established to identify spectral regions possessing high sensitivity to 
ozone variation (especially in the troposphere), while minimizing the impact of in- 
terfering species and maximizing the vertical extent of the O3 information content. 
In addition, based on instrument design considerations, it was desired to find these 
features in one narrow spectral region approximately 1 cm -1 wide. Although several 
regions had potential, the 1054.23-1055.23 cm -1 spectral region meets the desired 
instrument band criteria and was determined to be the most suitable from the can- 
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Upwellinq Radiance 



Figure 4.12: Upwelling radiance reaching and received by the instrument within the 
selected band region. 
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Energy Contribution Functions (ECFs.in mw/m 2 — sr— cm ’) 
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Figure 4.13: 03 energy contribution functions and how they are affected from in- 
cluding interfering species for the instrument band spectral region. 
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Figure 4.14: Effect of atmospheric ozone perturbations on the radiance and signal 
within the instrument band spectral region. The transmission through 
an unperturbed atmosphere and the ozone spectral line distribution are 
also illustrated. 
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Figure 4.15: Signal sensitivities to a 10 percent increase in tropospheric ozone and to 
a 1 percent increase in stratospheric ozone are shown with the selected 
instrument channels. 

didate regions. Figure 4.12 shows the upwelling radiance in the selected band region; 
the resultant smoothing effect of the instrument function (due to the measurement 
process) is also illustrated. Analogous to Figure 4.9, Figure 4.13 depicts the 0 3 
ECFs and the corresponding effect from including interfering species within the in- 
strument band spectral region. Additionally, Figure 4.14 is similar to Figure 4.5, 
however, it represents the instrument band region and includes the signal sensitivity 
to 0 3 changes superimposed over that for the upwelling radiance. 

The next step is to choose channel locations. Since the primary interferant effect 
within the band region is due to the spectrally near-uniform ff 2 0 vapor continuum 
(see Figure 4.13), there are no portions that must be avoided and it is desirable to 
sense the whole region. The channels are defined by the spatial distribution of detec- 
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Normalized 0 3 energy contribution functions 








Figure 4.16: Energy Contribution Functions (ECFs) for the selected channels nor- 
malized to their peak values for an ozone only atmosphere. 
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channel 

"Oc 

p(mb) 

[O^EC F]\ max 


(cm' 1 ) 

of peak layer 

(mw m~ 2 sr~ 1 /cm -1 ) 

i 

1054.240 

33.951 

3.523 

2 

1054.295 

33.951 

4.049 

3 

1054.336 

33.951 

3.985 

4 

1054.377 

983.597 

3.634 

5 

1054.418 

983.597 

5.703 

6 

1054.459 

983.597 

6.058 

7 

1054.513 

983.597 

4.845 

8 

1054.554 

25.416 

3.335 

9 

1054.595 

983.597 

4.908 

10 

1054.636 

983.597 

6.465 

11 

1054.677 

983.597 

5.812 

12 

1054.718 

25.416 

3.322 

13 

1054.759 

25.416 

4.134 

14 

1054.800 

983.597 

3.767 

15 

1054.841 

983.597 

5.135 

16 

. 1054.882 

25.416 

3.925 

17 

1054.923 

25.416 

4.880 

18 

1054.964 

25.416 

5.085 

19 

1055.016 

5.959 

4.116 

20 

1055.057 

983.597 

4.146 

21 

1055.098 

983.597 

5.119 

22 

1055.139 

983.597 

5.265 

23 

1055.180 

983.597 

4.432 

24 

1055.221 

33.951 

3.400 


Table 4.4: Characteristics of the 0 3 ECFs for the instrument channels. 
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Figure 4.17: Effect of including interfering species on ECFs from selected channels 
computed for an ozone only atmosphere. 
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Figure 4.18: Effect on O 3 ECFs within the measurement band from including inter- 
fering species. Note that the channel locations attempt to minimize im- 
pact from the CO 2 (~ 1053.92 and 1055.63 cm -1 ) and H^O (ss 1055.52 
cm -1 ) spectral lines and their associated wing effects. 
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tor elements within the linear detector array (see section 4.6). Since these detector 
elements have a finite size, detection design considerations impose a fundamental 
limit to the total number of channels possible within this region. Approximate chan- 
nel locations of interest were defined such that the measurements would provide 
knowledge from regions of peak sensitivity to O 3 variations, maximum vertical dis- 
crimination capability by having a range of O 3 ECFs, and the ability to reproduce 
the signal spectrum as a continuous function within the band region. The design 
specifications of the detector array were then varied until a “best fit” of the desired 
channel locations was achieved. Figure 4.15 shows the spectral positions for the 
selected set of 24 channels within the plots depicting the signal sensitivities to atmo- 
spheric O 3 perturbations. The corresponding vertical distributions of O 3 information 
content are illustrated in Figure 4.16 with plots of O 3 ECFs for all channels; note 
that the O 3 ECFs have been normalized to their peak values to enable inter-channel 
comparisons. These contribution functions demonstrate the ability to de-couple tro- 
pospheric O 3 signal contributions from those of the stratosphere using the chosen 
channels. Based upon the broadness of these ECFs, ozone profiling capability does 
not appear too likely. However, some vertical discrimination ability should be attain- 
able for the stratosphere along with integrated amounts for tropospheric O 3 . As will 
be discussed and demonstrated in Section 5.5, the vertical character of a retrieved 
O 3 profile for this measurement can be similar to that of the true O 3 profile under 
certain situations. Table 4.4 contains peak values for the O 3 ECFs as well as other 
channel specific information. As mentioned earlier, the H^O vapor continuum pro- 
vides the most significant impact from the interfering species within the instrument 
band region. Accordingly, interfering species primarily effect the O 3 ECFs in the 
tropospheric region; this is demonstrated in Figure 4.17 which emphasizes the lower 
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atmosphere. An attempt was made to minimize the potential impact of interfering 
species while positioning the instrument band. This can be appreciated upon view- 
ing the effect on 0 3 ECFs from including interfering species for a spectral interval 
slightly larger than the chosen channel region, as is depicted in Figure 4.18. Clearly, 
it is desirable to avoid overlap with the CO 2 (~ 1053.92 and 1055.63 cm” 1 ) and H 2 O 
(« 1055.52 cm -1 ) spectral lines and their associated wings. It should also be noted 
that the impact of interfering species would be equivalent to these calculations if and 
only if their existence was ignored; whereas, in reality the impact would be much less 
since significant knowledge pertaining to the interferants should exist. The potential 
effect from uncertainties in knowledge of the interferants will be addressed in the 
next chapter. 

4.6 Further Design Aspects and Signal-to-Noise Ratio Anal- 
ysis 

The light gathering power of an optical system is characterized by its etendue or 
throughput, which can be given (Wyatt, 1991) by 

U = An, (4.23) 

where A is the receiving area through some surface and fi is a projected solid angle; 
and the projected solid angle can be defined as the area of another surface which is 
intercepted (or viewed) by the receiving surface divided by the square of the distance 
separating the surfaces. Figure 4.19 presents a schematic of the instrument design 
configuration proposed for this research. Accordingly, the telescope throughput of 
this instrument can be represented by 

T T ^ J‘2 ^ n2 

U T - 1FOV ’ 


(4.24) 
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Figure 4.19: Schematic of proposed double-etalon instrument design configuration 
illustrating relative positioning of key optical elements. 
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with the instrument’s instantaneous field-of-view angle ( Ojfov ) obtainable from 

(4.25) 


,QlFOV x dfifF 

tan (— ) = 


dj is the telescope entrance diameter, d^F is the nadir footprint diameter, and z s 
corresponds to the satellite altitude. Similarly, an expression for the FPI throughput 
can be written as 


Ufpi = —DgClFPi, (4.26) 


with the Fabry-Perot interferometer solid angle CIfpi represented by 


FPI — 7rtan 


(4.27) 


De is the diameter of the etalons, and 6 max is the angle from the optical axis to the 
most outside fringe detected, which is determined by positioning of the linear array 
detector element farthest from the cone vertex. The angle 9 max can be given by 

tan 6 max = ^y£, (4.28) 

for fringe detection at a- maximum distance of r max from the cone vertex by an op- 
tical system having a fringe forming lens of focal length f. The projected solid angle 
can also be represented by fl = Trsin 2 a, where a is the beam convergence angle 
with respect to the optical axis at the throughput evaluation surface; this is a more 
useful expression when evaluating channel throughput. By assuming invariance of 
throughput for the optical system and evaluating the throughput at various points 
along the beam path, several key optical design parameters may be calculated. For 
example, by specifying the nadir footprint size, satellite orbit altitude, and telescope 
diameter, the telescope throughput may be computed using equation 4.24. This 
result combined with the invariance of throughput requirement can then be used 
with equation 4.26 to calculate the necessary diameter of the etalon working area. 
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However, evaluation of Qfpi necessitates knowledge of r mox which requires specifica- 
tion of the spectral extent of the measurement. Expressions relating these quantities 
can be developed by considering in further detail the geometry of the CLIO system 
configuration. Upon viewing the instrument’s optical and detection configuration in 
Figure 4.19, the following relationships can be realized: 

tan 6 C = y-, (4.29) 

tan dj = - - - , (4.30) 


tan 0j+ 1 = 



(4.31) 


6 C , 6j, and 0j+\ correspond to the angles between the optical axis and the detector 
center, starting, and ending positions, respectively; rj is the radial distance from the 
cone vertex to the detector center along the optical axis (or to the detector’s central 
spectral position in the fringe plane), A rj is the detector spatial length, and f is the 
focal length of the fringe forming lens as previously mentioned. Using equation A. 3 
we can write an expression relating the observation wavenumber and angle from the 
optical axis as 


;/ = - = — , (4.32) 

2 fit cos 9 cos0 

where Vo is the wavenumber for 9 = 0°. This equation enables a more useful rela- 
tionship for 9 max to be formulated as 


/ 

cos 6 max = ° , (4.33) 

Vo + Av 0 bs 

where Av 0 b s is the desired spectral extent of the measurement. Equation 4.32 also 
permits development of an expression for the wavenumber extent (or spectral reso- 
lution) of the linear detector array elements (or channels) to be given by 



_l 

cos 0j +i 


1 

cos Oj 


]• 


(4.34) 
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As can be seen from the above equations, the lengths (in the optical axis direction) 
of the detector array elements determine the spectral resolutions of the instrument 
channels. Whereas, the minimum necessary detector widths are determined by the 
detector plane spatial energy distribution perpendicular to the optical axis resulting 
from the CLIO system. The 90° kaleidoscope mirrors fold the fringe image (formed 
by the fringe forming lens) on top of itself, thus reducing the image area by a factor of 
4. As discussed in Section 3.3, the view from the fringe plane is composed of 3 virtual 
images along with the real image (see Figure 3.7); this increases the etalon effective 
diameter to that given by equation 3.5. As with throughput, the relative aperture or 
f-number (F # ) is a measure of the flux-gathering ability of an optical system. The 
f-number is given by the ratio of the effective focal length to the entrance aperture 
diameter; for this instrument we can write 

f* = 4r> < 4 - 35 > 


for the FPI f-number. However, the kaleidoscope mirrors cause the “system f- 
number” to be reduced since D e ff > De , as given by 


F#.v> = 


/ 


D 


(4.36) 


eJJ 


A smaller f-number usually implies a better system, however, the associated cost 
and complexity level usually increase as well. An F# syf = .5 can be considered as a 
theoretical minimum, while F# sys < 1 is generally difficult to achieve in practice. In 
addition, the expressions for throughput given above become less accurate for lower 
f-number systems; however, they are sufficient for most feasibility studies and will 
be used in the present study with F# syr « 1.2 set as a design goal. 

Following the kaleidoscope, the circle to line converter converts a 90° arc of the 
concentrated fringe image into a line directly below in the detector plane, while 
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leaving the radial (or spectral) direction unchanged. The length of this “line” which 
varies with distance from the cone vertex is the minimum width for the detector 
array element at a particular location and can be determined by assuming invariance 
of throughput. As indirectly mentioned earlier, the invariance of throughput can be 
stated as A 7 rsin 2 a = constant. We can apply this condition at the “imaginary” 
fringe plane (which would exist in absence of the circle to line converter) and at 
the plane containing the detector array elements. For a particular channel in the 
fringe plane, the area can be given by Af } = — where the factor 

of 4 accounts for the kaleidoscope impact; in addition, sin 2 a can be approximated 
by — . Alternatively, the area of the rectangular detector array elements in the 

#5J/5 

detector plane can be represented by Ad } = A rjWj, where Ar, is the length and Wj 
is the width; and a = j represents the half cone angle. Equating these throughputs 
and solving for the minimum necessary width, one obtains 

7 r / 2 [tan 2 9j + 1 — tan 2 9j] 




8 FI Ar, 


(4.37) 


The diffraction limitation must be considered in the design process for any optical 
system to ensure sufficient image quality. The quality of an optical system is usually 
expressed in terms of image size for a point source. The image of a point source is 
not a point but is usually referred to as an optical blur. The blur circle for a perfect 
optical system is called an “Airy disk,” and is composed of a central blur of light 
surrounded by alternating dark and light rings of rapidly decreasing intensity. The 
size of the blur circle is usually defined by the diameter of the first dark ring as 


Dbiur = 2.44A-/- = 2.44A (4.38) 

L>eff 

where A is the wavelength of radiation; f, D e jf and F# sys are the lens focal length, 
effective system diameter and f-number, respectively, as previously defined. The 


91 


Characteristics of Assumed Orbit 

Orbit: 

Sun-synchronous; altitude z s = 955 km 

Scan: 

±51° from nadir 

Scan rate (s): 

lscan | oc scan 

85ec * sec 

Instantaneous field-of-view (IFOV): 

3 °; 50 km at nadir 

^pixels / scan (n): 

1 + 12 «!2 - 35 

Step dwell time: 

t td = £ = 228.57 ms 

Dwell (integration) time: 

td 22 220 ms 


Table 4.5: Orbital characteristics and dwell time determination assuming synchro- 
nization with the Nimbus-7/T0MS orbit. 

angular size of the blur can be given by the ratio of the blur diameter to the effective 
focal length as 

2 44A 

(4.39) 

L) e }J 

In order to minimize the loss of energy due to diffraction, the angular width of the 
narrowest fringe 

= [*i + . - (4-40) 

A Or . 

must be larger than the blur circle. Therefore, the condition that o b ™ lTi ~ > 1 must 
be verified in the design process. 

Adequate signal-to-noise ratio (SNR) is essential for the successful design of any 
electro-optical system intended for information processing. The SNR can be evalu- 
ated from the radiometric performance equation, which for a uniform extended-area 
source can be expressed (Wyatt, 1991) as 

(4 - 41) 


D* (X 1.E+10 cm Hz ,/2 / W) 
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1054.23 1054.56 1054.90 1055.23 

Frequency (cm'’) 


Maximum achievable SNR 



1054.23 1054.56 1054.90 1055.23 

Frequency (cm' 1 ) 


Telescope diameter 
d T = 9.0 inches 

Nadir footprint size 
d NF = 50.0 km 

Integration (or dwell) time 
t d = y 220.0 Yns 


Objective lens focal length 
f- o9.0 cm 



Number of passbands= 1 
Optical ^ransmissivity 

Chopping factor /? = 2.50 


E f>°- ", 


diameter 

cm 


Effective (system) diameter 
D #ff = 32.D03 cm 

A0, = 8.554e-Q4 radian 

6 k, Tin = 7.229e-05 radian 

blur => ratio = 11.834 


Frinqe forming lens 
F|= 2.942 * 

Effective ^s^m) 

#»y» 

Spectral resolution 
FPI: 0.0540 cm- 1 
apertgre: Q.Q4Q4 cm— 1 
effective: 0.0680 cm-1 


Figure 4.20: Results from a signal-to-noise ratio (SNR) analysis using the U.S. Stan- 
dard (1976) atmosphere. Spectral locations for channels are indicated 
by *. 
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Channel 

j/o c (cm -1 ) 

9(radian ) 

r(cm) 

A r(cra) 

W(cra) 

SNR 

i 

1054.240 

1.07e-02 

0.42 

1.42e-01 

2.17e-01 

169.97 

2 

1054.295 

1.48e-02 

0.58 

1.01e-01 

3.03e-01 

156.44 

3 

1054.336 

1.72e-02 

0.67 

8.70e-02 

3.54e-01 

162.99 

4 

1054.377 

1.93e-02 

0.75 

7.73e-02 

3.98e-01 

206.50 

5 

1054.418 

2.12e-02 

0.83 

7.03e-02 

4.38e-01 

266.90 

6 

1054.459 

2.30e-02 

0.90 

6.49e-02 

4.74e-01 

289.88 

7 

1054.513 

2.51e-02 

0.98 

5.94e-02 

5.18e-01 

247.79 

8 

1054.554 

2.66e-02 

1.04 

5.61e-02 

5.49e-01 

211.52 

9 

1054.595 

2.81e-02 

1.09 

5.32e-02 

5.79e-01 

239.65 

10 

1054.636 

2.94e-02 

1.15 

5.08e-02 

6.07e-01 

272.47 

11 

1054.677 

3.07e-02 

1.20 

4.86e-02 

6.33e-01 

247.51 

12 

1054.718 

3.19e-02 

1.25 

4.68e-02 

6.59e-01 

191.67 

13 

1054.759 

3.31e-02 

1.29 

4.51e-02 

6.84e-01 

164.57 

14 

1054.800 

3.43e-02 

1.34 

4.36e-02 

7.07e-01 

180.02 

15 

1054.841 

3.54e-02 

1.38 

4.22e-02 

7.30e-01 

189.86 

16 

1054.882 

3.65e-02 

1.42 

4.09e-02 

7.53e-01 

163.44 

17 

1054.923 

3.75e-02 

1.46 

3.98e-02 

7.74e-01 

139.08 

18 

1054.964 

3.86e-02 

1.50 

3.87e-02 

7.96e-01 

133.20 

19 

1055.016 

3.98e-02 

1.55 

3.75e-02 

8.22e-01 

148.72 

20 

1055.057 

4.08e-02 

1.59 

3.66e-02 

8.41e-01 

179.91 

21 

1055.098 

4.17e-02 

1.63 

3.58e-02 

8.61e-01 

206.16 

22 

1055.139 

4.26e-02 

1.66 

3.50e-02 

8.80e-01 

215.83 

23 

1055.180 

4.35e-02 

1.70 

3.43e-02 

8.99e-01 

200.74 

24 

1055.221 

4.44e-02 

1.73 

3.36e-02 

9.17e-01 

183.83 


Table 4.6: Summary of detector element positioning and achievable SNR (assuming 
D* — 5.5 x \0™cmHz5 /W) for the chosen channels. 
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D* (D-star) is the detectivity, in cmH z* /W; r,- is the instrument peak optical trans- 
mittance; A 3 and Cl s are the area (cm 2 ) and solid angle (sr), respectively, of the 
source, so that represents the source throughput; /? is the chopping factor; Ad 

(= ArjWj) is the detector area, in cm 2 ; A / is the electrical information-noise band- 
width, in s -1 ; I(u) is the upwelling radiance (W cm ~ 2 sr~^ / cm~ l ) and 1Zi(t/) is the 
instrument relative spectral response, so that I(u)'R.i(v , )du represents the effective 
flux, in Wcm~ 2 sr~ l . Applying equation 4.41 to this instrument feasibility study, we 
can replace j v I(v)1li(v)dv with S(v) (equation 4.19), the FPI signal. By invariance 
of throughput, the source throughput can be replaced by the channel through- 

put U(v) = A^i^Trsin 2 (^). In addition, the information-noise bandwidth may be 
expressed (Wyatt, 1991) as A / = where t d is the scene integration or dwell 

time. We can now reformulate the radiometric performance equation (equation 4.41) 

for the feasibility study in this research as 

F ) * T iArf(j/)7rsin 2 ( Y) / l/h, » h 


„ , U TiAdxvmsm ( t ) r^^gh , , i i 

SNR(v) = . = -^~ / Td(u — v )Tf(v )I(v )dv . 


(4.42) 


At the expense of slightly additional mechanical complexity, the detectivity is de- 
sired to be as high as reasonably possible to enable achieving maximum SNR; ac- 
cordingly, this design will require detector cooling for improved detectivity, in addi- 
tion to cold stop limitation of unnecessary portions of the detector element FOV’s 
to block unwanted background radiation from being detected. Mercury Cadmium 
Telluride (HgCdTe) photoconductive (PC) detectors perform quite well in this spec- 
tral region; cryogenic operation at 77 K utilizing liquid nitrogen ( LN 2 ) cooling and 
60° FOV cold stop shielding will be assumed in this work. A detectivity value of 
D* = 5.5 x 10 10 cm//z2 jW will be implemented into the feasibility study calculations. 
Many detector manufacturers can meet or exceed this D-star value (i.e., the J15D 
series by EG&.G Judson) for the narrow bandwidth, small detector element, limited 
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FOV, linear detector array configuration being considered for this 10 /xm region. A 
chopping factor of fi = 2.5 will be used. The dwell time is defined as the time that 
the system dwells on a single resolution element and can be given approximately by 
t d « Jj, where n is the number of resolution elements per scan and s is the scan 
rate. Assuming a 50km nadir footprint and orbital parameters similar to those of 
the Total Ozone Mapping Spectrometer (TOMS) instrument aboard Nimbus-7 (Greg 
et al., 1991) to achieve continuous coverage at the equator, yields a dwell time of 
t d K 220 ms. The orbital characteristics and dwell time determination for this study 
are summarized in Table 4.5. Figure 4.20 represents the results of this SNR analysis. 
The top plot depicts the minimum D* necessary to achieve a SNR=100. Whereas, 
the maximum achievable SNR for D* = 5.5 x \0 lo cmH z? /W is illustrated in the 
bottom plot; these values will be assumed for the error analysis of the next chapter. 
Table 4.6 summarizes the design logistics for the chosen channels. Note that all 
suggested instrument design specifications are contained within Tables 4.3 and 4.6 
as well as Figure 4.20. 


CHAPTER V 


RETRIEVAL SIMULATION AND ERROR 

ANALYSIS 


5.1 Introduction 

Retrieval of atmospheric ozone will be performed through implementation of a 
nonlinear least squares iterative technique to achieve convergence of the “measured” 
and “calculated” forward model signals. In particular, the maximum likelihood 
method will be implemented in an approach similar to that discussed by Rodgers 
(1976, 1984); this is a statistically based technique which yields the most probable 
solution. For the error analysis, a procedure similar to Rodgers (1989, 1990) will be 
implemented. 

5.2 Ozone Retrieval Technique 

Retrieval theory enables the derivation of atmospheric parameters from measured 
values of radiation emerging from the top of the atmosphere, the intensity and spec- 
tral distribution of which depend on the state of the atmosphere according to known 
physics. The task is to deduce the best estimate of the state of the atmosphere or the 
particular parameters of interest from the radiation measurements. As discussed by 
Rodgers (1976), there are two distinct aspects to this problem, namely the “inverse” 
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problem and the “estimation” problem. The inverse problem involves inverting a 
known equation that expresses the measured radiation as a function of the atmo- 
spheric state, into a formulation for the atmospheric state in terms of the radiation; 
however, this result usually has no mathematically unique solution and is said to be 
“ill-conditioned.” Nonuniqueness in the solution arises from several factors. There 
are vertical regions of the atmosphere which provide minimal or no contribution to 
the measured upwelling radiances, and therefore their corresponding constituent con- 
centrations could be limitless; a priori constraints can be used to impose boundaries 
on their sizes. In addition, measurement noise in the observations implies a corre- 
sponding nonuniqueness in the solution; the magnitude of which can be estimated by 
mapping measurement error bars onto error bars of the solution. Nonuniqueness is 
also introduced into the solution from a finite vertical resolution, which results from 
a finite spectral resolution of the instrument as well as a consequence of the nadir- 
viewing geometry; note that even with infinite spectral resolution, finite vertical 
resolution would still exist when employing a nadir-viewing satellite configuration. 
The process of determining the best solution from all the possible ones constitutes 
an estimation problem. 

A constrained least squares solution, such that the sum of the squares of the 
differences between the measurements and the calculated signal values corresponding 
to the solution is minimized, is commonly implemented to solve an ill-conditioned 
problem. A priori constraints can refer to any type of constraints on a solution and 
they can be regarded as “virtual measurements” (Rodgers, 1976), since they contain 
information pertaining to the unknown profile similarly to the direct measurements. 
If nothing is known about the unknown function, one can obtain an inversion using 
some sort of objective constraint, such as minimization of variance (Twomey, 1977). 
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If statistics are known, then one can constrain accordingly. For example, from a 
statistical constraint approach, a climatological mean profile and its covariance can 
be considered as a measurement and its corresponding uncertainty; and a constrained 
solution can be found which minimizes the mean square departure from this mean. 
This method tends to push the indirectly sensed solution towards conformity with a 
body of past data possibly obtained by more direct methods. This formula is very 
useful when there is justification for choosing the statistical constraint. 

The atmospheric remote sounding retrieval problem is usually nonlinear; con- 
sider, for example, the extraction of an atmospheric constituent vertical profile from 
equation 4.14. In general, we want to know the value of an unknown vector x de- 
scribing the state of a particular atmospheric parameter; this is usually a vertical 
profile of some quantity given at a finite number n of levels (i.e., the O 3 mixing ratio 
profile). We measure a signal vector y containing m elements and possessing an 
m x m error covariance matrix S f , which is related to x by 

y = F(x), (5.1) 

where F is a nonlinear forward model. Often m <C n, making the inverse problem 
formally ill-posed; as is the case for this study with m = 24 (the number of instru- 
ment channels) and n = 46 (the number of levels used to represent the vertical O 3 
profile). Assuming we have a priori information Xo about x with an n x n error 
covariance matrix S x , the maximum likelihood condition (Rodgers, 1976) requires 
that x minimizes 

[F(x) - y] T S ( -1 [F(x) - y] + [x - x 0 ] T Sj 1 [x - x 0 ]. (5.2) 

Although usually difficult to implement, the solution can be found by finding x 
that minimizes equation 5.2. However, easier approaches are usually available if 
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the problem is not too nonlinear, or if good a priori information is available in the 
very nonlinear case. For example, nonlinear problems are frequently solved with 
some form of Newtonian iteration, which requires the forward model to be expanded 
as a Taylor series about a guessed value x t of the solution. This enables a direct 
measurement of y = F(x) to be expressed as 

y = F(x t ) + (x - x,)^| Xi , (5.3) 

where higher-order terms have been neglected by using a truncated Taylor series. 
This may be reformulated as 


y = y. + K*,(x - x,), 


(5.4) 


through defining the calculated signal for x, as y ; = F(x t ) and the forward model 
derivative by 


K _d_F 

Xi ~ <9x lx -’ 


(5.5) 


K x is an m x n matrix which is sometimes called the weighting function. This results 
in an under-constrained set of equations for x. Considering the a priori constraint 
Xo with covariance S x as a virtual measurement yields the following set of equations: 


K Xi x = y - y t + K Xi x, with cou(K Xi x) = S c , (5.6) 


and 

x = x 0 with cov(x) = S x . (5.7) 

Two independent measurements of a scalar quantity are usually combined by 
taking a weighted average, using the reciprocal of the variances as weights. Using 
the same principle applied to vector quantities, the direct measurement and virtual 
measurement of equations 5.6 and 5.7, respectively, can be combined to get the best 
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estimate of x (the O 3 mixing ratio profile). As shown by Rodgers (1976), this solution 
can be expressed by an iteration equation as 

x t+1 = Xo + S x Kl(K Xi S x Kl + S f ) _1 [y - y, - K x ,(x 0 - x,)], (5.8) 

with a covariance of the solution given by the n x n matrix 

S = (S;> +KjSr 1 K r )- 1 . (5.9) 

Equation 5.8 can be implemented in a sequential fashion when applied to actual 
satellite data, whereby the new estimate of x made at each measurement stage is 
used as the a priori vector for the next stage. The concept of sequential estimation 
allows one to use the continuity that exists along a subsatellite track to improve 
the accuracy of retrieval. However, in practice it can be difficult to estimate S x . 
If knowledge of how the profile should change horizontally is available (i.e., through 
analysis of global climatology), then this can be combined with the retrieval previous 
in time and therefore close in space to formulate a better a priori estimate using an 
extrapolative process. A best estimate of S x can be obtained by a statistical analysis 
of the a priori data set, for example, S x can be defined to be the statistical variance 
of the climatological average profile at the horizontal point of interest. As is obvious 
from equations 5.8 and 5.9, changing S x due to uncertainty in its appropriate value 
directly affects the estimated solution x and its corresponding covariance S. 

5.3 Method of Error Analysis 

A procedure similar to that discussed by Rodgers (1989, 1990) will be followed 
for the analysis of error in this research. The total error can be separated into 
components due to instrument measurement noise, portions of the vertical profile 
which cannot be seen by the observing system, and uncertainty in model parameters. 
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The “observing system” refers to a combination of both the measuring instrument 
and the retrieval technique; hence an actual atmospheric constituent profile is the 
input to the observing system, while a retrieved profile is its output. 

The measured signal vector in equation 5.1 can be re-expressed as 

y = F(x,b) + e y , (5.10) 

where b is a vector of forward model parameters (such as instrument calibration char- 
acteristics, spectral line data, surface characteristics, etc.) which are not perfectly 
known, and the error vector t y represents measurement uncertainty. The uncertainty 
in the forward model parameters can be a source of random or systematic differences 
between measured and calculated values of y. It is assumed that the forward model 
accurately represents the physics of the measurement process, so that forward model 
errors are due solely to uncertainty in forward model parameters; note that these 
model parameters can include quantities that are closely coupled to the state vector 
x itself (i.e., spectral line parameters of the target species) but are not being deduced 
in the retrieval process. 

The measurement error e y directly contributes to error D y e y in the retrieval pro- 
cess, where D y is an n x m matrix representing the contribution functions and may 
be written (Rodgers, 1976) as 

D y = (S; 1 + K T x s;'K x )-'K T x s;\ (5.11) 

This corresponds to an error covariance matrix of 

S m = DyS ( Dj\ (5.12) 

If the measurement channels (or elements of the signal vector y) are independent, 
then S ( is a diagonal matrix. As defined by Rodgers (1990), the “null-space error” 


102 


results from regions of profile space which cannot be measured by the observing 
system. Assuming the mean profile to be used for x 0 as an a priori constraint, the 
null-space error for any single retrieved profile can be given by (A x — I)(x — x); 
A x is the averaging kernel matrix defined by A x = D y K x , I is the identity matrix, 
and x represents the mean state so that x - x corresponds to the true profile’s 
deviation from the mean state. Knowledge of the statistics of x — x is necessary to 
reasonably estimate the null-space error; and, in particular, a covariance matrix S x 
of its variability is required. The null-space error can then be expressed by 

S n = (A x — I)S X (A X — I) T . (5.13) 

As previously mentioned, errors e b = b - b in the forward model parameters can 
be random or systematic in nature, depending on their source of origination. A 
vector b can be defined as a best estimate of the model parameters, having an error 
covariance matrix S b , which is used in the retrieval process; whereas, the measuring 
system uses exact values for the model parameters b. S b is an nb x nb matrix, where 
nb is the number of uncertain forward model parameters which are being considered 
as contributors to the overall forward model error; in addition, the vectors b and b 
both contain nb elements. This introduces a model parameter component into the 
retrieval error of D y Kf,£6, with a covariance matrix given by 

S s = AfcSfcAjf; (5.14) 

the measurement sensitivity to forward model parameters is represented by the mxnb 
matrix K(, and is defined by 

K b = — ; (5.15) 

b db' y 

and the n x nb matrix A b = D y K t corresponds to the state parameter sensitivity to 
model parameters b. Note that S m , S„, and S s are all n x n covariance matrices. As 
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Atmospheric model tot 03 (DU) trop 03 (DU) 

1) tropicol (15N onnuol ovg.) 281.7 36.6 

2) mid-latitude summer (45N July) 334.0 56.2 

3) mid— latitude winter (45N Jan) 377.7 42.1 

sub-arctic summer (60N July) 347.3 40.5 

5) sub-arctic winter (60N Jan) 375.1 24.4 

6; U.S. standard (1976) 343.8 36.9 


Figure 5.1: Atmospheric models from the 1976 Standard Atmosphere database. 


discussed above, the overall error in the retrieved profiles results from a combination 
of three components: measurement error, null-space error, and forward model error. 
The total random error can then be formulated by summing the null-space and 
measurement errors, in addition to any random components of the forward model 
error. 


5.4 Application to Atmospheric Cases 


An atmospheric database containing temperature, pressure, density, and gas vol- 
ume mixing ratios (Anderson et ah, 1986) of the necessary chemical species as a 
function of altitude for six reference atmospheres will be utilized for the retrieval 
and error analysis simulations in this research. The vertical profiles were formulated 
by a combination of measurements and theoretical predictions, and dayside estimates 
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for diurnally varying species (i.e., O3) were incorporated. The database includes only 
those twenty-eight minor and trace atmospheric gases which were part of the AFGL 
molecular line parameter compilation (Rothman et al., 1987) at the time of publi- 
cation; note that the most recent HITRAN database (Rothman, 1993) contains line 
parameters for thirty-one molecular species, however, as mentioned in Section 4.3 
these additional species are relatively insignificant for the radiative transfer within 
the 9.6 nm spectral region of the current study. The reference atmospheric model 
profiles consist of the following types of atmospheres: 1) tropical (15N annual av- 
erage), 2) mid-latitude summer (45N July), 3) mid-latitude winter (45N January), 
4) subarctic summer (60N July), 5) subarctic winter (60N January), and 6) U.S. 
standard (1976). Vertical profiles of temperatures and ozone volume mixing ratios 
for these atmospheric types are illustrated in Figure 5.1. 

The diagonal elements of a covariance matrix are the error variances of the el- 
ements of the corresponding vector which it represents; whereas, the off-diagonal 
elements indicate inter-level correlations existing between the variances at different 
levels. If the a priori information vector x 0 can be derived from an independent data 
set of high spatial resolution, then it should be possible to compute an accurate co- 
variance matrix S x ; however, such a data set is not always available and other means 
must be explored to arrive upon a reasonable estimate for S x . Realistic values for the 
variances at each level (or diagonal elements of S x ) are usually obtainable, whereas 
the inter-level covariances (or off-diagonal elements of S x ) may be somewhat harder 
to find. If the inter-level covariances are not available, then the variances of the 
different levels can be assumed independent and the off-diagonal elements of S x set 
to zero; alternatively, these off-diagonal elements ( S(j , 1 ^ j) can be estimated. As 
discussed by Rodgers (1990), by assuming a vertical length scale / of correlations, a 
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20 40 60 

Percent (%) 



80 100 180 200 220 240 260 280 300 

Degrees (K) 

Atmospheric model 

1) tropical (15N annual avg.) 

2) mid-latitude summer (45N July) 

3) mid — latitude winter (45N Jan) 

4) sub-arctic summer (60N July) 

5) sub-arctic winter (60N Jan) 

6) U.S. standard (1976) 


Figure 5.2: Ozone standard deviation profiles used in generating the a priori co- 
variance matrices corresponding to atmospheric models within the 1976 
Standard Atmosphere database. Note that atmospheric models 3) and 4) 
have similar tropopause heights (~ 12 km), making their corresponding 
standard deviation profiles overlap and indistinguishable (appearing as 
one very long-dashed profile). 
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formulation such as 

Si, = (SiiSjtf exp ^ -] (5.16) 

can be used to generate the elements of an artificial covariance matrix; z,- and Zj 
correspond to altitudes (in the same units as /) of levels i and j, respectively. Equa- 
tion 5.16 will be used to estimate the a priori covariance matrices in this study by 
assuming a vertical length scale of correlations 1 = 5 km for the O 3 covariances. 
The ozone variances at each level were arrived upon as follows. The uncertainty (in 
percent) of ozone concentrations in the Standard Atmosphere database (Anderson 
et al., 1986) is approximately constant throughout much of the stratosphere. At 
higher altitudes, the uncertainty increases primarily due to lack of measurements. 
And at altitudes below, local maxima exist in the lower stratosphere/upper tropo- 
sphere and lower troposphere regions due to uncertainty in tropopause height and 
surface ozone variability, respectively. Based on these criteria, smooth curves were 
generated to approximate these ozone variances for use as diagonal elements in the 
a priori covariance matrices in the different atmospheric simulation cases. Figure 5.2 
shows the resulting a priori ozone standard deviation profiles implemented within 
this analysis. Note the conservative approach used in the lower (than tropical at- 
mosphere) tropopause atmospheric cases, where the standard deviations for lower 
atmosphere regions have been arbitrarily increased in an effort not to under estimate 
tropospheric errors for these cases possessing shallow tropospheres. It should be men- 
tioned that while a more elaborate technique for S x generation is certainly possible, 
the author believes that the method implemented enables the a priori uncertainty 
to be adequately represented for a feasibility study at this stage in an instrument 
design process. 

The measurement channels will be assumed independent for this analysis. Ac- 
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Signal-to-noise ratio 

ranges 


Atmospheric model 

SNRm t n 

SNR max 

1) tropical (15N annual average) 

176.84 

392.00 

2) mid-latitude summer (45N July) 

159.01 

327.69 

3) mid-latitude winter (45N Jan) 

104.63 

206.79 

4) subarctic summer (60N July) 

159.97 

285.34 

5) subarctic winter (60N Jan) 

89.95 

158.23 

6) U.S. standard (1976) 

133.20 

289.88 


Table 5.1: Summary of signal-to-noise ratio range estimates corresponding to the 
instrumentation of this research applied to different atmospheric types. 
Values of SNR for each spectral channel are necessary for formulation of 
the measurement covariance matrices. 

cordingly, the measurement covariance matrix S f will be a diagonal matrix with 
elements equal to the variance of the corresponding channels which they represent. 
Assuming the measurement error to be normally distributed and having a zero mean, 
the diagonal elements of the measurement covariance matrix can be estimated for 
each channel i by 

S,„ = (^)V. (5-.T) 

in units of the measured signal squared (yf)', the signal-to-noise ratio for channel i is 
represented by SNR t . Values of SNRi , y t , and thus S (ii were adjusted to correspond 
to the particular type of atmosphere being simulated. Table 5.1 summarizes the 
range in signal-to-noise ratios corresponding to the six atmospheres simulated in this 
research project; see Figure 4.20 for the spectral dependence of the SNR pertaining 
to the U.S. standard (1976) atmosphere. 
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The partial derivatives composing K x and Kf, (equations 5.5 and 5.15, respec- 
tively) were evaluated numerically by perturbing the forward model and using finite 
forward differences on the results. For example, all rows in column p of K x can be 
generated from 


K 


*(.= 1-24, .,=p) 


F(x 0 + Sx 0 ) - F(x 0 ) 


Sx 


Op 


(5.18) 


where <5xo represents the perturbation to Xo and Sx o p corresponds to the perturbation 
magnitude at level p (the only non-zero element in 6xo). The matrix K x is necessary 
for both the retrieval process and error analysis, whereas Kj is needed only for the 
error analysis. In the error analysis, K x and K(, are required to be evaluated only 
once for each atmospheric state about which the error budget is being considered. 
Whereas K x could be updated at each stage of the iteration process for each retrieval, 
if desired; although this step is rarely necessary, since convergence can usually be 
achieved in just a few iterations without derivative updates by using a reasonable first 
guess. As mentioned earlier, the vector b (needed for Kb evaluation) includes those 
forward model parameters for which the effect of their uncertainty on the retrieved 
ozone amounts is desired. It is beyond the scope of this present study to include 
all possible parameters, however, an attempt has been made to consider those most 
important. In particular, contributions to the overall error will be assessed (in Sec- 
tion 5.6) for the following forward model parameters: water vapor continuum model 
coefficients, surface temperature and emissivity, instrument radiometric channel cal- 
ibrations, vertical profiles of atmospheric temperature and interfering species, and 
spectral line parameters. 


Altitude (km) 
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Ozone 



Atmospheric type tot 

1) tropical (15N annual avg.) 

2) BZV sounding (8/5/92) 


Temperature 



180 200 220 240 260 280 300 
Degrees (K) 


03 (DU) trop 03 (DU) 

281.7 36.6 

293.7 51.9 


Figure 5.3: 


Ozone and temperature profiles from the 1976 Standard Atmosphere 
(tropical) along with a sounding from Brazzaville, Congo, on August 5, 


1992. 


5.5 Retrieval of 0 3 for a Perturbed Tropical Atmosphere 


The retrieval technique implemented into this study incorporates measurement 
uncertainty and null-space error for one retrieved profile, whereas the error analysis 
estimates how good the retrieved profiles will be in a statistical sense (encompassing 
many such retrievals) while including the effects of measurement, null-space, and for- 
ward model errors. Additionally, the mechanics of the retrieval process are the same 
for different atmospheric situations; therefore, the retrieval ability will be demon- 
strated for only one case, a perturbed tropical atmosphere. A tropical atmosphere 
has been chosen due to the significant recent interest in the south Atlantic ocean 
region where elevated tropospheric ozone has been observed (Fishman et al., 1990). 


Pressure (mb) 
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Sounding data from Brazzaville, Congo, on August 5, 1992, provided by Dr. 
Jack Fishman (NASA Langley Research Center) and Mr. Vincent Brackett (SAIC 
Corporation, Hampton, Virginia), were used for generation of the “true” profile. 
Although the ozone sounding profile was at high vertical resolution, the vertical levels 
did not coincide exactly with all of those needed for input into the forward model; 
the necessary interpolation was handled using cubic spline interpolation (see, for 
example, Press et al., 1990). In addition, the sounding profile stops below 1 1 mb (~30 
km). Specification of the true 0 3 profile for altitudes above 30 km was performed by 
using a modified average tropical atmosphere profile, and connecting the two profiles 
with cubic spline interpolation; the final modified profile was then formulated by 
adjusting the average tropical atmosphere profile with a multiplicative factor until 
the total O 3 column abundance matched the value provided with the sounding data 
(which was obtained by an independent method). Note that the uncertainty of 
this independently obtained total ozone value is not currently known to the author; 
however, even if quite erroneous, it serves its purpose in this retrieval simulation as 
it separates the true and first guess profiles. The average tropical atmosphere O 3 
profile was then used for a first guess in the retrieval process. Figure 5.3 illustrates 
vertical profiles of ozone and temperature for the average tropical atmosphere and 
the Brazzaville (BZV) sounding. 

As discussed in Section 4.5, the O 3 energy contribution functions ( 0 3 ECFs) 
indicate the vertical regions from which information content arises in the measured 
signal; and their broadness near the maxima dictate approximate vertical length 
scales of resolution about the corresponding peak contributions. By examining the 
broad O 3 ECFs for this measurement (Figure 4.16), one can see that any hopes for 
vertical profiling capability are quickly reduced to expectations of slight layering or 
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simply integrated layer amounts, especially in the troposphere. However, as this 
retrieval simulation will demonstrate, in vertical regions where the character of the 
true profile is well represented by the a priori information, one can also expect good 
agreement between the retrieved and true profile shapes. This is made possible 
by the number of levels used to specify the O 3 input profile in the forward model 
(n = 46) being large, and actually much greater than simply the small number of 
levels that one might expect to achieve from viewing the O 3 ECFs; note that this 
can lead to undesirable oscillations in the retrieved profiles for bad choices of the a 
priori constraint, but should not significantly degrade the retrieved integrated ozone 
amounts. 

Figure 5.4 summarizes the results of the retrieval simulation for this tropical at- 
mosphere case with elevated tropospheric ozone. After two iterations, the retrieved 
integrated ozone amounts appear closest to the truth. However, the result of the fifth 
iteration is defined as the retrieved profile for this simulation, since convergence of 
signal residuals (i.e., root-sum-square (rss) of the signal residuals) has been used for 
determining retrieval completion. And as the rss signal residuals show, the outcome 
does not change with continued iteration since the signals have already been matched 
to within the measurement uncertainty (equation 5.6) while satisfying the a priori 
constraint (equation 5.7). As Figure 5.4 demonstrates, the retrieved profile matches 
the true profile quite well for this particular simulation in both profile character as 
well as integrated ozone amounts. The inability to match the profile variability in the 
lower troposphere is an expected consequence of the broad nature of the O 3 ECFs 
in this vertical region. This simulation serves to exemplify the successful implemen- 
tation of the maximum likelihood retrieval technique applied to this research and 
to a certain extent the capabilities of this particular measurement. Note that while 
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Figure 5.4: Retrieval simulation for a tropical atmosphere having enhanced tropo- 
spheric ozone. Sounding data from Brazzaville, Congo, on August 5, 
1992, used for true profile; the 1976 Standard Atmosphere (tropical) 
used for first guess. 
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these results are very encouraging, they pertain to one particular retrieval and are 
not representative of “average” expectations. One would expect these results to be 
better than the average due to the likeness in character between the first guess and 
true profiles. In addition, while measurement and null-space errors are incorporated 
into this retrieval, an error analysis must be performed to consider the impact of 
forward model uncertainties. Section 5.6 will give a more complete assessment of the 
overall capabilities of this measurement. 

5.6 Results of Error Analysis 

This section will address the effect of errors on the retrieved state of atmospheric 
ozone corresponding to measurement with the simulated instrument of this research. 
The error analysis technique (Section 5.3) will be implemented as discussed in Sec- 
tions 5.3 and 5.4. The total random error will be estimated by combining the null- 
space and measurement errors with random components of the forward model error; 
this can be thought of as the instrument measurement precision. A measurement 
accuracy will also be estimated by combining the effects of forward model systematic 
uncertainties. 

Accuracy can be defined as closeness to the true value, and it is generally limited 
by systematic errors; however, it can be thought of as depending on random errors 
to some extent, since a given accuracy implies a precision at least as good. Precision 
can be interpreted as how exactly a result is determined (regardless of what the result 
means), or simply as the reproducibility of a result; the precision is generally limited 
by random errors. Variability of the retrieved ozone solution will be determined by 
random independent errors as defined by measurement error, null-space error, and the 
random part of the forward model uncertainties (i.e., model parameter errors which 
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can vary randomly from one measurement to the next). The systematic components 
of the forward model error are constant and do not contribute to variability of the 
ozone solution vector x, they just make a non-zero ensemble mean of x true — x retrieved 
which will loosely be considered as accuracy in this analysis. 

Although all possible forward model errors have not been simulated, an attempt 
has been made to include those most significant. The following uncertainties will 
be considered as random components of the forward model error: surface emissivity 
e, ± 1%, surface temperature T 3 ± 1 K, temperature profile T(p ) ± 2 K (one level 
at a time), along with the interfering species H 2 0(p) ± 30% and C0 2 (p) ± 5%. 
While the systematic components of the forward model error will be represented 
by: water vapor continuum model self-broadening coefficient C s ± 10% and nitrogen- 
broadening coefficient C n ± 10%, radiometric signal calibration ±1%; sensitivity to 
uncertainties in spectral line parameters will be assessed with H 2 0 halfwidths ±5%, 
O 3 halfwidths ±5%, H 2 0 line strengths ±10%, and 0 3 line strengths ±1%. The 
magnitudes of these uncertainties were chosen based on anticipated knowledge of the 
forward model parameters, presumably obtained from some independent sources (i.e., 
laboratory measurements and calibrations, climatological data, independent sensors, 
etc.). Note that while the ± signs are shown with the errors, the effect of either 
+ or — each uncertainty will be demonstrated; however, since the actual signs are 
unknown and may be assumed to vary randomly among the individual errors, total 
errors will be estimated by taking a root-sum-square (rss) of the component errors 
(even for those systematic in nature). Some uncertainties included in the random 
error compilation (i.e., emissivity and temperature of the surface) might also possess 
systematic errors; their corresponding effects on retrieved ozone will be demonstrated 
but not included in the overall error estimates, which assume individual errors to be 


115 


either random or systematic. In addition, effects from other systematic offsets will 
be calculated for demonstration purposes. For example, a case showing the impact 
of completely ignoring the water vapor continuum will be simulated; thus illustrating 
its importance, especially in the tropical troposphere. Also, the impact of systematic 
offsets in input profiles (at all levels simultaneously) of temperature and interfering 
species will be evaluated. 

The retrieved ozone sensitivity to each uncertainty will be calculated as a func- 
tion of altitude, and will then be used in computing the impact on integrated ozone 
amounts for both the troposphere and total column abundance. In computing the 
integrated effect, the impact on ozone amount for each layer can be simply summed 
up for the systematic errors, whereas they may be combined by taking an rss (and 
thus reducing the impact on integrated amounts) for random errors since the signs 
may vary independently between errors at different levels (i.e., in the nominal tem- 
perature profile). 

The annual average tropical atmosphere used as a first guess in the retrieval sim- 
ulation of Section 5.5 will be the first case considered in this error analysis. We 
will start by examining some of the random components to the forward model error. 
Equation 5.14 can be used to estimate their combined effect on the ozone retrieval. 
In order to accurately determine the sensitivities to random uncertainties in these 
input parameter profiles, knowledge of their profile covariances (variances at each 
level as well as correlations among the different levels) is necessary. Note that the 
impact of random perturbations (of fixed magnitude) is a function of the number of 
levels at which the input profiles are specified in the forward model. Consequently, 
perturbation amplitudes deduced from a statistical analysis of typical input profiles 
would require adjustment if the number of vertical levels requiring specification in 
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Figure 5.5: Retrieved ozone sensitivity to perturbations of interfering species at spe- 
cific altitudes for a tropical atmosphere. 
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T+2K at selected altitudes 




Figure 5.6: Retrieved ozone sensitivity to temperature profile perturbations and sys- 
tematic offsets in some interfering species for a tropical atmosphere. 
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the forward model were different from that of the database from which the input 
profile covariances were derived. For this study, perturbation amplitudes were arbi- 
trarily chosen to represent forward model input profile uncertainties at the various 
model input levels without including inter-level correlations. Figures 5.5 a), 5.5 b), 
and 5.6 a) show the retrieval sensitivity to perturbations in knowledge of H 2 0, C0 2 , 
and temperature, respectively, at selected altitudes. Note that a perturbation at a 
particular altitude also affects the retrieved ozone amount at other altitudes due to 
the broadness of the O 3 ECFs. In addition, local maxima of perturbation effects 
exist in the stratosphere and lower troposphere due to peaks in the O 3 ECFs at 
these altitudes; whereas the local maximum in the upper troposphere region results 
from this vertical region having low O 3 information content (minima in O 3 ECFs) 
and high a priori information uncertainty (due to tropopause height variability). 
The sensitivities to errors in interfering species illustrated in Figure 5.5 show max- 
imum error in O 3 near the surface to result from perturbations slightly above the 
surface (~715 mb). The near-surface sensitivity is greater for perturbations in this 
above-surface region since they affect the amount of upwelling radiation (reaching 
the nadir-viewing satellite sensor) originating from this altitude and below; and the 
sensitivity decreases with further elevation in the troposphere as the temperature 
and concentrations of interfering species decrease. Note that this decrease happens 
sooner in the H 2 0 perturbation case, since the H 2 0 concentration decreases with 
increasing altitude more sharply than C0 2 • This feature is not exhibited in the tem- 
perature perturbation case because perturbations at the higher temperature (near 
the surface) have a larger impact, since both the Planck function and transmittance 
(see equation 4.14) are affected. These plots demonstrate the importance of knowl- 
edge pertaining to the H 2 0 and temperature profiles. The uncertainty in retrieved 
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Results of Error Analysis: 1 ) tropical ( 15 N annual average) 


Effect on retrieved integrated 0 3 amounts 



Accuracy (%) 

Perturbation 

(systematic errors) 


Troposphere 

Total 

H 2 O continuum model 

without H 2 O continuum 

71.73 

12.21 

Cs & Cn +10% 

-6.485 

-1.104 

surface emissivity 

e. -1% 

-12.99 

-2.187 

surface temperature 

T a + IK 

22.32 

3.755 

temperature profile 

T+1K all levels 

10.65 

-1.267 

T+ 1 K in troposphere & 

16.75 

5.531 

-IK in regions above 



interfering species 

CO 2 + 1% all levels 

-0.04610 

-0.01887 

CH 4 +10% all levels 

-0.0008151 

-0.0001381 

tf 2 O+ 10 % all levels 
#20+30% all levels 

- 12.21 

-0.001551 

-2.152 

-0.0002209 

502+30% all levels 

-0.0006869 

-0.0001187 


Table 5 . 2 : The effect of certain systematic perturbations on retrieved integrated ozone 
amounts for a tropical atmosphere. 

ozone due to a random error in an input profile can be estimated by combining (using 
an rss) the error profiles resulting from perturbations (of the desired error magni- 
tude) at each altitude of which the forward model input profile is specified. When 
examining these error profiles, one must keep in mind the broad nature of the 0 3 
ECFs pertaining to this measurement, which makes the integrated effects (effects on 
integrated ozone amounts) more meaningful quantities; these integrated effects will 
be presented in tabulated form later in this section. 

Figure 5.6 b) illustrates the effect on retrieved ozone from systematic offsets in 
some of the interfering species; this shows error contribution associated with lack 
of knowledge regarding C// 4 , A^O, * 9 ( 92 , and probably CO 2 to be insignificant at 
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Figure 5.7: Retrieved ozone sensitivity to perturbations in various forward model 
parameters for a tropical atmosphere. 
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Figure 5.8: Retrieved ozone sensitivity to perturbations in spectral line parameters 
and radiometric signal calibration errors for a tropical atmosphere. 
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Results of Error Analysis: 1) tropical (15N annual average) 


Effect on retrieved integrated O3 amounts 



Accuracy (%) 

Perturbation 


(systematic errors) 


Troposphere 

Total 

signal calibration 



all channels +1% 

17.97 

1.963 

channel # 1 + 1% 

-4.948 

2.609 

channel #2+1% 

-4.450 

0.7598 

channel #3+1% 

-5.347 

0.7472 

channel #4+1% 

-2.030 

-0.001413 

channel #5+1% 

4.852 

0.2224 

channel #6+1% 

6.412 

1.729 

channel #7+1% 

2.034 

0.3676 

channel #8+1% 

-0.2036 

0.004968 

channel #9+1% 

2.616 

-0.3632 

channel #10+1% 

9.612 

-0.6286 

channel # 11 +1% 

6.309 

-0.4318 

channel #12+1% 

1.078 

-0.03715 

channel #13+1% 

-2.385 

1.181 

channel #14+1% 

-0.7961 

0.7265 

channel #15+1% 

1.677 

-0.2140 

channel #16+1% 

-0.5309 

-0.7034 

channel #17+1% 

0.1599 

-1.718 

channel #18+1% 

1.337 

-2.109 

channel #19+1% 

2.715 

-2.252 

channel #20+1% 

-0.1729 

-0.6613 

channel #21+1% 

1.117 

0.6545 

channel #22+1% 

1.729 

0.5754 

channel #23+1% 

-0.4567 

0.2314 

channel #24+1% 

-2.361 

1.274 


Table 5.3: The effect of systematic perturbations of the radiometric signal calibration 
on retrieved integrated ozone amounts for a tropical atmosphere. 
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the uncertainty levels simulated, due to avoidance of line overlap in the channel se- 
lection process (see Section 4.5). The effects resulting from uncertainties in some 
other forward model parameters are exhibited in Figure 5.7. Uncertainties in surface 
emissivity or temperature, such as those shown in Figure 5.7 a), can contribute to 
systematic or random errors (or both) depending on their source of origination; ac- 
cordingly, how they contribute to error in integrated O 3 or the method of combining 
with other errors varies. The radiometric signal calibration error (which has been 
applied to all channels simultaneously) and the temperature profile offsets would con- 
tribute to the systematic error, if present. Note that the case containing both positive 
and negative temperature biases for the troposphere and stratosphere (and regions 
above), respectively, corresponds to what is believed to be a worst case temperature 
error scenario. Such anticorrelated offsets could arise from a temperature retrieval 
scheme trying to compensate for a temperature error in one particular vertical re- 
gion by introducing an error of opposite sign in another region, although the regions 
would probably be of smaller vertical length scale than shown here. The magnitudes 
of all potential errors in Figure 5.7 a) are significant and, accordingly, effort should 
be made to minimize chances of their existence. Figure 5.7 b) illustrates retrieval 
sensitivity to systematic offsets in forward model parameters related to water vapor. 
The case neglecting the H 2 0 continuum demonstrates its importance to the radia- 
tive transfer within this spectral region for a (moist) tropical atmosphere. However, 
since this is known, a water vapor continuum model has been incorporated into the 
forward model; therefore, neglecting it is an unrealistic source of error. Alterna- 
tively, the sensitivity to H 2 0 continuum model coefficients and the actual amount of 
H 2 0 must be considered. The self- and nitrogen-broadening coefficients are generally 
thought of as being good to within 10% (see, for example, Zhao, 1992). Figure 5.7 
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b) illustrates the impact of such an uncertainty in these coefficients, individually and 
for both simultaneously. The effect of an offset in the H 2 0 profile is also shown. 

The effect on integrated O 3 amounts of each of the mentioned errors acting indi- 
vidually will be presented, and for a selected subset of these errors a total error impact 
will be estimated. Table 5.2 shows how some of the previously discussed systematic 
errors would affect the accuracy, if they were to exist. As already mentioned, the 
simulation without the water vapor continuum is merely for demonstration purposes. 
The simultaneous perturbation of the self- and nitrogen-broadening coefficients also 
fits into this demonstration category, since it is probably more likely that uncer- 
tainty in these coefficients would vary independently of each other. Uncertainties in 
the surface emissivity, surface temperature, temperature profile, and H 2 0 profile will 
be assumed random in nature for compiling the total error estimates; however, as 
represented in Table 5.2, systematic offsets in these forward model parameters can 
severely impact the accuracy of the retrieved ozone amounts. This table also shows 
reasonable uncertainties on some other interfering species to be insignificant, as was 
discussed earlier regarding Figure 5.6 b). Figure 5.8 a) illustrates the retrieved ozone 
sensitivity to errors in spectral line parameters. In particular, sensitivities to typical 
uncertainties in halfwidths and line strengths for both O 3 and H 2 0 are demonstrated. 
The impacts on retrieved ozone are negative (in sign) since a decrease in O 3 con- 
centration is required to offset the increased absorption which results from each one 
of these spectral line parameter errors shown. The radiance calibration error shown 
earlier (Figure 5.7 a)) assuming an error in all channels simultaneously is shown in 
Figure 5.8 b) along with the impact of the same +1% perturbation acting alone for 
each of the 24 channels. Note that average pressures from layers which peak energy 
contributions originate for the instrument channels are contained within Table 4.4. 
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Table 5.3 summarizes the effect on accuracy of retrieved integrated ozone amounts 
with these radiometric signal calibration errors included as systematic offsets. The 
maxima and minima error effects occur at different altitudes for the different chan- 
nels perturbed separately, since the vertical distribution of 0 3 information content 
varies with the various channels. Considering these radiance calibration errors to act 
independently would lessen their net error impact, however, a more conservative (and 
realistic) approach is to assume a systematic offset for all channels simultaneously. 
Note that the erroneous increase in retrieved tropospheric ozone from the positive 
perturbation of all signal channels simultaneously results from the algorithm increas- 
ing absorption by ozone (of the large surface signal contribution) to counteract the 
artificial signal increase. 

As discussed in Section 5.3, the total random error is composed of the random 
parts of the forward model error, along with the measurement noise and null-space 
error. Figure 5.9 illustrates the effects of these main random error components and 
the total (rss) random error as a function of altitude for this tropical atmosphere 
case. The measurement noise and null-space error components were computed from 
equations 5.12 and 5.13, respectively. The random components of the forward model 
error being considered in this analysis have already been discussed in this section. 
Figure 5.9 shows the net forward model random error, formulated by taking an rss 
of its components. This total forward model random error is the most significant 
random error impacting the troposphere for this tropical case. 

The next step of this error analysis is to estimate the impact of these assumed 
systematic and random errors on the overall accuracy and precision, respectively, 
for this measurement. Table 5.4 summarizes the component and composite errors 
on retrieved integrated ozone amounts for the tropical atmosphere case under eval- 
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Profile random errors 



Figure 5.9: Effects of component and total random errors on retrieved ozone for a 
tropical atmosphere. 


Pressure (mb) 
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Results of Error Analysis: 1/ 

tropical (15N annual average) 

Effect on retrieved integrated 0 3 amounts 


Accuracy (%) 

Precision (%) 

Perturbation 

(systematic 

errors) 

(random errors) 


Troposphere 

Total 

Troposphere 

Total 

H 2 0 continuum model 





Cs +10% 

-5.186 

-0.8828 



Cn +10% 

-1.314 

-0.2239 



surface emissivity 





€5 ± 1% 



3.734 

0.5333 

signal calibration 





all channels + 1% 

17.97 

1.963 



surface temperature 





T s ± IK 



6.415 

0.9162 

temperature profile 





T ± 2k one level at a time 



2.574 

1.006 

interfering species 





H 2 0 ± 30% one level at a time 



3.470 

0.5087 

C0 2 ± 5% one level at a time 



0.02772 

0.005789 

spectral line parameters 





H 2 O halfwidths +5% 

-0.1448 

-0.03795 



O 3 halfwidths +5% 

-8.379 

-3.103 



H 2 O line strengths +10% 

-0.2331 

-0.1155 



0 3 line strengths +1% 

-1.245 

-0.9663 



Error impact (rss) 





forward model uncertainties 

20.61 

4.024 

8.589 

1.548 

measurement noise 



1.978 

1.206 

null-space error 



5.514 

2.648 

total (rss) error 

20.61 

4.024 

10.40 

3.298 


Table 5.4: Overall errors in retrieved integrated ozone amounts corresponding to ap- 
plication of the proposed instrumentation to a tropical atmosphere. 
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uation. Therefore, based upon the uncertainties assumed and the error compilation 
procedure implemented, one would expect this instrument to enable (tropical) tropo- 
spheric ozone determination to within ~ 10% and ~ 20% for precision and accuracy, 
respectively; and knowledge of total column ozone abundance to within ~ 3% and 
~ 4% for precision and accuracy, respectively. Note that this analysis has been per- 
formed for the annual average tropical atmosphere and that tropical atmospheric 
cases of most interest would probably contain elevated tropospheric ozone, for which 
the sensitivity to these uncertainties and the overall error in retrieved ozone amounts 
would be reduced. Also, it may be possible to perform the signal calibration to better 
than 1% and thus significantly improve the accuracy of this measurement. 

This error analysis procedure was repeated for the remaining five atmospheric 
types introduced in Section 5.4. A brief discussion of these final results will now 
follow. Summary tables similar to Tables 5.2 and 5.4 have been formulated for each 
of the additional atmospheres and are included in Appendices B and C, respec- 
tively. Upon examining the sensitivities to systematic error perturbations of the 
various atmospheres presented in Appendix B, one can observe some expected re- 
sults. For example, the water vapor continuum drastically decreases in importance 
as the focus goes from a warm, moist tropical atmosphere all the way to a cold, 
dry subarctic winter atmosphere. Also, the significance of uncertainties in surface 
emissivity and temperature decreases as one considers atmospheres less sensitive to 
surface contribution variations (i.e., a winter atmosphere). Appendix C presents the 
total error budget (as done for the tropical atmosphere) for the other atmospheric 
cases. Some notable observations are as follows. The sensitivity to signal calibration 
errors increases for atmospheric cases with larger surface signal contributions (i.e., 
warm atmospheres), for which a perturbation of total signal becomes a larger part 
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Results of Error Analysis: total (rss) errors 



Effect on retrieved integrated O 3 amounts 


Accuracy (%) 

Precision (%) 

Atmospheric type 

(systematic errors) 

(random errors) 


Troposphere 

Total 

Troposphere 

Total 

1) tropical (15N annual average) 

20.61 

4.024 

10.40 

3.298 

2) mid-latitude summer (45N July) 

15.12 

3.278 

7.730 

3.311 

3) mid-latitude winter (45N Jan) 

11.94 

1.987 

8.627 

3.681 

4) subarctic summer (60N July) 

15.31 

2.932 

8.224 

3.317 

5) subarctic winter (60N Jan) 

7.549 

2.235 

10.32 

3.945 

6) U.S. standard (1976) 

16.38 

2.779 

9.656 

3.421 

mid-latitude [2) & 3)] average 

13.53 

2.633 

8.179 

3.496 

subarctic [4) 5)] average 

11.43 

2.583 

9.272 

3.631 

overall [!)— >5)] average 

14.11 

2.891 

9.060 

3.510 


Table 5.5: Summary of total (rss) errors in retrieved integrated 0 3 amounts for all 
atmospheric cases considered. 

of the “ozone signal.” The null-space error results from vertical regions that cannot 
be measured by the observing system, so this error is largest for atmospheric cases 
having the least variation in thermal structure (i.e., the subarctic winter case). One 
can also notice that the measurement noise is smallest for the tropical atmosphere, 
which possesses the highest signal-to-noise ratio. 

Table 5.5 summarizes the effect of all assumed errors on retrieved integrated O 3 
amounts for each of the various reference atmospheres considered. Approximate an- 


130 


nual means have been computed by averaging the July and January results for both 
the mid-latitude and subarctic cases. In addition, an overall (northern hemisphere) 
average has been formulated using the results from atmospheric types 1) through 
5). The accuracies have been estimated with knowledge of the retrieval sensitivity 
to forward model systematic uncertainties which, in general, become more sensitive 
(and thus less accurate) when the atmosphere has more moisture or higher tempera- 
tures. The corresponding precisions vary for these atmospheric cases in a less obvious 
manner due to the different behavior exhibited by its three components of random 
error. In terms of magnitudes, the author finds these results to be very encouraging; 
in fact, the tropospheric values are better than originally anticipated. In addition, 
these predictions compare quite favorably with capabilities of the limited number of 
instruments/techniques currently enabling similar observations. For example, mea- 
surements from the High-Resolution Infrared Sounder (HIRS), which is part of a set 
of three instruments defined as the TIROS Operational Vertical Sounder (TOVS), 
can be used to determine total ozone to within 5 — > 7% and 5 — ► 9% for preci- 
sion and accuracy, respectively (Ma et al., 1984; Planet et ah, 1984); the ozone 9.6 
/xm band is also sensed by HIRS, however, the retrieval process utilizes correlative 
data from the UV spectral region (i.e., Dobson spectrophotometer measurements). 
Data from the Total Ozone Mapping Spectrometer (TOMS) have been improving 
over its years of operation, mainly from improved characterization of instrumen- 
tal uncertainties and modifications to the data reduction software; accordingly, the 
current estimates of total ozone can be good to within ~ 2% and ~ 5% for preci- 
sion and absolute accuracy, respectively, depending on solar zenith angle (Bhartia, 
private communication, 1993). Although tropospheric ozone is not currently mea- 
sured directly from any satellite-based instrument, tropospheric ozone residuals can 
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be deduced from a technique utilizing independent measurements from two satellite 
instruments. Fishman et al. (1990) have used measurements from the Stratosphere 
Aerosol and Gas Experiment (SAGE) and TOMS instruments to infer tropospheric 
ozone integrated amounts, and report deriving the climatology of tropospheric ozone 
to within an accuracy of ~ 10% in the tropics and ~ 15% at mid-latitudes; note 
that these accuracies refer to annual mean values of tropospheric O 3 residuals, and 
that the monthly means may typically be accurate to within ~ 30% for this resid- 
ual technique. As a final comparison comment, it is worth mentioning again that 
the estimated tropospheric errors associated with this study would be smaller for 
atmospheres containing more tropospheric ozone; this would be the case for current 
typical atmospheres (since nominal amounts of tropospheric ozone have increased 
since formulation of the reference atmosphere database used in this study), as well 
as for atmospheric scenarios of particular interest possessing enhanced tropospheric 
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CHAPTER VI 


SUMMARY AND CONCLUSIONS 


The ability to monitor the global distribution of tropospheric ozone (0 3 ) would 
be beneficial for enhancement of basic scientific understanding as well as to possibly 
enable minimizing the occurrences of ill-health effects associated with exposure to 
elevated amounts of 0 3 in the lower atmosphere. Such a measurement capability can 
be achieved using a satellite-based device performing at high spectral resolution with 
a large signal-to-noise ratio (SNR), to enable observation in the pressure-broadened 
wings of strong O3 lines while minimizing the impact of associated undesirable signal 
contributions introduced from the terrestrial surface, interfering species, clouds, etc. 
The Fabry-Perot interferometer (FPI) provides high spectral resolution and high 
throughput capabilities that are fundamental to the success of this measurement 
task. An FPI optimized for tropospheric 0 3 measurements can simultaneously ob- 
serve a stratospheric component and thus the total 0 3 column abundance, through 
proper spectral design of the instrument channels. Concentrations of stratospheric 
O 3 are known to be decreasing; this may lead to an increase in biologically harm- 
ful solar ultraviolet (UV) radiation reaching the earth’s surface, which would cause 
detrimental health effects. 

A review of the necessary theory of Fabry-Perot interferometry has been pre- 
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sented in this thesis, including characteristics unique to ideal etalons, real etalons, 
and instrument configurations involving multiple etalons in series. In light of the 
intended design, emphasis has been placed on a double-etalon series configuration; 
accordingly, an expression for a double-etalon FPI instrumental transfer function has 
been formulated, including the associated broadening effects and incoherent reflec- 
tions between etalons. A conceptual double-etalon instrument design for measuring 
atmospheric ozone has been developed in this project. A double-etalon configura- 
tion in conjunction with an ultra narrow bandpass filter enable single-order passband 
transmission with high spectral resolution to achieve the desired spectral isolation. 
The central order full-width at half-maximum (FWHM) (or resolution) is ~ .054 
cm -1 for the etalons alone and an overall resolution of ~ .068 cm" 1 results after 
inclusion of aperture broadening introduced from the detection process. The new 
Circle to Line Interferometer Optical (CLIO) system has been incorporated into this 
design to efficiently handle the IR-FPI detection task while enabling multichannel 
operation. The CLIO system converts circular fringes of an FPI into a linear pattern 
which can be detected by a conventional linear array detector. Other significant in- 
strument design specifications include: a bandpass filter FWHM « 5.27 cm" 1 , defect 
finesse of « 37.5 for each etalon, plate coating reflectivities of .92 (corresponding to 
reflectivity finesses of % 37.7), an aperture finesse of « 37.7; and free spectral ranges 
of 1.52 cm" 1 and 5.46 cm" 1 for the high-resolution etalon and low-resolution etalon, 
respectively. Optical element and system transmittanc.es are as follows: a peak 
transmittance of .65 is assumed for the bandpass filter, the transmission through 
each defect-broadened etalon is ~ .73, .95 transmittance of the medium between 
etalons (which is so placed for reduction of inter-etalon reflections), and a value of 
.85 is assumed for the transmittance of the remaining optical elements not consid- 
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ered; this corresponds to an instrument peak transmittance of ~ .28. This design 
implies an integrated white light amount of .037 in the optical system, and an over- 
all effective system finesse of ~ 22.4. It is also worth noting that these instrument 
design specifications are all obtainable within capabilities of current technology. 

A method has been presented for calculating the upwelling radiance reaching 
a satellite sensor; the target gas ozone and all significant interfering species have 
been included. The water vapor continuum absorption provides the most significant 
impact of all interfering species (i.e., interfering spectral lines can be avoided in se- 
lection of the channel spectral locations and thus minimize their effects) in this 10 
/xm spectral region; consequently, modeling of its inclusion into the radiative trans- 
fer has been discussed in much detail. The desired instrument spectral resolution 
and approximate placement for the measurement band was determined by results 
of simulations assessing the radiance sensitivity to atmospheric ozone perturbations. 
Several criteria were then considered in determining the exact band placement as well 
as spectral locations for individual channels: ozone energy contribution functions (0 3 
ECFs) needed to demonstrate the existence of sufficient information content (while 
maxima were desired) in both the troposphere and stratosphere, effects of interfering 
species should be minimized, and a maximum sensitivity to perturbations in atmo- 
spheric ozone amount (especially for the troposphere) was desired. This phase of the 
research resulted in 24 channels being chosen within the 1054.23 to 1055.23 cm' 1 
spectral region. 

A signal-to-noise ratio (SNR) analysis has also been performed, which required 
definition of other design parameters. Assuming a 50 km nadir footprint along with 
orbital and instrument scan parameters similar to those of the Total Ozone Mapping 
Spectrometer (TOMS) instrument aboard Nimbus-7 (achieving continuous coverage 
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at the equator) yields a dwell time of ~ 220 ms for each target scene. A conservative 
telescope size of 9 inches in diameter has been assumed which then implies etalon 
plates of diameter ~ 13.3 cm. A design goal of an effective system f-number > 1.2 
has been reached. In order to achieve the maximum attainable SNR, the detector 
detectivity (D-star) values are desired as high as possible. Consequently, this design 
incorporates mercury cadmium telluride (HgCdTe) photoconductive (PC) detector 
elements within the linear array with cryogenic operation at 77 K, utilizing liquid 
nitrogen ( LN 2 ) cooling and 60° FOV cold stop shielding; this enabled a detectivity 
D * = 5.5 x 10 l0 cmH /W to be assumed for each detector element. The resulting 
SNR values are a function of the thermal characteristics of the atmosphere and 
surface below (in addition to the spatial distributions of radiatively active gases 
within the atmosphere); and for the U.S. Standard (1976) atmosphere correspond to 
a minimum-maximum range of approximately 133.2 to 289.9 for the various channels. 

The maximum likelihood method has been implemented for an ozone retrieval 
technique in this research. The ability to successfully retrieve integrated ozone 
amounts for both the troposphere and total column abundance, in the presence 
of measurement noise and null-space error, has been demonstrated for this con- 
ceptual measurement. It has also been shown that the retrieved profile shape can 
match the true profile shape quite well in some situations (i.e., when the character 
of the true profile is well represented by the a priori information). An error anal- 
ysis was performed to ascertain how good the retrieved ozone amounts would be 
in a statistical sense (encompassing many retrievals) while including the effects of 
measurement, null-space, and forward model errors. This was applied to six different 
reference atmospheres to acquire a feeling for the latitudinal and seasonal dependence 
of errors in retrieved ozone. Results show the proposed instrumentation to enable 
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a good measurement of actual ozone amounts and even a better determination of 
relative changes, with a range of accuracy to within 7.55 to 20.6% for integrated 
tropospheric amounts (and 1.99 to 4.02% for total 0 3 column abundance) and a 
corresponding range in precision to within 7.73 to 10.4% (and 3.30 to 3.95% for 
total 0 3 column abundance), for various atmospheric types considered. Note that 
these are predicted ozone errors based upon assumed uncertainties in the observing 
system (measurement process and retrieval technique) and consequently may vary 
slightly. Nevertheless, these predictions compare quite favorably with capabilities of 
instruments currently providing total ozone observations or techniques enabling de- 
duction of integrated tropospheric ozone amounts. It is also worth mentioning that 
the estimated tropospheric errors associated with this study would be smaller for 
atmospheres containing more tropospheric ozone; this would be the case for current 
typical atmospheres (since nominal amounts of tropospheric ozone have increased 
since formulation of the reference atmosphere database used in this study), as well 
as for atmospheric scenarios of particular interest possessing enhanced tropospheric 
ozone. Other significant findings include: the importance of including a water vapor 
continuum model, especially for a warm, moist tropical atmosphere; the significance 
of uncertainties in surface contribution characteristics (i.e., emissivity and tempera- 
ture); and great care must be taken in performing the radiometric signal calibration, 
as retrieved ozone amounts are quite sensitive to such errors. 

In conclusion, this thesis project has demonstrated that the Fabry-Perot inter- 
ferometer used in a double-etalon series configuration in conjunction with the new 
Circle to Line Interferometer Optical (CLIO) system is a feasible new approach for 
a satellite-based atmospheric sensor of tropospheric and total ozone. The FPI has 
a proven success record in previous space applications. While the CLIO subsystem 
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has only been demonstrated in the laboratory thus far, it is a desirable satellite in- 
strument component as it enables “spatial scanning” through multichannel operation 
without moving parts. In addition, all suggested design parameters are achievable 
based on current technological capabilities or expectations, making the proposed 
instrumentation even more attractive. 

As with any feasibility study, additional effort is still needed toward this partic- 
ular application. From a science perspective, it would be beneficial to determine the 
effect on retrieved ozone amounts of uncertainties in other forward model parame- 
ters, even for those model parameters which have not been included in the current 
study. For example, the impact of various cloud types and their spatial distributions 
when cloud cover is present within the instrument field-of-view (FOV), the sensi- 
tivity to variations in bandpass filter spectral positioning and shape, the sensitivity 
to variations in surface topography, and the impact of atmospheric aerosols and 
chlorofluorocarbons (CFC’s). In addition, it would be interesting to assess this in- 
strument’s capabilities when applied to a depleted polar stratospheric ozone scenario. 
Regarding engineering aspects, continued laboratory studies of the CLIO system are 
recommended, with particular emphasis toward IR applications. It would be helpful 
to explore pre-spaceflight etalon mounting techniques which minimize surface dis- 
tortions and thus avoid defect finesse degradation. In addition, it is suggested that 
ray-trace simulations be performed using optical design software to characterize ex- 
act energy distributions at the detector plane; this would enable a more accurate 
determination of detector element widths and minimum necessary acceptance angles 
for designing the detector cold stop FOV’s, both of which would improve the SNR 
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APPENDIX A 

THEORY OF FABRY-PEROT 
INTERFEROMETRY 


A.l Introduction 

Although early work done by George Airy with interference of radiation can 
be traced back to 1831, the birth of the Fabry-Perot interferometer is attributed 
to articles by Charles Fabry and Alfred Perot in 1897 and 1899 (Vaughan, 1989). 
The multiple beam Fabry-Perot interferometer was developed primarily as an im- 
provement on the two-beam Michelson interferometer. The two-beam interferometer 
produces fringes which are relatively diffuse whereas the multiple beam instrument 
yields much sharper fringes; this was a great advantage when the eye was being used 
as a detector (Chantry, 1982). The Fabry-Perot interferometer (FPI) is one of the 
most flexible and versatile of the multiple-beam interferential devices. Because of the 
large luminosity-resolving power product (LRP), it has been the instrument of choice 
for many high-resolution spectroscopic investigations, where the sources are usually 
intrinsically weak. It is commonly used as a high-resolution spectrometer, a resonant 
cavity in laser physics, a tuning selector in very high-resolution spectroscopy, etc., 
for a wide range of applications. Depending upon the application, various versions 
of the same instrument are implemented; accordingly, the corresponding name may 
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vary but is usually referred to as one of the following: Fabry-Perot interferometer, 
Fabry-Perot etalon, tunable etalon, confocal etalon, solid etalon, Piezoelectric Trans- 
ducer (PZT) scanning resonant cavity interferometer, fixed air gap etalon, or pressure 
tuned Fabry-Perot. However, this instrument is most commonly known as simply a 
Fabry-Perot interferometer. The FPI is an extremely powerful tool that has been 
employed in various forms for different applications over a wide range of wavelengths 
during the last several decades. As shown by Jacquinot (1954), the principle advan- 
tage of the FPI is that its throughput is much larger than that of a prism or grating 
instrument having the same spectral resolution. This appendix contains an overview 
of Fabry-Perot interferometry concepts essential to this research project, beginning 
with characteristics of an ideal FPI (one with perfectly flat, parallel plates) and then 
discussing the impact of accounting for realistic plate imperfections. 

A. 2 Characteristics of an Ideal FPI 

A Fabry-Perot interferometer is constructed with two partially transmitting mir- 
rors which are parallel to each other and may be flat (piano FPI) or curved (spherical 
or confocal FPI); this forms the FPI cavity. The “radiused” mirrors in a spherical 
FPI are separated by their radius of curvature for the special case of a confocal 
FPI. In a solid etalon, the two surfaces of one substrate are highly parallel and have 
partially transmitting coatings. The plane Fabry-Perot is comprised of two parallel 
flat, transparent plates coated with thin films of high reflectivity and low absorption, 
and is frequently called a Fabry-Perot etalon when the separation between the plates 
is fixed. An etalon has two mirrors, or plates as they are frequently called, which 
have partially transmitting coatings on their first surfaces and anti-reflection coat- 
ings on their second surfaces. The plates are arranged with their high reflectivity 
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Figure A.l: Schematic diagram of monochromatic ray propagation through a simple 
Fabry-Perot interferometer system. Light from the source is collected by 
the collimating lens. The energy undergoes successive reflections within 
the etalon cavity. The multiple waves are then superposed by the fringe 
forming lens which forms fringes in the focal plane. 

first surfaces facing each other. In the piano (or plane) FPI, the second surface is 
slightly wedged with respect to the first surface with a wedge angle « 30' (Vaughan, 
1989) to avoid forming additional (resonant) cavities due to reflections from the low 
reflectivity surfaces. The plane Fabry-Perot is the simplest and most versatile of the 
FPI instruments. Etalons of this kind will be the type of FPI’s simulated in this 
research project. 

The space between the coated surfaces of the etalon forms a cavity which is res- 
onant at certain wavelengths dictated by the optical properties (i.e., thickness t and 
refractive index fi) of the gap as well as the incidence angle 0 of radiation within 
the cavity. Due to the low transmission coating, only a small portion of the incident 
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Reflectivity 


Figure A. 2: Ideal FPI throughput versus reflectivity for several values of plate ab- 
sorption. 
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Figure A. 3: Transmission through ideal etalon for different plate reflectivities. 

radiation enters the cavity where it then undergoes multiple reflections. A simpli- 
fied schematic of ray propagation through an FPI system is shown in Figure A.l. 
The amplitude and phase of radiation leaving the cavity depend upon wavelength. 
Constructive interference produces transmission maxima at resonant wavelengths; 
conversely, at wavelengths in between these maxima, destructive interference causes 
most of the incident energy to be reflected thus producing transmission minima. The 
general equation for transmission through an ideal FPI (assuming parallel, perfectly 
flat plates with no surface defects; thus having broadening caused only from reflec- 
tivity of the surface coatings) is given by the Airy formula for transmission (Born 
and Wolf, 1975): 


AW) = r a 


(1 ~ A ) 2 

1 — 2 R cos 4- R 2 ’ 


(A.l) 


where r a is a transmission (or maximum throughput) coefficient that includes the 
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effects of absorption and/or scattering by the reflective coatings, R is the plate re- 
flectivity, 4* = 47r/xfi/cos0 + 2 <j>, ji is the index of refraction of the medium between 
the plates, t is the spacing (or gap) between the plates, (f> is the phase change upon 
internal reflection, v is the wavenumber of incident radiation, and 6 is the angle of 
incidence of the incident radiation. The transmission coefficient r a corresponds to 
the throughput of the ideal FPI at resonance (or its peak transmission) and for an 
etalon of plate reflectivity R is given by 

r - = ( 1 -rrfl) 2 ’ < A - 2 > 

where A corresponds to losses due to scatter and absorption for a single plate. The 
quantity (equation A.l normalized to peak transmission) is defined as the Airy 
shape function or simply “Airy function,” however, «4('I') alone is commonly referred 
to as the Airy function as well. Figure A. 2 shows the dependence of throughput on 
reflectivity for several values of A; A is typically < .2%/surface for modern coatings 
(Burleigh). It can be seen from equation A.l that the phase changes upon reflection 
have no effect on the shape of the Airy transmission function distribution, but serve 
merely to offset the fringes; and for higher orders of interference (> 100) they can be 
ignored for most practical purposes (Vaughan, 1989). Since the phase contributions 
can be safely ignored for this application, equation A.l implies that the FPI has 
a “comb” of transmission resonances with the peaks in transmission occurring in 
narrow spectral regions at each of a series of wavenumbers given by: 

= 47r/xf// cos 0 = 27TU, (A. 3) 

where n is an integer known as the order of interference. It is obvious from equa- 
tions A.l and A. 3 that if the Airy function is plotted as a function of any one of the 
quantities n, t, v, or cos 0, one would obtain a regular series of peaks occurring wher- 
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ever 2/i/z/cos0 = n, the integer order. This is illustrated in Figure A. 3 for several 
values of the plate reflectivity. 

Taking the total derivative of equation A. 3, an expression for the local change in 
order An can be derived giving 


Av A t A/z 

An = n — n 0 = 2/z 0 Mo cos0 o ( h - — I 0oA9), 

Vo to ^ 0 


(A.4) 


where the zero subscript denotes some nearby reference state. For a fixed gap etalon 
with a fixed (or small) angle of incidence (i.e., for a particular “channel”), equa- 
tion A.4 reduces to 

An = 2/x 0 t()A;/ cos 6 0 , (A. 5) 


where Ai/ = u — uq. 

In discussing an axial-fringe FPI, the Airy function is commonly considered as 
a function of v only, with the other parameters fixed. In this case, the spectral 
spacing between consecutive passbands at a fixed angle of incidence 6 , the inter- 
order distance, is known as the free spectral range vj sr (FSR). The free spectral 
range can be obtained from equation A. 5 by solving for the wavenumber difference 
Ai / induced from consecutive integer values of the order n; it is given by: 


Vjsr = 


(A. 6) 


2 lit cos 6 

However, for most axial-fringe instrument applications cos 6 is approximately equal 
to 1 and the FSR may be accurately represented by: 



(A.7) 


this is commonly referred to as the “nominal” free spectral range. 

The half-intensity width 6p, which is sometimes loosely called the resolving limit, 
is the full width at half-maximum (FWHM) of a transmission peak in the FPI instru- 


Normalized Transmission 
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Change From Central Order (An) 


Figure A. 4: Transmission through an ideal etalon with the corresponding FWHM 
(Sjr) and FSR (I'/sr) depicted. 
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mental function. For an ideal etalon having perfectly flat plates, with only reflectivity 
broadening existing, this is commonly expressed (i.e., Vaughan, 1989) as: 

The FWHM and FSR of an Airy transmission function are illustrated in Figure A. 4. 
The theoretical spectral resolution, or resolving power, of an FPI is defined by 

ftsf. (A-9) 

OF 

The finesse N is the key measure of the interferometer’s spectral resolving capa- 
bility. The finesse is defined as the ratio of the FSR near the resonance frequency to 
the FWHM of an etalon transmission peak, and is given by 

N= V -¥~. (A. 10) 

of 

Thus higher finesse implies higher resolution, for a fixed FSR. The finesse can also 
be thought of as the number of passband widths over which the FPI can be scanned 
(or tuned) without overlapping of orders. For an ideal etalon, N is a function of only 
the reflectance of the plate coatings and is referred to as the “reflectivity finesse,” 
Nr, which from equations A. 8 and A. 10 can be directly represented by 

_7tVR_ (A. 11) 

R (1 -R) v ’ 

The dependence of reflectivity finesse on plate reflectance is illustrated in Figure A. 5. 

Another useful parameter is the contrast C which is defined as the ratio of the 
maximum to minimum FPI transmittance. For an ideal etalon with perfectly flat 
plates, the contrast is simply the ratio of maximum to minimum Airy function trans- 
mission, and can be easily shown to be given by 

c _ U + fl ) 2 

(1 — R) 2 ' 


(A.12) 
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Figure A. 5: Reflectivity finesse versus reflectivity for an ideal etalon. 

A. 3 The Real Etalon 

The theory presented in the last section for an ideal etalon assumes parallel, 
perfectly flat plates. In reality, the Fabry-Perot plates are neither flat nor perfectly 
parallel; and any deviation in either one (or both) of these quantities will broaden the 
interferometer passband, requiring modification to some of the equations developed 
in Section A. 2. 

Departures from flatness of the interferometer plates, errors in their parallel align- 
ment, and nonuniformities in the reflective coatings cause the effective plate separa- 
tion t to vary across the surface of the etalon; thus making the optical pathlength 
vary over the aperture of the interferometer. Consequently, the Airy function (equa- 
tion A.l) does not coincide for all surface elements within the aperture, and the 
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resulting fringes are broadened. Broader fringes imply that the same energy must 
be distributed over a larger area in the focal plane, thus reducing the peak intensity. 
Therefore, the peak transmission or throughput of a “real” etalon is less than that 
given for an “ideal” etalon in equation A. 2. 

The transmission profile of a broadened fringe is necessary in order to simulate 
a real etalon. This can be determined with knowledge of two functions: the Airy 
function (equation A.l) representing the fringe profile for an ideal etalon and 

the surface defect distribution function D(y>). In these functions, ^ is the “phase 
advance” of A defined by equation A.l, and <p corresponds to the phase decoupling 
or a perturbation in 'k. For random surface defects, D(i^) can be considered as a 
probability density function describing the distribution of phase advance ^ errors <^>. 
For surface defects of a nonrandom nature, D(<y?)d< / ? corresponds to the proportion 
of the FPI aperture for which the phase advance differs from some value of 'k by the 
amount p due to plate defects. Therefore D(i^) has the attributes of a probability 
density function with the area under the defect function curve being normalized to 
unity, denoted mathematically by 


Accordingly, the mean //£> and variance ajj of the defect function can be defined, 
respectively, by 


r°° 

/ D((p)d<p = 1. 

J — oo 


(A.13) 



(A- 14) 


/ OO 

D(p)(ip - p D )d<p. 

-OO 


(A. 15) 


Following an approach similar to Sloggett (1984) and Wilksch (1985), the fractional 
contribution to the total transmission of the etalon due to those parts of the FPI 
aperture having phase defect p can be given by ^(ty — <p)D((p)d<p. The total etalon 
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transmission at 'I' can be obtained upon integration, yielding 

/ oo 

A(* - <p)D((p)dip; (A. 16) 

-OO 

E('I') can be called the etalon transmission function (ETF). Equation A. 16 can be 
more conveniently represented by 

£('I') = .4('I') *£>(*), (A. 17) 

where * denotes the mathematical process of convolution (see, for example, Churchill, 
1972). 

The Airy function .A('I') describes symmetrical fringes. If the mean (or central 
maximum) occurs at ha-, then from a standard property of convolutions the mean of 
the defect broadened fringes he can be expressed (Sloggett, 1984) as he = /M + HD • 
Therefore, the mean of the defect function hd is the amount of displacement of 
the Airy function peaks found in the broadened fringe profile E('k), due to the 
plate defects. In addition, the peaks (or fringes) described by E('k) will retain their 
symmetry as long as the defect distribution function is symmetrical about its mean. 

Analytical expressions for defect functions representing many common plate im- 
perfections have had many appearances in the literature [e.g., Del Piano and Que- 
sada (1965), Hernandez (1966, 1986), Sloggett (1984), Wilksch (1985), and Vaughan 
(1989)]. In general, the defect functions are developed by assuming that all imper- 
fections are on one of the reflecting plates which is positioned parallel to the surface 
of a flawless second plate. Figure A. 6 illustrates some of the more typically observed 
plate imperfections along with their corresponding defect distribution functions most 
commonly used to represent them. Figure A. 6(a) is a situation of flat plates with 
random defects of small scale and can correspond to plates with poor coatings. The 
Gaussian distribution function can be used to represent the defect function as this 
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Plate Imperfection 


Defect Function 
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(b) 


(c) 
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Figure A. 6: Common etalon plate imperfections and their corresponding defect dis- 
tribution functions plotted as a function of phase error normalized to 
the standard deviation ctd for each distribution. The spacing between 
plates (or gap) has a mean of t and a range t ± A (except points greater 
than one crp from t in the Gaussian distribution), (a) shows a random 
surface irregularity which can be represented by the Gaussian distribu- 
tion function, (b) is a case where the plates have a spherical defect and 
a rectangular function can be used for the defect function, (c) illustrates 
a pair of flat plates slightly out of alignment, for which a semielliptical 
function can represent. 
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Distribution Type 

Analytical Expression 

DM 

Standard Deviation 

Gaussian 

{V2tt cr) 1 exp[- 

cr 

Rectangular 

(2\/3cr ) -1 

tymax 

VT 

Semielliptical 

1 /i _ 

7T (7 V 1 4<7 2 

Vmax 

2 


Table A.l: Analytical expressions used to represents defect distribution functions for 
some common etalon plate defects. Note that all distribution functions 
except for the Gaussian vanish for |y> — hd\ > max , where v> max is the 
value of 9 ? corresponding to the limit of defect function nonzero extent. 

case is that of a small scale random surface irregularity. Figure A. 6 (b) shows plates 
with a spherical defect. This type of defect might arise from plates usually being of 
slightly poorer quality toward the edges compared with the center; in addition, the 
resultant strain from faulty mounting of the plates can induce this type of curvature 
effect. For circular plates (or a circular aperture), a rectangular function can be 
used to represent the corresponding defect function. Figure A. 6 (c) illustrates a pair 
of flat plates slightly out of alignment. A semielliptical function can represent its 
defect distribution for circular plates; and a hemispherical distribution may also be 
used as a good approximation (Vaughan, 1989). Table A.l summarizes the above de- 
scribed surface defect distribution functions which are illustrated in Figure A. 6 . The 
expressions given for the standard deviations can be obtained through evaluation of 
equation A. 15. 

The magnitude of a plate defect is commonly expressed as a fraction of wave- 
length. This usually corresponds to a maximum excursion from a perfectly plane, 
parallel surface, and is denoted by 2 A in Figure A. 6 . For situations with random 
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Defect Type 

Defect Magnitude 

Defect Finesse 

Distribution Function 
Spectral Extent 

Random surface 
imperfections 

(«$)* = A 

Nd g = 7 ? 

T — U l ST 

D 2\/21n2 N Dg 

Spherical curvature 

2A * = A 

Nd s = 

U - "I” 

'pmax 2 N Dq 

Parallelism error 

= 


Vfsr 

^ max ~ 75 Nd~ p 


Table A. 2: Plate defect characteristics. 

surface irregularities, a root-mean-square (rms) deviation ( 6^)2 is normally given to 
denote microscopic flatness imperfections. A spherical curvature of sagitta is often 
used as a measure of macroscopic deviations from flatness. The broadening effect 
of defects on the FPI transmission function can be quantified by a defect finesse 
for each defect. And the net finesse is obtained by treating the component finesses 
as parallel impedances. An overall finesse for the above described defects occurring 
simultaneously can be given by 


1 


+ 


+ 


1 


(A. 18) 


A j d 2 Nd g 2 Nd s 2 Nd p 2 

The spectral extent of the defect distribution function D(y>) can be characterized by 
a standard deviation gd and/or y> ma x, where < p max is the value of tp corresponding to 
the limit of defect function nonzero extent. 

The broadening effect from spherical defects and surface irregularities has been 
extensively covered in the literature; for example, a thorough discussion of these 
types of broadening has been given by Hernandez (1966, 1986) and by Hays and 
Roble (1971). A treatment of misalignment broadening due to a lack of parallelism of 
the etalon plates has been done by Gupta and Prasad (1991). Atherton et al. (1981) 
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cover characterization of these defects including their corresponding defect finesses. 
These plate defect characteristics are summarized in Table A. 2. The expression for 
standard deviation of the Gaussian defect distribution can be easily arrived upon 
by combining the definition of standard deviation for a Gaussian (a = where 

6h =half-width at half-maximum (HWHM)) with the HWHM expression for a defect 
(6h d = where the defect finesse Nd = Np a for this case) given by Hays 

and Roble (1971). The formulation for <p max for a spherical defect can be obtained 
by noting that tp max = Sh d for a rectangular distribution. And, finally, the <p max 
relationship given for the parallelism error is obvious from the development presented 
in Gupta and Prasad (1991). 

A particular defect distribution function may apply to several very different types 
of surface defects. For example, Sloggett (1984) discusses how a rectangular distri- 
bution may also be used to characterize uniformly distributed random defects, a 
parabolic defect, and even a parallelism error in plates used with a square aperture. 
In addition, other types of plate defects may exist which are best described using 
distribution functions other than those given above. For example, a surface defect of 
sinusoidal form might result from certain polishing techniques (Vaughan, 1989); its 
corresponding defect distribution function would be an inverted semiellipse having 
peaks at the extremes of phase difference. 

Equations A. 16 and A. 17 show how the defect distribution function D(y>) can be 
used to arrive upon the etalon transmission function E('I'). In particular, the integral 
in equation A. 16 may be solved analytically for several forms of the defect function 
(Hernandez, 1966, 1986; Hays and Roble, 1971), or it may simply be evaluated nu- 
merically. Alternatively, the Airy function can be numerically convolved with the 
defect distribution function (equation A. 17) using the discrete Fast Fourier Trans- 
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Figure A. 7: Etalon transmission characteristics in the presence of plate defects. Plots 
on the left side of the Figure show passband transmission variation near 
the central order normalized to the peak transmission for several values 
of plate defect finesse; whereas the plots on the right depict the peak 
transmission dependence on plate defect finesse. The solid curves repre- 
sent plates without defects, (a) represents a Gaussian defect distribution 
function for a random surface imperfection, (b) corresponds to a rectan- 
gular defect distribution function for a spherical plate defect. A semiel- 
liptical defect function was used in (c) to represent a plate parallelism 
error. And (d) shows the combined effect of the defects represented in 
(a) through (c) occurring simultaneously. 
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form (FFT) (Churchill, 1972). This later approach has been implemented in this 
research upon its validation through comparisons among the different techniques. 
The Convolution Theorem (Churchill, 1972; Press et al., 1990) may be interpreted 
as saying that the Fourier transform of a convolution is equal to the product of the 
individual Fourier transforms. Accordingly, equation A. 17 may be expressed as 

E(V) = A(V) * D(V) = FFT- 1 {FFT{A('H)}FFT{D{'i!)}}. (A. 19) 

For the realistic case of several plate imperfections existing simultaneously, the defect 
distribution function can be formulated by convolving all of the individual defect 
functions together. This can be expressed as 

D{V) = D a {*) * D S {V) * DpiV) (A. 20) 

for the three main defects discussed earlier. Figure A. 7 illustrates the effects of 
plate imperfections on the characteristics of the etalon transmission function. The 
reflectivity and absorption/scattering of the plates are assumed to be .94 (Nr ~ 50) 
and .002, respectively. The solid curves represent plates without defects; for ex- 
ample, the transmission function with reflectivity broadening only (Airy function, 
equation A.l), and peak transmission reduced only from absorption/scattering of 
the plates(equation A. 2). Figures A.7(a)-(c) show how Gaussian, rectangular, and 
semielliptical defect distribution functions, respectively, affect the etalon transmis- 
sion function when acting alone for different values of defect finesse. Whereas, Fig- 
ure A. 7(d) depicts the combined effect of the first three defects occurring simultane- 
ously, with equal individual defect finesses of same value as the corresponding curves 
in (a) through (c) of the Figure (and thus yielding an overall defect finesse given 
by equation A. 18). As the plate quality gets better (Np — ► Np Nr ), the etalon 
transfer function approaches the Airy function; conversely, as the plate quality gets 
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poorer ( No —*■ No <C Nr), the etalon transfer function gets broadened into a shape 
resembling the defect function itself. 

In addition to plate imperfections, an aperture in the optical system will also 
produce broadening of the FPI transfer function. The radiation passing through 
an etalon will cover some finite solid angle, which is usually defined by an aperture 
placed in the fringe plane. The most common configuration involves a fixed, circular 
aperture centered on the optical axis. And, as with plate defects, the effect will be 
to broaden the FPI instrumental transmission function. This can be represented 
mathematically by 

7W = £W*ZW, (A. 21) 


where //(ty) is the aperture function, £’('!') is the etalon transmission function as 
given by equation A. 17, and T’('I') corresponds to the FPI instrumental transmission 
function; it expresses the transmission of the interferometer with imperfect etalon 
plates and a finite sized aperture. As with the function used to represent a spherical 
defect of etalon plates, the aperture function can also be expressed as a rectangular 
function (Hernandez, 1986; Hays and Roble, 1971). However, the half-width at half- 
maximum (HWHM) of the aperture function is related to the optical geometry of 
the FPI; and for a circular aperture in the fringe plane it can be given by 

ud L 2 


8Hl (4/) 2 ’ 


(A. 22) 


where v is the wavenumber of incident radiation, d i is the diameter of the circular 
aperture, and f is the focal length of the fringe forming lens. As mentioned by Hays 
and Roble (1971), a true finesse cannot be assigned to aperture function broadening 
since (unlike the Airy and defect functions) its HWHM 6r l is not a function of the 
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plate separation t. However, it is still useful to define a pseudo aperture finesse by 


N l 


Vjjr_ 

2 <V 


(A. 23) 


which for the circular aperture system can be given by Nl = ~ • The 

effective resolution (or full-width at half-maximum (FWHM)) of the FPI is a function 
of the reflectivity, defect, and aperture finesses; this can be given approximately by 


f>F B = 


Vfsr 

~Ne' 


where Ne is an “effective finesse” that can be defined by 

1 _ i 1 1 

N e 2 “ N r 2 + N d 2 + Nl 2 ' 


(A. 24) 


(A. 25) 


The luminosity-resolution product (LRP) (Jacquinot, 1954) is usually optimum when 
Nr ~ Ne « Nl (Atherton et al., 1981). This is usually accomplished by first obtain- 
ing plates with a minimum of imperfections (maximizing No), then applying a reflec- 
tivity coating which equates Nr to Nd, and then implementing an aperture/detector 
system such that Nl ~ No- However, depending upon the application, it may be 
desirable to alter this finesse selection procedure and allow Nr, Nd, and Nl to have 
different values. For example, if the expected signal level is of concern and if the 
corresponding resolution is higher than necessary, then Nr and Nl can be made less 
than N d • Conversely, to increase resolution at the expense of lowering the signal, one 
can specify Nr and Nl to be greater than Nd- However, it should be noted that as 
Nr ~ N d — < ► Nr > Nd there would be significantly more reduction in throughput 
than would be gained in overall instrument finesse. 
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APPENDIX B 

ERROR ANALYSIS: SENSITIVITIES TO 
SYSTEMATIC PERTURBATIONS 


This appendix contains supplemental material for the error analysis discussed 
within Section 5.6. The following tables summarize the potential impact on re- 
trieved integrated ozone amounts from systematic offsets of certain forward model 
parameters for reference atmospheres used in this study. 
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Results of Error Analys 

sis: 2) mid-latitude summer (45N July) 

Effect on retrieved integrated O 3 amounts 

Perturbation 

( 

Troposphere 

Accuracy (%) 
systematic errors) 
Total 

H 2 O continuum model 

without H 2 O continuum 
Cs & Cn +10% 

24.43 

-2.318 

3.856 

-0.3659 

surface emissivity 

€ s -1% 

-9.704 

-1.495 

surface temperature 

T s + IK 

17.30 

2.666 

temperature profile 

T+1K all levels 

4.838 

-3.607 

interfering species 

CO 2 +l% all levels 
C# 4 + 10% all levels 
# 20 + 10 % afl levels 
JV 2 O+30% all levels 
SO 2 +30% all levels 

-0.02581 

0.0001292 

-4.049 

-0.0005321 

-0.0001197 

-0.01027 

-0.0001588 

-0.7317 

6.873e-05 

1.531e-05 


Table B.l: The effect of certain systematic perturbations on retrieved integrated 
ozone amounts for a mid-latitude summer atmosphere. 
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Results of Error Analys 

sis: 3 ) mid-latitude winter ( 45 N Jan) 

Effect on retrieved integrated 0 3 amounts 

Perturbation 

Accuracy (%) 
(systematic errors) 


Troposphere 

Total 

H 2 O continuum model 



without H 2 O continuum 

3.666 

0.1076 

Cs k Cn +10% 

-0.3633 

-0.01062 

surface emissivity 



€ s - 1 % 

-9.641 

-0.2660 

surface temperature 



t,+ik 

20.08 

0.5540 

temperature profile 



T+1K all levels 

-3.995 

-0.8588 

interfering species 



CO 2+IV0 all levels 

-0.01034 

-0.003279 

CH 4 + 10% all levels 

0.0002633 

-6.712e-05 

H 20 +lO%o all levels 

-0.5218 

-0.02885 

N 2 O+ 30% all levels 

-0.0002198 

5.256e-05 

502+30% all levels 

2.605e-05 

2.116e-06 


Table B.2: The effect of certain systematic perturbations on retrieved integrated 
ozone amounts for a mid-latitude winter atmosphere. 
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Results of Error Analys 

>is: 4) subarctic summer (60N July) 
Effect on retrieved integrated O 3 amounts 

Perturbation 

(s; 

Troposphere 

Accuracy (%) 
ystematic errors) 
Total 

H 2 O continuum model 

without H 2 O continuum 
Cs & Cn +10% 

22.04 

-2.126 

1.908 

-0.1841 

surface emissivity 

€ s -1% 

-11.28 

-0.9499 

surface temperature 

T S + 1 K 

21.10 

1.777 

temperature profile 

T+ 1 K all levels 

0.8645 

-3.828 

interfering species 
CO 2 + 1 % all levels 
C //4 4 T 0 % all levels 
H 2 O+ 10 % all levels 
■^20+30% all levels 
502+30% all levels 

-0.02320 

0.0004479 

-3.604 

-0.0001661 

-1.263e-05 

-0.007093 

-0.0001548 

-0.3845 

4.002e-05 

1.959e-05 


Table B.3: The effect of certain systematic perturbations on retrieved integrated 
ozone amounts for a subarctic summer atmosphere. 
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Results of Error Analys 

sis: 5) subarctic winter (60N Jan) 
Effect on retrieved integrated O3 amounts 

Perturbation 

Accuracy (%) 
(systematic errors) 


Troposphere 

Total 

H 2 O continuum model 



without H 2 O continuum 

0.4382 

0.0009927 

Cs & Cn -+"10% 

-0.04363 

-0.0001334 

surface emissivity 



c, -1% 

-6.913 

0.02194 

surface temperature 



T S +1K 

16.14 

-0.05098 

temperature profile 



T+1K all levels 

-5.797 

3.596 

interfering species 



CO 2 + 1 % all levels 

-0.001987 

-0.001434 

CH 4 +10% all levels 

0.0001559 

-5.254e-05 

H 20 + 10 % all levels 

-0.04769 

-0.002093 

^0+30% all levels 

2.226e-05 

-2.203e-05 

502+30% all levels 

1.546e-06 

-7.845e-06 


Table B.4: The effect of certain systematic perturbations on retrieved integrated 
ozone amounts for a subarctic winter atmosphere. 
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Results of Error Analyj 

sis: 6) U.S. standard (1976) 

Effect on retrieved integrated O 3 amounts 

Perturbation 

(s 

Troposphere 

Accuracy (%) 
ystematic errors) 
Total 

H 2 0 continuum model 

without H 2 0 continuum 
Cs k Cn +10% 

15.31 

-1.503 

1.470 

-0.1441 

surface emissivity 

e s -1% 

-13.16 

-1.241 

surface temperature 

t s +ik 

24.44 

2.305 

temperature profile 

T+1K all levels 

-2.057 

-4.307 

interfering species 

C0 2 +1% all levels 
CH 4 +\Q% all levels 
# 20 + 10 % all levels 
N 2 O+ 30% all levels 
SO 2 +30% all levels 

-0.02481 

0.0003253 

-2.450 

0.0001823 

0.0001796 

-0.008861 

-0.0001757 

-0.2948 

2.717e-05 

-6.827e-06 


Table B.5: The effect of certain systematic perturbations on retrieved integrated 
ozone amounts for the U.S. standard (1976) atmosphere. 
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APPENDIX C 

ERROR ANALYSIS: OVERALL ERROR 

BUDGETS 


This appendix contains supplemental material for the error analysis discussed 
within Section 5.6. The following tables summarize estimates of overall errors in 
retrieved integrated ozone amounts based upon the assumed systematic and random 
uncertainties for reference atmospheres used in this study. 
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Results of Error Analysis: 2) mid-latitude summer 

(45N July) 



Effect on retrieved integrated O 3 amounts 


Accuracy (%) 

Precision (%) 

Perturbation 

(systematic 

errors) 

(random errors) 


Troposphere 

Total 

Troposphere 

Total 

H 2 O continuum model 





Cs +10% 

-1.675 

-0.2638 



Cn - 1 - 10 % 

-0.6466 

-0.1027 



surface emissivity 





€ s ± 1 % 



2.215 

0.4450 

signal calibration 





all channels - 1 - 1 % 

11.95 

0.2213 



surface temperature 





T s ± IK 



3.948 

0.7932 

temperature profile 





T ± 2 k one level at a time 



1.238 

0.9294 

interfering species 





H 2 O ± 30% one level at a time 



0.9955 

0.2027 

CO 2 ± 5% one level at a time 



0.01317 

0.003535 

spectral line parameters 





H 2 O halfwidths +5% 

-0.01318 

-0.02858 



O 3 halfwidths +5% 

-8.885 

-3.000 



H 2 O line strengths + 10 % 

0.03848 

-0.08861 



0 3 line strengths + 1 % 

-1.349 

-0.8957 



Error impact (rss) 





forward model uncertainties 

15.12 

3.278 

4.798 

1.316 

measurement noise 



2.328 

1.218 

null-space error 



5.582 

2.781 

total (rss) errors 

15.12 

3.278 

7.730 

3.311 


Table C.l: Overall errors in retrieved integrated ozone amounts corresponding to ap- 
plication of the proposed instrumentation to a mid-latitude summer at- 
mosphere. 
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Results of Error Analysis: 3^ 

mid-latitude winter (45N Jan) 

Effect on retrieved integrated O3 amounts 


Accuracy (%) 

Precision (%) 

Perturbation 

(systematic errors) 

(random 

errors) 


Troposphere 

Total 

Troposphere 

Total 

H 2 0 continuum model 





Cs +10% 

-0.1671 

-0.004771 



Cn +10% 

-0.1966 

-0.005808 



surface emissivity 





( s ± 1% 



2.213 

0.5462 

signal calibration 





all channels +1% 

8.445 

-0.1085 



surface temperature 





T, ± IK 



4.609 

1.137 

temperature profile 





T ± 2k one level at a time 



0.7924 

0.5695 

interfering species 





H 2 0 ± 30% one level at a time 



0.1337 

0.03042 

C0 2 ± 5% one level at a time 



0.004020 

0.0009045 

spectral line parameters 





H 2 0 halfwidths +5% 

0.0006142 

-0.002767 



O3 halfwidths +5% 

-8.273 

-1.716 



H 2 0 line strengths +10% 

0.02382 

-0.01300 



0 3 line strengths +1% 

-1.189 

-0.7040 



Error impact (rss) 





forward model uncertainties 

11.94 

1.987 

5.175 

1.385 

measurement noise 



2.517 

1.284 

null-space error 



6.426 

3.154 

total (rss) errors 

11.94 

1.987 

8.627 

3.681 


Table C.2: Overall errors in retrieved integrated ozone amounts corresponding to ap- 
plication of the proposed instrumentation to a mid-latitude winter atmo- 
sphere. 
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Results of Error Analysis: 4] 

1 subarctic summer (60N July) 

Effect on retrieved integrated O 3 amounts 


Accuracy (%) 

Precision (%) 

Perturbation 

(systematic 

errors) 

(random errors) 


Troposphere 

Total 

Troposphere 

Total 

H 2 O continuum model 





Cs + 10 % 

-1.436 

-0.1239 



Cn +10% 

-0.6929 

-0.06028 



surface emissivity 





€3 ± 1 % 



2.363 

0.4212 

signal calibration 





all channels + 1 % 

11.98 

-0.5733 



surface temperature 





T s ± IK 



4.420 

0.7877 

temperature profile 





T ± 2 k one level at a time 



0.9182 

0.8352 

interfering species 





H 2 O ± 30% one level at a time 



0.8793 

0.1518 

CO 2 ± 5 % one level at a time 



0.008623 

0.002129 

spectral line parameters 





H 2 O halfwidths +5% 

-4.776e-05 

-0.01980 



O 3 halfwidths 4-5% 

-9.227 

-2.630 



H 2 O line strengths 4 - 10 % 

0.07047 

-0.06774 



0 3 line strengths 4 - 1 % 

-1.260 

-0.8150 



Error impact (rss) 





forward model uncertainties 

15.31 

2.932 

5.171 

1.232 

measurement noise 



2.360 

1.225 

null-space error 



5.939 

2.822 

total (rss) errors 

15.31 

2.932 

8.224 

3.317 


Table C. 3 : Overall errors in retrieved integrated ozone amounts corresponding to ap- 
plication of the proposed instrumentation to a subarctic summer atmo- 
sphere. 
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Results of Error Analysis: 5] 

1 subarctic winter (60N Jan) 

Effect on retrieved integrated O3 amounts 


Accuracy (%) 

Precision (%) 

Perturbation 

(systematic errors) 

(random 

errors) 


Troposphere 

Total 

Troposphere 

Total 

H 2 0 continuum model 





Cs +10% 

-0.008697 

-3.080e-05 



Cn +10% 

-0.03506 

-6.291e-05 



surface emissivity 





0) 

H- 

25 



2.083 

0.5924 

signal calibration 





all channels +1% 

5.007 

0.9888 



surface temperature 





T, ± IK 



4.863 

1.383 

temperature profile 





T ± 2k one level at a time 



1.026 

0.8166 

interfering species 





H 2 0 ± 30% one level at a time 



0.02406 

0.006576 

CO 2 ± 5% one level at a time 



0.001325 

0.0004652 

spectral line parameters 





H 2 O halfwidths +5% 

0.0001325 

-0.0003161 



O3 halfwidths +5% 

-5.522 

-1.763 



H 2 O line strengths +10% 

0.004202 

-0.001945 



O3 line strengths +1% 

-0.8424 

-0.6742 



Error impact (rss) 





forward model uncertainties 

7.549 

2.235 

5.389 

1.712 

measurement noise 



2.893 

1.290 

null-space error 



8.313 

3.303 

total (rss) errors 

7.549 

2.235 

10.32 

3.945 


Table C.4: Overall errors in retrieved integrated ozone amounts corresponding to ap- 
plication of the proposed instrumentation to a subarctic winter atmo- 
sphere. 
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Results of Error Analysis: 6 ) U.S. standard (1976) 




Effect on retrieved integrated O 3 amounts 


Accuracy (%) 

Precision (%) 

Perturbation 

(systematic 

errors) 

(random errors) 


Troposphere 

Total 

Troposphere 

Total 

H 2 O continuum model 





Cs +10% 

-0.8781 

-0.08415 



Cn +10% 

-0.6261 

-0.06031 



surface emissivity 





€5 ± 1 % 



2.894 

0.4772 

signal calibration 





all channels + 1 % 

12.76 

-0.2557 



surface temperature 





T s ± IK 



5.376 

0.8864 

temperature profile 





T ± 2 k one level at a time 



1.116 

0.9058 

interfering species 





H 2 O ± 30% one level at a time 



0.5966 

0.09579 

CO 2 ± 5 % one level at a time 



0.01057 

0.002926 

spectral line parameters 





H 2 O halfwidths +5% 

-0.01611 

-0.01684 



O 3 halfwidths -f 5% 

-10.03 

-2.503 



H 2 O line strengths + 10 % 

0.02241 

-0.05835 



O 3 line strengths + 1 % 

-1.358 

-0.8295 



Error impact (rss) 





forward model uncertainties 

16.38 

2.779 

6.236 

1.358 

measurement noise 



2.793 

1.272 

null-space error 



6.808 

2.868 

total (rss) errors 

16.38 

2.779 

9.656 

3.421 


Table C. 5 : Overall errors in retrieved integrated ozone amounts corresponding to ap- 
plication of the proposed instrumentation to the U.S. standard (1976) 
atmosphere. 
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